Nucleate boiling of FC-87/FC-72 zeotropic mixtures on a horizontal copper disc

This paper presents nucleate boiling experimental results, at atmospheric pressure, for heat fluxes q ≤ 40 kW/m², for FC-87/FC-72 binary mixtures in molar fractions of 0/100, 25/75, 50/50, 75/25, 85/15 and 100/0, at saturation temperatures for pure fluids and bubble points for mixtures. The test section was an upward facing copper disc of 12 mm diameter and 1 mm thickness. The experimental heat transfer coefficient was compared with the correlations of Rohsenow (1952), as reported by Rohsenow et al. (Handbook of heat transfer, McGraw-Hill, New York, 1998), Stephan and Abdelsalam (Int J Heat Mass Transfer 23;73–78, 1978) and Cooper (Int Chem Eng Symp Ser 86:785–792, 1984) for pure fluids and the semi-empirical models of Stephan and Körner (Chem Ing Tech Jahrg 7:409–484, 1969), Thome (J Heat Transfer 104:474–478, 1982), Fujita et al. (1996), as reported by Rohsenow et al. (Handbook of heat transfer, McGraw-Hill, New York, 1998), Fujita and Tsutsui (Int J Heat Mass Transfer 37(1):291–302, 1994) and Calus and Leonidopoulos (Int J Heat Mass Transfer 17:249–256, 1973) for mixtures.     

Necessary and sufficient conditions for the absolute stability of discrete type lurie control system

ＮＥＣＥＳＳＡＲＹＡＮＤＳＵＦＦＩＣＩＥＮＴＣＯＮＤＩＴＩＯＮＳＦＯＲＴＨＥＡＢＳＯＬＵＴＥＳＴＡＢＩＬＩＴＹＯＦＤＩＳＣＲＥＴＥＴＹＰＥＬＵＲＩＥＣＯＮＴＲＯＬＳＹＳＴＥＭＺｈａｎｇＪｉｙｅ（张继业）（ＲｅｃｅｉｘｅｄＯｃｔ．５．１９９４Ｃｏｍｍｕ...     

Singular perturbation of nonlinear vector boundary value problem

In this paper we study the perturbed boundary value problem of the form dx/dt=f(x,y,t;ε), εdy/dt=g(x,y,t;ε), a_1(ε)x(0,ε)+a_2(ε)y(0,ε)=a(ε) b_1(ε)x(1,ε)+εb_2(ε)y(1,ε)=β(ε)in whichx,f,β∈E~m, y,g,a∈E~n, 0<ε(?)1and a_1(ε), a_2(ε), b_2(ε)and b_2(ε) are matrices of the appropriate size. Under the condition that g_y(t) is nonsingular and other suitable restrictions, the existence of the solution is proved, the asymptotic expansion of solution of order n is constructed, and the remainder term is estimated.     

高温气体热化学反应的DSMC微观模型分析

热化学耦合的非平衡现象一直是高温气体热化学问题研究的难点, 制约了诸如爆轰波胞格结构、低温点火速率等现象的分析. 本文以高温氮气离解和氢氧燃烧中的链式置换反应为例, 从微观反应概率、振动态指定的反应速率、热力学非平衡态的宏观反应速率、碰撞后的能量再分配等角度, 分析了直接蒙特卡罗模拟中的典型化学反应模型(TCE, VFD, QK模型)的微观动力学性质. 研究发现, 无论是高活化能的高温离解反应还是低活化能的链式置换反应, 实际参与反应的分子的振动能概率分布都偏离了平衡态的Boltzmann分布, 包含较强振动能额外影响的VFD模型可以很好地模拟高温离解反应, 而TCE (VFD的一个特例)和QK模型对活化能较低的链式置换反应的预测效果相对更好. 此外, 化学反应碰撞后的能量再分配应遵循微观细致平衡原理, 细微的偏差都可能造成平动能和振动能难以达到最终的平衡状态. 直接蒙特卡罗模拟的应用评估结果表明, 化学反应的振动倾向对热化学耦合过程产生了明显的影响, 特别是由于高振动能分子更多地参与了化学反应, 气体平均振动能的下降将影响后续化学反应的进行.      

