A family of lower- and higher-order transversal linearization techniques in non-linear stochastic engineering dynamics

Sample pathwise numerical integration of noise-driven engineering dynamical systems cannot generally be performed beyond a limited level of accuracy, especially when the noise processes are modelled using (filtered) white noises. Recently, a locally transversal linearization (LTL) strategy has been proposed by the author (Proc Roy Soc London A 2001; 457 :539–566) for direct integration of deterministic and stochastic non-linear dynamical systems. The present effort is focussed on a host of extensions along with detailed theoretical error analyses of the linearization approach, especially as applied to problems in non-linear stochastic engineering dynamics. Thus, to begin with, estimates of local and global error orders in the basic LTL scheme are obtained separately for the displacement and velocity vectors when the system is driven either by a set of additive noises or by an arbitrary combination of (independently evolving) additive and multiplicative noises. Following this, a new family of higher-order LTL schemes is proposed in order to improve upon the basic LTL method and the associated error orders are established. A stepwise implementation of the lower- and higher-order versions of the LTL method, along with certain computational aspects, is also outlined. The proposed schemes are numerically illustrated, to a limited extent, for a single degree-of-freedom (SDOF) and a two degree-of-freedom (TDOF) non-linear engineering systems under additive and/or multiplicative white noise excitations. Copyright © 2004 John Wiley & Sons, Ltd.     

ON AVERAGED DIFFUSION EQUATIONS

Methods of obtaining averaged diffusion equations are considered in case of non-uniform profile of the velocity in a channel. With the flow of Couette as an example, the comparison of exact and approximate solutions (obtained by means of perturbation method) has been carried out. The peculiarities of function of residence time distribution of liquid in the flows with non-uniform velocity field are noted. It is shown that the distribution function moments values including the zero and first moments values would depend on the degree of the velocity profile irregularity, on efficiency of radial mixing in a system as well as on the averaging method. The averaged diffusion equations which have been found by means of perturbation method are the most general of proposed ones at present in the appropriate literature. The Taylor's model and Goldstein's hyperbolic equations are included in said averaged equations as particular cases. The table of the numerical values of first three moments of RTD-function necessary for determining of the model parameters is given. The problems of application of the obtained averaged equation for calculating real chemical apparatuses, e.g. reactors, are discussed.     

Mathematical analysis of viscometric (polymer) flow fields in capillaries: Taylor dispersion revisited

