泥石流固液分相流速计算方法研究

泥石流固液分相流速是泥石流对岸坡、防治结构冲击、磨损机理的核心问题．将泥石流体简化为具有相同粒径的固相和具有相同力学性质的液相,基于泥石流体为沿流动方向的一维两相流体,运用两相流理论建立了泥石流固液分相流速控制方程．构建了泥石流平均压力、彻体力及平均表面力的计算方法,尤其通过浆体的Binhanm体流变方程、Bagnold颗粒相互作用试验成果建立了控制体平均表面力计算方法;建立了固液两相流速比例系数,以及理论固相流速与实际流速的比例系数．据此求解控制方程得到了固液分相流速计算方法,该方法既可同时适用于粘性泥石流和稀性泥石流,也可在泥石流爆发以后通过现场采集沉积物分析反求泥石流爆发期间的分相流速．工程实例分析显示,该方法计算结果与实测结果吻合较好．     

A comparison of numerical models for one-dimensional Stefan problems

In this paper, we present a critical comparison of the suitability of several numerical methods, level set, moving grid and phase field model, to address two well-known Stefan problems in phase transformation studies: melting of a pure phase and diffusional solid-state phase transformations in a binary system. Similarity solutions are applied to verify the numerical results. The comparison shows that the type of phase transformation considered determines the convenience of the numerical techniques. Finally, it is shown both numerically and analytically that the solid-solid phase transformation is a limiting case of the solid–liquid transformation.     

CFD methodology for sedimentation tanks: The effect of secondary phase on fluid phase using DPM coupled calculations

Computational Fluid Dynamics (CFD) methods are employed in order to simulate the 3D hydrodynamics and flow behaviour in a sedimentation tank. Unlike most of the previous numerical investigations, in the present paper the momentum exchange between the primary and the secondary phase is taken into account, using a Lagrangian method (discrete phase model) with two-way coupled calculations. By computing particle trajectories the proposed numerical model can track the momentum gained or lost by the particle stream that follows that trajectory and these quantities can be incorporated in the subsequent continuous phase calculations. Thus, while the continuous phase always impacts the discrete phase, the effect of the discrete phase trajectories on the continuum can be incorporated. This interchange affects fluid velocity, especially in the case of large particles sizes, which have a greater relaxation time in relation to the characteristic time of the tank. The present investigation compares a series of numerical simulations for a sedimentation tank with varying particle diameters and volume fractions, in order to identify the influence of the secondary phase to the primary phase and vice-versa and the way that this influence affects the efficiency of the tank.     

Differential phase space reconstructed for chaotic time series

A new numerical differential filter is built to estimate the numerical differential for a chaotic time series and then a differential phase space for the chaotic time series is reconstructed. Correlation dimensions, Lyapunov exponents and forecasting are discussed for the chaotic time series on the reconstructed differential phase space and on the delay phase space, respectively. Comparison results show that the numerical results on the differential phase space are better than that on the delay phase space.     

Glass transition and crystallization process of hard magnetic bulk Nd60Al10Fe20Co10 metallic glass

Glass transition and crystallization process of bulk Nd60Al10Fe20Co10 metallic glass were investigated by means of dynamic mechanical thermal analysis (DMTA), differential scanning calorimetry (DSC), X-ray diffraction (XRD) and scanning electronic microscopy (SEM). It is shown that the glass transition and onset crystallization temperature determined by DMTA at a heating rate of 0.167 K/s are 480 and 588 K respectively. The crystallization process of the metallic glass is concluded as follows: amorphous α→α′＋metastable FeNdAl novel phase →α′＋primary δ phase→primary δ phase＋eutectic δ phase Nd3Al phase＋Nd3Co phase. The appearance of hard magnetism in this alloy is ascribed to the presence of amorphous phase with highly relaxed structure. The hard magnetism disappeared after the eutectic crystallization of amorphous phase.     

Phase space Feynman path integrals via piecewise bicharacteristic paths and their semiclassical approximations

We give a fairly general class of functionals for which the phase space Feynman path integrals have a mathematically rigorous meaning. More precisely, for any functional belonging to our class, the time slicing approximation of the phase space path integral converges uniformly on compact subsets of the phase space. Our class of functionals is rich because it is closed under addition and multiplication. The interchange of the order with the Riemann integrals, the interchange of the order with a limit and the perturbation expansion formula hold in the phase space path integrals. The use of piecewise bicharacteristic paths naturally leads us to the semiclassical approximation on the phase space.     

一类非最小相位随机系统的离散自适应滑动模态控制