三向受压混凝土的动态特性研究

为探究混凝土在三向受压状态下的动态特性,利用自行研制的大型多功能三轴材料试验机,进行不同应变速率(10^-5/s、10^-4/s、10^-3/s、10^-2/s)下混凝土不同定侧压比(1:1、2:1、3:1、4:1)的动态真三轴抗压试验,研究了混凝土在动态抗压下的强度和变形特性。结果表明:混凝土在三向受压状态下表现出明显的应变率效应,峰值应力随着应变速率的增加先减小后增大,峰值应力减小最大幅度为5.42%而后开始增大,最大增幅为18.22%。而峰值应力随着侧压比的增大而增大,到达3:1时应力速率敏感性降低。混凝土的峰值应变在低侧压比(1:1和2:1)时随着应变速率的增加而减小,而在侧压比较高(3:1和4:1)时先减小再增大。随着应变速率的增加,裂纹开始从骨料处产生,塑性应变增大,导致峰值应变增大;在动态加载条件下,峰值应变随着侧压比的增大有先增大后减小的趋势。较高侧压比时混凝土的弹性模量随着应变速率的增加具有增大的趋势,最大增幅为240.66%。应力-应变曲线上升段随着应变速率增大越来越陡峭。试件的破坏形态随着侧压比的增大从柱状破坏变为斜剪破坏。     

PMMA膨胀环动态拉伸碎裂实验研究

在强动载作用下, 脆性材料的碎裂问题是一个重要的研究课题, 而脆性材料在冲击拉伸载荷下的力学行为的实验研究相对较匮乏. 提出了一种动态拉伸断(碎)裂的液压膨胀环实验技术, 可用于准脆性/脆性材料的动态拉伸. 利用该技术对有机玻璃(PMMA)圆环试件进行了不同膨胀速度下的动态碎裂实验研究. 从回收碎片的断口形貌和碎片内部残余裂纹观察可知试件的破碎由环向拉伸应力造成, 碎片断口处发出的稀疏波会将周围的拉伸应力卸载, 从而抑制其他裂纹的进一步发展. 利用超高速相机记录了试件的膨胀碎裂过程, 利用DISAR激光速度干涉仪获得了试件外表面粒子的径向膨胀速度历史, 通过试件上的应变片获得了试件的应变历史和断裂应变. 实验结果表明: 在拉伸应变率150\mathrm{~}500s^{-1}范围, 材料的动态断裂应变低于准静态加载下的断裂应变, 体现出“动脆”现象; 随着加载应变率的提高, PMMA 材料的碎片尺寸减小; 无量纲化的PMMA圆环的平均碎片尺寸介于韧性碎裂模型和脆性碎裂模型的预测数值之间, 反映出材料的准脆性特性.      

Investigation of tensile deformation behavior of PC,ABS, and PC/ABS blends from low to high strain rates

This paper experimentally studies the tensile deformation behavior of polycarbonate (PC), acrylonitrile-butadiene-styrene (ABS), and PC/ABS blends (with the blending ratios of PC to ABS being 80:20, 60:40, 50:50, and 40:60) from low to high strain rates. Using the universal MTS-810 machine and the split Hopkinson tension bar (SHTB) testing system, the quasi-static and impact tension tests are carried out at the room temperature. The curves of the true stress and the true strain are obtained. The deformation behaviors of PC, ABS, and PC/ABS blends are characterized in detail. The linear relationship between the strain rate and the yielding stress is given.     

Kinetics of ultrafine coating of SiO2 nanoparticles, suspended in non-thermal CH4/C2H6 plasma

