TO PROMOTE STUDENTS' INNOVATIVE ABILITY BY USING CLASSROOM DEMONSTRATION EXPERIMENTS1)

如何激发学生的创新意识、培养创新能力是当前教学的中心任务之一。本文笔者根据在流体力学课程教学中将简单易行、灵活多样的演示实验带进课堂,引发学生兴趣,启发创新思维方面的教学实践,总结了几点经验和体会。     

Effects of the method of identification of the diffusion coefficient on accuracy of modeling bound water transfer in wood

An alternative approach to determining the bound water diffusion coefficient is proposed. It comprises a method for solving the inverse diffusion problem, an improved algorithm for the bound-constrained optimization as well as an alternative submodel for the diffusion coefficient’s dependency on the bound water content. Identification of the diffusion coefficient for Scots pine wood (Pinus sylvestris L.) using the proposed inverse approach is presented. The accuracy of predicting the diffusion process with the use of the coefficient values determined by traditional sorption methods as well as by the inverse modeling approach is quantified. The similarity approach is used and the local and global relative errors are calculated. The results show that the inverse method provides valuable data on the bound water diffusion coefficient as well as on the boundary condition. The results of the identification can significantly improve the accuracy of mass transfer modeling as studied for drying processes in wood.     

各向同性弹性损伤的双标量描述

损伤状态的描述是损伤力学中仍未完善解决的基本问题．我们旨在对此问题就最简单的一种情形——各向同性弹性损伤，进行较为全面的研究．首先，指出了古典各向同性损伤理论中，基于应变等效假设，用单个标量损伤变量描述损伤状态的局限性．然后，建立了一个用两个标量损伤变量描述的各向同性弹性损伤模型．此模型解除了古典理论的局限，能完全描述各向同性弹性损伤，并且得到本文数值实验的验证．最后，将本文模型与现有细观力学结果连接，给出了宏细观损伤变量之间的关系，使得细观量可以通过宏观量来反映，建立了一个用细观损伤材料常数描述细观缺陷特征的损伤本构模型     

基于光声气体滤波的早期火灾烟气测量方法

火灾烟气中气体和烟雾浓度测量对火灾研究有着重要意义。本文提出一种利用光声技术同时在线测量CO浓度和烟雾消光的非接触式方法。根据气体滤波和粒子散射原理,采用两个光声腔分别作为CO测量腔和烟雾参比腔,在其与光源之间的开放路径中进行火灾产物的测量,通过测量两个光声腔信号,依据信号变化量计算CO浓度和烟雾消光。基于标准火的实验证明该方法可以实现烟雾颗粒和CO气体的复合测量。     

Rheological properties of branched polystyrenes: linear viscoelastic behavior

Oscillatory shear measurements on a series of branched graft polystyrenes (PS) synthesized by the macromonomer technique are presented. The graft PS have similar molar masses (M _w between 1.3×10⁵ g/mol and 2.4×10⁵ g/mol) and a polydispersity M _w /M _n around 2. The molar masses of the grafted side chains M _w,br range from 6.8×10³ g/mol to 5.8×10⁴ g/mol, which are well below and above the critical entanglement molar mass M _c of linear polystyrene. The average number of side chains per molecule ranges from 0.6 to 6.7. The oscillatory measurements follow the time–temperature superposition principle. The shift factors do not depend on the number of branches. The zero-shear viscosities of all graft PS are lower than those of linear PS with the same molar mass, which can be attributed to the smaller coil size of the branched molecules. It is shown that the influence of branching on the frequency dependence of the dynamic moduli is weak for all graft PS that were investigated, which can be explained by the low entanglement density.Electronic Supplementary Material Supplementary material is available for this article at This article has already been published online first (DOI: ). Due to an oversight at the publisher's, this version contained several mistakes. The article is herewith republished in its entirety as a "publisher's erratum".     

米箭沟尾矿坝料动力特性试验研究

本文对米箭沟尾矿坝料的组成、结构特征和物理力学性质进行分析,运用振动三轴试验确定尾矿坝料的动力特性参数,分析研究动力特性指标模型,得出动应力、动模量、阻尼比、动剪切摸量随动应变的变化规律,以及动应力与破坏振次,孔压比与振次比的变化规律,并提出符合试验结果的孔压模型,对该尾矿坝的地震动力分析及抗震稳定性评价提供一定科学依据和技术指导。     

Numerical and experimental investigation of wave dynamic processes in high-speed train/tunnels