This paper presents the discrete adaptive sliding mode control of input-output non-minimum phase system in the presence of the stochastic disturbance. The non-minimum phase system can be transformed into a minimum phase system by a operator. According to the minimum phase system, the controller and the adaptive algorithm we designed ensures the stability of system and holds that the mean-square deviation from the sliding surface is minimized.     

估计的和理论的保证以及启动Ⅲ期试验的概率

预期的Ⅱ期试验通常会导致Ⅲ期试验失败.对于随机分为两个治疗组的随机对照的Ⅱ期和Ⅲ期试验,在假设正态分布响应的方差已知的情况下,我们解析性地获得了这三种情况下的估计的和理论的保证(无、加性和乘性偏差调整).在一些较小的假设下,我们证明了对这三种情况下的估计保证分别是Ⅱ期试验的每组患者数和Ⅱ期试验观察到的治疗效果的增加函数...     

Numerical simulation two phase flows of casting filling process using SOLA particle level set method

Mould filling process is a typical gas–liquid metal two phase flow phenomenon. Numerical simulation of the two phase flows of mould filling process can be used to properly predicate the back pressure effect, the gas entrapment defects, and better understand the complex motions of the gas phase and the liquid phase. In this paper, a novel sharp interface incompressible two phase numerical model for mould filling process is presented. A simple ghost fluid method like discretization method and a density evaluation method at face centers of finite difference staggered grid are proposed to overcome the difficulties when solving two phase Navier–Stokes equations with large-density ratio and large-viscosity ratio. A new mass conservation particle level set method is developed to capture the gas–liquid metal phase interface. The classical pressure-correction based SOLA algorithm is modified to solve the two phase Navier–Stokes equations. Two numerical tests including the Zalesak disk problem and the broken dam problem are used to demonstrate the accuracy of the present method. The numerical method is then adopted to simulate three mould filling examples including two high speed CCD camera imaging water filling experiments and an in situ X-ray imaging experiment of pure aluminum filling. The simulation results are in good agreement with the experiments.     

考虑晶粒几何形状的尺度效应

复合材料方法中的混合律方法在研究纳米晶体材料力学性能时得到了广泛的使用,准确得出各相的体积分数对该方法结果的准确性具有十分重要的影响．将纳米晶体看成由晶界、晶粒和三叉晶三相组成的复合材料,根据晶粒具有多面体的几何特点,用二维的三相复合的正多边形模型来研究纳米晶体力学性能的尺度效应,对于不同几何形状的晶粒采用对应的正多边形模型,这样我们就可以更加准确地得到各相的体积分数,从而更好地预测纳米晶体材料的力学性能．     

A Hybrid Metaheuristic for the Quadratic Assignment Problem

The quadratic assignment problem (QAP) is known to be NP-hard. We propose a hybrid metaheuristic called ANGEL to solve QAP. ANGEL combines the ant colony optimization (ACO), the genetic algorithm (GA) and a local search method (LS). There are two major phases in ANGEL, namely ACO phase and GA phase. Instead of starting from a population that consists of randomly generated chromosomes, GA has an initial population constructed by ACO in order to provide a good start. Pheromone acts as a feedback mechanism from GA phase to ACO phase. When GA phase reaches the termination criterion, control is transferred back to ACO phase. Then ACO utilizes pheromone updated by GA phase to explore solution space and produces a promising population for the next run of GA phase. The local search method is applied to improve the solutions obtained by ACO and GA. We also propose a new concept called the eugenic strategy intended to guide the genetic algorithm to evolve toward a better direction. We report the results of a comprehensive testing of ANGEL in solving QAP. Over a hundred instances of QAP benchmarks were tested and the results show that ANGEL is able to obtain the optimal solution with a high success rate of 90%. This work was supported in part by the National Science Council, R.O.C., under Contract NSC 91-2213-E-005-017.     

Singularities of the thermodynamic potential under second-order phase transitions due to the interband electron-phonon interaction

We consider phase transitions in crystals with a strong interband electron-phonon interaction. We investigate the thermodynamic potential of the system using the method of temperature Green’s functions, which takes quantum and thermal fluctuations into account. We show that in the absence of striction, these phase transitions are realized as a sequence of second-order phase transitions in each of which the thermodynamic potential has a logarithmic singularity, as in the Onsager model. This suggests that this singularity is characteristic of all second-order phase transitions. We show that the energy preference of the transition to the ordered phase is ensured by the electron coupling to coherent displacements of ions along normal coordinates of the phonon modes. We calculate the limit value of the energy decrease in the ordered phase compared with the symmetric phase as T → 0 K. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 157, No. 2, pp. 273–285, November, 2008.     