Capacitively coupled RF discharges in methane and ethane (1000–3000 Pa, 5–15 W/cm³) were used for the in-flight coating of SiO₂ nanopowders (NP) with an a-C:H layer [A. Kouprine, F. Gitzhofer, M. Boulos, A. Fridman, Polymer-like C:H thin film coating of nanopowders in capacitively coupled RF discharge, Plasma Chemistry and Plasma Processing 24 (2) (2004) 189–215]. In this present work a model of the chemical kinetics of this plasma processing has been developed, based on the GRI-Mech 3.0 mechanism [G.P. Smith, D.M. Golden, M. Frenklach, Gas Research Institute, Detailed chemical reaction mechanism. Available from: ], and calculations are performed, using the Chemkin-II code. The model includes the formation of two solid phases: an amorphous C:H coating, deposited on the suspended NP, and soot, incepted in the gas phase. Non-equilibrium plasma effects are taken into account by the rates of radical species’ production, using the reference data on the “G-value”. The CH₃ and H species appear to have key roles in surface activation and the nanofilm growth. The results demonstrate the competitive character of the formation of the two solid phases: the C:H coating is dominant under T_gas < 1400–1600 K, depending on the conditions, and soot commencing to dominate at higher temperatures. This division corresponds to the passage from non-thermal plasma conditions, with the dominance of single-valent radicals and solid phase growth on the substrate only, to the plasma thermalisation with the formation of two-valent radicals and v. inception of the solid phase. The results are validated experimentally by the observation of soot formation as a function of the gas temperature, by the mass spectrograph data on the gas composition after plasma reforming, and by the reference data on coating growth rates and induction times for soot formation.     

Coherent fine scale eddies in turbulence transition of spatially-developing mixing layer

To investigate the relationship between characteristics of the coherent fine scale eddy and a laminar–turbulent transition, a direct numerical simulation (DNS) of a spatially-developing turbulent mixing layer with Re_ω,0 = 700 was conducted. On the onset of the transition, strong coherent fine scale eddies appears in the mixing layer. The most expected value of maximum azimuthal velocity of the eddy is 2.0 times Kolmogorov velocity (u_k), and decreases to 1.2u_k, which is an asymptotic value in the fully-developed state, through the transition. The energy dissipation rate around the eddy is twice as high compared with that in the fully-developed state. However, the most expected diameter and eigenvalues ratio of strain rate acting on the coherent fine scale eddy are maintained to be 8 times Kolmogorov length (η) and :β:γ = −5:1:4 in the transition process. In addition to Kelvin–Helmholtz rollers, rib structures do not disappear in the transition process and are composed of lots of coherent fine scale eddies in the fully-developed state instead of a single eddy observed in early stage of the transition or in laminar flow.     

From crack deflection to lattice vibrations—macro to atomistic examination of dynamic cleavage fracture

A set of cleavage experiments with strip-shaped single-crystal silicon specimens subjected to three-point bending is reported. The experiments enabled examination of the relationships between the dynamic energy release rate, the velocity, the orientation-dependent cleavage energy, and the cleavage plane of propagation.Dynamic crack propagation experiments show that when a [0 0 1] silicon single crystal is fractured under three-point bending at ‘parallel’ velocity (directly measured at the bottom surface of the specimen) of up to , it prefers to cleave along the vertical (1 1 0) plane, while when the specimen is fractured under the same conditions but at a velocity higher than , it cleaves along the inclined (1 1 1) plane. At intermediate velocities, the crack will deflect from the (1 1 0) plane to the (1 1 1) plane. Crack velocity was determined by the initial notch length. The local (calculated) velocity of deflection between the cleavage planes ranges from , for a crack propagating on the (1 1 0) plane in the direction, to about , for a crack on the (1 1 0) plane, but in the [0 0 1] direction.It is suggested that the cause of the deflection phenomenon is the anisotropic, velocity-dependent cleavage energy, resulted phonon radiation caused by anisotropic, velocity-dependent lattice vibrations. We have studied the effect of material properties and propose selection criteria to explain the deflection phenomenon: the crack will deflect to the plane of least-energy, for which G−Γ_i(V)=max, or to the plane with maximum crack tip velocity, V_i(Γ)=max.     

An experimental study on developing air-water two-phase flow along a large vertical pipe: effect of air injection method