A general formulation for capillary flow of two miscible fluids – one a dilute plug of polymer fluid inserted into a fully developed Poiseuille flow of the other, a Newtonian stream – is examined for its long time behavior. Phenomenologically, the system evolves from an initial state, that of a plug within the boundaries of sharp, well defined fronts inside a Newtonian stream, to a more homogenized state in the very long time scale. This problem was addressed by G.I. Taylor but with regard to a system of two Newtonian fluids, leading to the well-known results commonly described as `Taylor axial dispersion'. In this paper, a general and systematic perturbation analysis is presented from which Taylor's result is recovered as a special case of a more general solution which applies to fluids incorporating elastic properties. In particular, the influence of viscoelasticity and (polymer) diffusivity on the observed pressure profile in the capillary conduit is examined. This effect is clearly separated out for small Peclet number flows using asymptotic and numerical analysis. The results identify the influence of fluid viscosity, elasticity, and diffusivity on the observed pressure profile and form the basis for the improved characterization of polymeric elasticity using capillaries – a finding that is of significant scientific and commercial interest. These results were obtained by the authors as a class of observations resulting from the perturbation analysis of forced-flow capillary devices in viscoelastic fluid property investigation.     

两个同心施转圆柱间的Taylor问题

本文推导了在修正双极坐标系下的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和讨论了它的基本流的扰动解     

Development of the meshless finite volume particle method with exact and efficient calculation of interparticle area

The Finite Volume Particle Method (FVPM) is a meshless method based on a definition of interparticle area which is closely analogous to cell face area in the classical finite volume method. In previous work, the interparticle area has been computed by numerical integration, which is a source of error and is extremely expensive. We show that if the particle weight or kernel function is defined as a discontinuous top-hat function, the particle interaction vectors may be evaluated exactly and efficiently. The new formulation reduces overall computational time by a factor between 6.4 and 8.2. In numerical experiments on a viscous flow with an analytical solution, the method converges under all conditions. Significantly, in contrast with standard FVPM and SPH, error depends on particle size but not on particle overlap (as long as the computational domain is completely covered by particles). The new method is shown to be superior to standard FVPM for shock tube flow and inviscid steady transonic flow. In benchmarking on a viscous multiphase flow application, FVPM with exact interparticle area is shown to be competitive with a mesh-based volume-of-fluid solver in terms of computational time required to resolve the structure of an interface.     

Thermal and kinetic properties of poly(lactic acid) and transglutaminase‐crosslinked wheat gluten blends

Wheat gluten (10 g) was crosslinked (XL) using 10 units of transglutaminase. Different blends of XL gluten and poly(lactic acid) (PLA) were mixed in a Brabender mixer at 180°C for 10 min. Neat PLA and blends were analyzed using modulated DSC (MDSC). Neat PLA displayed a glass transition (T_g) and exothermic (Cry) followed by endothermic (Mel) transitions. The profile showed a T_g of 0.46 J/g/°C, Cry with 29.9 J/g, whereas Mel exhibited 28.7 J/g. XL wheat gluten displayed one T_g with 0.45 J/g/°C. Samples were subjected to repeated heating and cooling cycles to show the level of compatibility between the two polymers. The activation energy (E_a) and pre‐exponential factor (Z) were determined according to Borchardt and Daniels (B/D) kinetics approach. The blends showed increased E_a values with an increase in the amount of XL gluten. In the presence of 5 and 20% XL gluten, the E_a of PLA increased from 150 to 200 kJ/mol, respectively. A higher number of cycles caused an increase in E_a. The T_g temperature of different PLA/XL gluten blends can be predicted by Gordon–Taylor equation and its modified forms. © 2007 Wiley Periodicals, Inc. J Appl Polym Sci 2007     

Mass transfer from Taylor bubbles rising in single capillaries

Gas-liquid mass transfer from Taylor bubbles rising in 1, 2 and 3 mm diameter capillaries of circular and square cross-sections was investigated for air-water system. The liquid-phase volumetric mass transfer coefficient k_La was obtained from experimental oxygen absorption dynamics. The experimental k_La values are in good agreement with the model developed by van Baten and Krishna (2004. Chemical Engineering Science 59, 2535-2545), with the additional assumption that the dominant mass transfer contribution is to the film surrounding the bubble.     

基于相关系数的决策树优化算法

通过分析ID3算法的基本原理及其多值偏向问题,提出了一种基于相关系数的决策树优化算法。首先通过引进相关系数对ID3算法进行改进,从而克服其多值偏向问题,然后运用数学中泰勒公式和麦克劳林公式的性质,对信息增益公式进行近似简化。通过具体数据的实例验证,说明优化后的ID3算法能够解决多值偏向问题。标准数据集UCI上的实验结果表明,在构建决策树的过程中,既提高了平均分类准确率,又降低了构建决策树的复杂度,从而还缩短了决策树的生成时间,当数据集中的样本数较大时,优化后的ID3算法的效率得到了明显的提高。     

自动微分在隐式曲线绘制中的应用

自动微分是用于计算多变量函数的导数和偏导数的一种微分技术,在给定一个多变量光滑函数值的程序代码后,可以很容易地利用自动微分来实现有关导数和偏导数的精确计算。将自动微分技术与泰勒方法相结合应用到计算机图形学领域隐式函数曲线绘制的细分算法中,并与未使用自动微分技术前的隐式曲线绘制方法作比较和分析,展示了自动微分方法在绘制隐式曲线方面的优势。     

一种高速正弦信号发生器的设计

文章提出了基于高速DSP芯片TMS320C5402实现正弦信号发生器的设计原理与方法,介绍了所设计的正弦信号发生器的硬件结构电路图和软件程序结构图。结合DSP硬件特性,通过使用泰勒级数展开法得到设定参数的正弦波波形输出,从而达到设计目的。