An optimal phase demodulator

This paper presents the results of investigation on the optimal performance of the phase demodulator using Bucy's non-linear filtering theory. A cyclically optimal non-linear phase estimation technique has been developed for modular phase tracking systems.

This paper considers a first-order Markov or Brownian motion phase process, because of faster computation in digital simulations and for the extensive research results available in Weiper optimum phase-lock loop for comparison purposes. However, the technique developed in this paper is equally applicable to higher-order phase process models.

The actual performance of the non-linear filter was investigated using Monte Carlo techniqes. The classical phase-lock loop was also simulated as a reference against which the optimal non-linear filter was compared. The cyclic phase non-linear filter appears to be about 0.7 dB better than the phase lock with respect to noise alone where the maximum conceivable improvement is 2.2 dB, or the excess noise relative to the ideal is about one-third less in dB than that of the phase-lock loop.     

Attractor Bifurcation and Phase Transition for Liquid <Superscript>4</Superscript>He

Experiments show that the liquid helium-4 has either superfluid phase or solid phase when temperature decreases or the pressure of the system changes. In this paper, we discuss the equations which govern these phase transitions and derive the criterions for these phase changes. Meanwhile, we give related approximate solutions and draw the phase diagram. In addition, we also prove that the liquid helium-4 bifurcates from the trivial solution to an attractor as parameters cross certain critical value. The topological structure of the bifurcated attractor is also illustrated.     

Renormalization group study of the melting of a two-dimensional system of collapsing hard disks

We consider the melting of a two-dimensional system of collapsing hard disks (a system with a hard-disk potential to which a repulsive step is added) for different values of the repulsive-step width. We calculate the system phase diagram by the method of the density functional in crystallization theory using equations of the Berezinskii–Kosterlitz–Thouless–Halperin–Nelson–Young theory to determine the lines of stability with respect to the dissociation of dislocation pairs, which corresponds to the continuous transition from the solid to the hexatic phase. We show that the crystal phase can melt via a continuous transition at low densities (the transition to the hexatic phase) with a subsequent transition from the hexatic phase to the isotropic liquid and via a first-order transition. Using the solution of renormalization group equations with the presence of singular defects (dislocations) in the system taken into account, we consider the influence of the renormalization of the elastic moduli on the form of the phase diagram.     

饱和多孔介质中骨架的应变局部化萌生条件

应用饱和多孔介质控制方程和Liapunov稳定理论,导出了固相应力和有效应力描述的多孔介质骨架应变局部化的萌生条件．不同应力形式表达的多孔介质基体的控制方程,相应的应变局部化萌生条件的表达形式也不尽相同,其原因源于骨架本构中固液两相之间相互作用的不同描述．应用得出的Terzaghi有效应力描述的应变局部化萌生条件,可以理论解释多孔介质中固、液两相不同相对运动出现的破坏方式,如管涌、滑坡和泥石流．应用简单算例说明了应变局部化条件的具体实施方法．     

Effect of the phase on the dynamics of a perturbed bouncing ball system

The phase control method is a non-feedback control technique which has been used for different purposes in continuous periodically driven dynamical systems. One of the main goals of this paper is to apply this control technique to the bouncing ball system, which can be seen as a paradigmatic periodically driven discrete dynamical system, and has a rather simple physical interpretation. The main idea is to apply a periodic control signal including a phase difference with respect to the periodic forcing of the initial system and to analyze its effect on the dynamics of the bouncing ball system. The numerical simulations we have carried out clearly show the strong effect of the phase of the control signal in suppressing or generating chaotic behavior and in changing the period of a periodic orbit. We have also analyzed the effect of the phase in the phenomenon of the crisis-induced intermittency, showing how the phase enhances the size of the attractor near a crisis and can induce the intermittent behavior. Finally we have analyzed the scaling behavior of the crisis by varying the phase difference between the perturbation and the external forcing.     

A Multi-phase-field model including inelastic deformation for solid state transformations

A multi phase field model is presented in order to take the plastic deformation during a solid state transformation into account and to investigate its effect on the transformation kinetics and morphology in a multi phase material. The model is formulated consistently with the multi phase field model for diffusional and surface driven phase transformations [1]. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

线性逻辑相空间的分层结构

本文研究了相空间的代数结构,定义了多层相空间,给出了一个相空间的最小子空间的性质.并据此讨论了线性重言式的某些特点.