The flow structure in a developing air-water two-phase flow was investigated experimentally along a large vertical pipe (inner diameter, D_h: 0.48 m, ratio of length of flow path L to D_h: about 4.2). Two air injection methods (porous sinter injection and nozzle injection) were adopted to realize an extremely different flow structure in the developing region. The flow rate condition in the test section was as follows: superficial air velocity: 0.02–0.87 m/s (at atmospheric pressure) and superficial water velocity: 0.01–0.2 0.01–0.2 m/s, which covers the range of bubbly to slug flow in a small-scale pipe (D_h about 0.05 m).

No air slugs occupying the flow path were recognized in this experiment regardless of the air injection methods even under the condition where slug flow is realized in the small-scale pipe. In the lower half of the test section, the axial distribution of sectional differential pressure and the radial distribution of local void fraction showed peculiar distributions depending on the air injection methods. However, in the upper half of the test section, the effects of the air injection methods are small in respect of the shapes of the differential pressure distribution and the phase distribution. The comparison of sectional void fraction near the top of the test section with Kataoka's correlation indicated that the distribution parameter of the drift-flux model should be modeled including the effect of D_h and the bubble size distribution is affected by the air injection methods. The bubble size distribution is considered to be affected also by L/D_h based on comparison of results with Hills' correlation.     

Changes in some mechanical properties of a loamy soil under the influence of mechanized forest exploitation in a beech forest of central Belgium

Modification of some soil mechanical properties (penetration resistance and consolidation pressure) induced by vehicle compaction during mechanized forest exploitation was studied in an acid and loamy leached forest soil of the loessic belt of central Belgium. In situ penetration tests and laboratory Bishop–Wesley cell tests were undertaken for the two main soil horizons of a beech high-forest, i.e. the eluvial E horizon (5–30 cm depth) and the underlying clay-enriched B_t horizon (30–60 cm depth). Both undisturbed and wheel-rutted soil areas were studied (E and B_t horizons vs. E_g and B_tg horizons).

Results show that: The experimental overconsolidation pressure of the eluvial reference horizon (E) is about 50 kPa higher than the value calculated from soil overburden pressure; this probably results from suction action during dry periods. The clay-enriched reference horizon (B_t) shows the same trends. In wheel-rutted areas, seven years after logging operations, the E_g horizon memorizes only 14.5% of the wheel induced stress due to forest machinery.

In the compacted B_tg horizon, the experimental overconsolidation pressure represents 96% of the exerted theoretical stresses due to harvesting actions. The good recording of the exerted stresses, after seven years, can be explained by: (1) The B_tg depth which keeps it from seasonal variations i.e. from desiccation–moistening or freeze–thaw cycling; (2) amorphous and free iron accumulation inducing a “glue” effect of the B_tg soil matrix, which could stabilize the soil structure and prevent recovery to initial conditions. These results provide clear evidence that on loessic materials, soil compaction due to logging operations leads to modifications in both physical (bulk density, total porosity) and mechanical (penetration resistance and consolidation pressure) soil properties.     

An evolutionary optimization of diffuser shapes based on CFD simulations

An efficient and robust algorithm is presented for the optimum design of plane symmetric diffusers handling incompressible turbulent flow. The indigenously developed algorithm uses the CFD software: Fluent for the hydrodynamic analysis and employs a genetic algorithm (GA) for optimization. For a prescribed inlet velocity and outlet pressure, pressure recovery coefficient C (the objective function) is estimated computationally for various design options. The CFD software and the GA have been combined in a monolithic platform for a fully automated operation using some special control commands. Based on the developed algorithm, an extensive exercise has been made to optimize the diffuser shape. Different methodologies have been adopted to create a large number of design options. Interestingly, not much difference has been noted in the optimum C values obtained through different approaches. However, in all the approaches, a better design has been obtained through a proper selection of the number of design variables. Finally, the effect of diffuser length on the optimum shape has also been studied. Copyright © 2009 John Wiley & Sons, Ltd.     

Discrete dislocation plasticity analysis of the wedge indentation of films