Numerical and experimental investigation on wave dynamic processes induced by high-speed trains entering railway tunnels are presented. Experiments were conducted by using a 1:250 scaled train-tunnel simulator. Numerical simulations were carried out by solving the axisymmetric Euler equations with the dispersion-controlled scheme implemented with moving boundary conditions. Pressure histories at various positions inside the train-tunnel simulator at different distance measured from the entrance of the simulator are recorded both numerically and experimentally, and then compared with each other for two train speeds. After the validation of nonlinear wave phenomena, detailed numerical simulations were then conducted to account for the generation of compression waves near the entrance, the propagation of these waves along the train tunnel, and their gradual development into a weak shock wave. Four wave dynamic processes observed are interpreted by combining numerical results with experiments. They are: high-speed trains moving over a free terrain before entering railway tunnels; the abrupt-entering of high-speed trains into railway tunnels; the abrupt-entering of the tail of high-speed trains into railway tunnels; and the interaction of compression and expansion waves ahead of high-speed trains. The effects of train-tunnel configuration, such as the train length and the train-tunnel blockage ratio, on these wave processes have been investigated as well.     

Coherent structures and riblets

This work deals with the effect of the riblets on the coherent structures near the wall. The emphasis is put on the genesis of the quasi-streamwise vortices in the presence of the riblets. The quasi-streamwise vortices regenerate by the tilting of wall normal vorticity induced by prevailing structures. This requires a mechanism which leads to a temporal streamwise dependence near the elongated flow structures and to a subsequent formation of new wall normal vorticity. It is suggested here that the action of existing quasi-streamwise vortices on the sidewalls of wall normal vorticity may create a local, streamwise dependent spanwise velocity and therefore, a secondary wall normal vorticity field. A preliminary analysis of the set-up and the time and space development of this secondary three-dimensional flow associated with the regeneration mechanism, is given. An attempt is made, in order to explain the drag reduction performed by the riblets through an intermittent model, based on the protrusion height. Logical estimates of the amount of drag reduction are obtained. The differences between the mechanism suggested here and those based on forced control experiments are also discussed.     

Two-phase flow in heterogeneous porous media II: Numerical experiments for flow perpendicular to a stratified system

In this work, we make use of numerical experiments to explore our original theoretical analysis of two-phase flow in heterogeneous porous media (Quintard and Whitaker, 1988). The calculations were carried out with a two-region model of a stratified system, and the parameters were chosen be consistent with practical problems associated with groundwater flows and petroleum reservoir recovery processes. The comparison between theory (the large-scaled averaged equations) and experiment (numerical solution of the local volume averaged equations) has allowed us to identify conditions for which the quasi-static theory is acceptable and conditions for which a dynamic theory must be used. Byquasi-static we mean the following: (1) The local capillary pressure,everywhere in the averaging volume, can be set equal to the large-scale capillary pressure evaluated at the centroid of the averaging volume and (2) the large-scale capillary pressure is given by the difference between the large-scale pressures in the two immiscible phases, and is therefore independent of gravitational effects, flow effects and transient effects. Bydynamic, we simply mean a significant departure from the quasi-static condition, thus dynamic effects can be associated with gravitational effects, flow effects and transient effects. To be more precise about the quasi-static condition we need to refer to the relation between the local capillary pressure and the large-scale capillary pressure derived in Part I (Quintard and Whitaker, 1990). Herep _c ¦_y represents the local capillary pressure evaluated at a positiony relative to the centroid of the large-scale averaging volume, and {p _c }¦_x represents the large-scale capillary pressure evaluated at the centroid.In addition to{p _c } ^c being evaluated at the centroid, all averaged terms on the right-hand side of Equation (1) are evaluated at the centroid. We can now write the equations describing the quasi-static condition as , , This means that the fluids within an averaging volume are distributed according to the capillary pressure-saturation relationwith the capillary pressure held constant. It also means that the large-scale capillary pressure is devoid of any dynamic effects. Both of these conditions represent approximations (see Section 6 in Part I) and one of our main objectives in this paper is to learn something about the efficacy of these approximations. As a secondary objective we want to explore the influence of dynamic effects in terms of our original theory. In that development only the first four terms on the right hand side of Equation (1) appeared in the representation for the local capillary pressure. However, those terms will provide an indication of the influence of dynamic effects on the large-scale capillary pressure and the large-scale permeability tensor, and that information provides valuable guidance for future studies based on the theory presented in Part I.Roman Letters A scalar that maps {}*/t onto - A scalar that maps {}*/t onto - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - A interfacial area between the -region and the -region contained within, m² - a vector that maps ({}*/t) onto , m - a vector that maps ({}*/t) onto , m - b vector that maps ({p}– g) onto , m - b vector that maps ({p}– g) onto , m - B second order tensor that maps ({p}– g) onto , m² - B second order tensor that maps ({p}– g) onto , m² - c vector that maps ({}*/t) onto , m - c vector that maps ({}*/t) onto , m - C second order tensor that maps ({}*/t) onto , m² - C second order tensor that maps ({}*/t) onto . m² - D third order tensor that maps ( ) onto , m - D third order tensor that maps ( ) onto , m - D second order tensor that maps ( ) onto , m² - D second order tensor that maps ( ) onto , m² - E third order tensor that maps () onto , m - E third order tensor that maps () onto , m - E second order tensor that maps () onto - E second order tensor that maps () onto - p _c =(), capillary pressure relationship in the-region - p _c =(), capillary pressure relationship in the-region - g gravitational vector, m/s² - largest of either or - - - i unit base vector in thex-direction - I unit tensor - K local volume-averaged-phase permeability, m² - K local volume-averaged-phase permeability in the-region, m² - K local volume-averaged-phase permeability in the-region, m² - {K } large-scale intrinsic phase average permeability for the-phase, m² - K –{K }, large-scale spatial deviation for the-phase permeability, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K –{K }, large-scale spatial deviation for the-phase permeability in the-region, m² - K ^* large-scale permeability for the-phase, m² - L characteristic length associated with local volume-averaged quantities, m - characteristic length associated with large-scale averaged quantities, m - I _i i = 1, 2, 3, lattice vectors for a unit cell, m - l characteristic length associated with the-region, m - ; characteristic length associated with the-region, m - l _H characteristic length associated with a local heterogeneity, m - - n unit normal vector pointing from the-region toward the-region (n =–n ) - n unit normal vector pointing from the-region toward the-region (n =–n ) - p pressure in the-phase, N/m² - p local volume-averaged intrinsic phase average pressure in the-phase, N/m² - {p } large-scale intrinsic phase average pressure in the capillary region of the-phase, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p local volume-averaged intrinsic phase average pressure for the-phase in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - p –{p }, large scale spatial deviation for the-phase pressure in the-region, N/m² - P _c p –{p }, capillary pressure, N/m² - {p_c}^c large-scale capillary pressure, N/m² - r ₀ radius of the local averaging volume, m - R ₀ radius of the large-scale averaging volume, m - r position vector, m - , m - S /, local volume-averaged saturation for the-phase - S ^* {}*{}*, large-scale average saturation for the-phaset time, s - t time, s - u , m - U , m² - v -phase velocity vector, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - v local volume-averaged phase average velocity for the-phase in the-region, m/s - {v } large-scale intrinsic phase average velocity for the-phase in the capillary region of the-phase, m/s - {v } large-scale phase average velocity for the-phase in the capillary region of the-phase, m/s - v –{v }, large-scale spatial deviation for the-phase velocity, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - v –{v }, large-scale spatial deviation for the-phase velocity in the-region, m/s - V local averaging volume, m³ - V volume of the-phase in, m³ - V large-scale averaging volume, m³ - V capillary region for the-phase within, m³ - V capillary region for the-phase within, m³ - V _c intersection of m³ - V volume of the-region within, m³ - V volume of the-region within, m³ - V () capillary region for the-phase within the-region, m³ - V () capillary region for the-phase within the-region, m³ - V () , region in which the-phase is trapped at the irreducible saturation, m³ - y position vector relative to the centroid of the large-scale averaging volume, m Greek Letters local volume-averaged porosity - local volume-averaged volume fraction for the-phase - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region - local volume-averaged volume fraction for the-phase in the-region (This is directly related to the irreducible saturation.) - {} large-scale intrinsic phase average volume fraction for the-phase - {} large-scale phase average volume fraction for the-phase - {}* large-scale spatial average volume fraction for the-phase - –{}, large-scale spatial deviation for the-phase volume fraction - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - –{}, large-scale spatial deviation for the-phase volume fraction in the-region - a generic local volume-averaged quantity associated with the-phase - mass density of the-phase, kg/m³ - mass density of the-phase, kg/m³ - viscosity of the-phase, N s/m² - viscosity of the-phase, N s/m² - interfacial tension of the - phase system, N/m - , N/m - , volume fraction of the-phase capillary (active) region - , volume fraction of the-phase capillary (active) region - , volume fraction of the-region ( + =1) - , volume fraction of the-region ( + =1) - {p }– g, N/m³ - {p }– g, N/m³     

Experimentelle Bestimmung Binghamscher Stoffparameter

An experiment is described to determine the two Bingham material constants (yield stress _f and differential viscosity ) of viscoplastic fluids. The principle of the experiment is based on the falling-ball technique, where the stationary velocities of balls with different diameters and densities are measured. The required theory to calculate the Bingham material constants is illustrated. Experimental results of aqueous Carbopol 941-solutions are reported. These are listed in dependence of concentration in tables and diagrams.

---

Zusammenfassung Es wird ein Versuchsaufbau beschrieben, der es ermöglicht, die beiden Binghamschen Stoffparameter (Fließspannung _f und differentielle Viskosität ) einer viskoplastischen Flüssigkeit zu bestimmen. Der Versuch basiert auf dem Kugelfallprinzip, bei dem in einem Zylinder die stationäre Sinkgeschwindigkeit von Kugeln im Schwerefeld gemessen werden. Neben der Geschwindigkeit gehen das spezifische Gewicht der Flüssigkeit sowie die Geometrie und das spezifische Gewicht verschieden großer Kugeln in die Berechnung der Stoffparameter ein. Die zugehörige Theorie wird kurz erläutert. Im experimentellen Teil werden wäßrige Carbopol 941-Lösungen untersucht. Die Ergebnisse sind in Abhängigkeit der Konzentration tabellarisch angegeben und graphisch dargestellt.