The plane strain indentation of single crystal films on a rigid substrate by a rigid wedge indenter is analyzed using discrete dislocation plasticity. The crystals have three slip systems at ±35.3^° and 90^° with respect to the indentation direction. The analyses are carried out for three values of the film thickness, 2, 10 and , and with the dislocations all of edge character modeled as line singularities in a linear elastic material. The lattice resistance to dislocation motion, dislocation nucleation, dislocation interaction with obstacles and dislocation annihilation are incorporated through a set of constitutive rules. Over the range of indentation depths considered, the indentation pressure for the 10 and thick films decreases with increasing contact size and attains a contact size-independent value for contact lengths . On the other hand, for the films, the indentation pressure first decreases with increasing contact size and subsequently increases as the plastic zone reaches the rigid substrate. For the 10 and thick films sink-in occurs around the indenter, while pile-up occurs in the film when the plastic zone reaches the substrate. Comparisons are made with predictions obtained from other formulations: (i) the contact size-independent indentation pressure is compared with that given by continuum crystal plasticity; (ii) the scaling of the indentation pressure with indentation depth is compared with the relation proposed by Nix and Gao [1998. Indentation size effects in crystalline materials: a law for strain gradient plasticity. J. Mech. Phys. Solids 43, 411-423]; and (iii) the computed contact area is compared with that obtained from the estimation procedure of Oliver and Pharr [1992. An improved technique for determining hardness and elastic-modulus using load and displacement sensing indentation experiments, J. Mater. Res. 7, 1564-1583].     

Numerical Investigations of the Steady State Relative Permeability of a Simplified Porous Medium

The purpose of this paper is to investigate, by flow simulations in a uniform pore-space geometry, how the co and countercurrent steady state relative permeabilities depend on the following parameters: phase saturation, wettability, driving force and viscosity ratio. The main results are as follows: (i) with few exceptions, relative permeabilities are convex functions of saturation; (ii) the cocurrent relative permeabilities increase while the countercurrent ones decrease with the driving force; (iii) with one exception, phase 2 relative permeabilities decrease and phase 1 relative permeabilities increase with the viscosity ratio M=₁/₂; (iv) the countercurrent relative permeabilities are always less than the cocurrent ones, the difference being partly due to the opposing effect of the viscous coupling, and partly to different levels of capillary forces; (v) the pore-level saturation distribution, and hence the size of the viscous coupling, can be very different between the cocurrent and the countercurrent cases so that it is in general incorrect to estimate the full mobility tensor from cocurrent and countercurrent steady state experiments, as suggested by Bentsen and Manai (1993).(Now at AS Norske Shell, Norway.) e-mail:     

Comparison between the dynamic moduli of fully non-linear and quasilinear viscoelastic materials

This paper deals with the non-linear viscoelastic model with multiple hereditary integrals (MHI) in the frequency domain and the conditions that it reduces to single hereditary integral or the quasilinear viscoelastic (QLV) model. It is shown that when the higher order complex moduli are related to the first-order modulus as the MHI model reduces to the QLV model. The coefficients of quasilinearity should be real and independent of amplitude and frequency.     

SiO2-coated pure anatase TiO2 catalysts for enhanced photo-oxidation of naphthalene and anthracene

Our current efforts reveal the preparation of SiO₂@TiO₂ nanocomposites having different thicknesses of silica shell and the relationship to photocatalytic activity (PCA) for the photo-oxidation of naphthalene and anthracene. The presence of SiO₂ coating over TiO₂ surface was demonstrated by FT-IR analysis, with peaks corresponding to SiOSi (1081 cm⁻¹) and SiOTi (950 cm⁻¹) bonds observed. High-resolution transmission electron microscopy analysis confirmed the presence of SiO₂ in the as-prepared nanocomposites and the amount of Si, Ti, and O was determined by energy dispersive X-ray spectroscopy analysis. Increasing the SiO₂ shell thickness increases the surface area of the nanocomposites (69–235 m²/g), which enhances naphthalene/anthracene adsorption. However, the observed PCA trend presents an inverse correlation to the adsorption studies, where the as-prepared samples possessing the highest surface areas exhibited the least PCA, while catalysts having lower surface areas (among silica coated samples) displayed the highest PCA in the degradation of naphthalene and anthracene to CO₂. Despite complete degradation of naphthalene and anthracene, incomplete mineralization occurred, ascribed to the formation of various intermediates, identified by GC–MS analysis.