Electromagnetohydrodynamic flows and mass transport in curved rectangular microchannels

Curved microchannels are often encountered in lab-on-chip systems because the effective axial channel lengths of such channels are often larger than those of straight microchannels for a given per unit chip length. In this paper, the effective diffusivity of a neutral solute in an oscillating electromagnetohydrodynamic(EMHD)flow through a curved rectangular microchannel is investigated theoretically. The flow is assumed as a creeping flow due to the extremely low Reynolds number in such microflow systems. Through the theoretical analysis, we find that the effective diffusivity primarily depends on five dimensionless parameters, i.e., the curvature ratio of the curved channel, the Schmidt number, the tidal displacement, the angular Reynolds number, and the dimensionless electric field strength parameter. Based on the obtained results, we can precisely control the mass transfer characteristics of the EMHD flow in a curved rectangular microchannel by appropriately altering the corresponding parameter values.     

Homogeneous deformation of a two-fluid plasma with allowance for electron inertia

A method of finding the homogeneous deformations of a two-fluid plasma with allowance for the electron inertia is proposed. By homogeneous deformation is meant an axisymmetric plasma flow with a linear dependence of the radial velocity on the radius. Three families of homogeneous deformations are found using this method. One of these families, consisting of deformations with an arbitrary law of variation of the total current, is of particular interest with reference to plasma column dynamics. The method proposed is based on the reduction of the equations of two-fluid plasma dynamics to single-fluid equations of the hydrodynamic type (the equations of electromagnetic hydrodynamics (EMHD)) with a non-diagonal internal stress tensor, three-parameter thermodynamics, and a nonlocal form of the generalized Ohm’s law. Possible applications of the exact solutions found to the analysis of the data obtained using certain experimental apparatus are discussed.     

Electromagnetohydrodynamic (EMHD) flow of fractional viscoelastic fluids in a microchannel

This study investigates the electromagnetohydrodynamic (EMHD) flow of fractional viscoelastic fluids through a microchannel under the Navier slip boundary condition. The flow is driven by the pressure gradient and electromagnetic force where the electric field is applied horizontally, and the magnetic field is vertically (upward or downward). When the electric field direction is consistent with the pressure gradient direction, the changes of the steady flow rate and velocity with the Hartmann number Ha are irrelevant to the direction of the magnetic field (upward or downward). The steady flow rate decreases monotonically to zero with the increase in Ha. In contrast, when the direction of the electric field differs from the pressure gradient direction, the flow behavior depends on the direction of the magnetic field, i.e., symmetry breaking occurs. Specifically, when the magnetic field is vertically upward, the steady flow rate increases first and then decreases with Ha. When the magnetic field is reversed, the steady flow rate first reduces to zero as Ha increases from zero. As Ha continues to increase, the steady flow rate (velocity) increases in the opposite direction and then decreases, and finally drops to zero for larger Ha. The increase in the fractional calculus parameter α or Deborah number De makes it take longer for the flow rate (velocity) to reach the steady state. In addition, the increase in the strength of the magnetic field or electric field, or in the pressure gradient tends to accelerate the slip velocity at the walls. On the other hand, the increase in the thickness of the electric double-layer tends to reduce it.

     

基于突变理论的液化评估方法

作为众多突变类型的一种形式,尖点突变模型由于其优良的性质而广泛应用于生物学和社会科学,近年来又有越来越多的学者用尖点突变模型来研究工程实际问题。本文基于地震液化的基本机理和突变理论的特点,以土体的抗液化能力和土体得到的地震能量为控制变量,以孔隙水压力发展为状态变量,建立了地震诱导的孔隙水压力发展模型。宏观震害调查结果和本文的计算结果对比表明,本文建议的方法方便、可靠,是一种理想的液化评估方法。     

Real-time Measurement of Ultrasonic Waves at Bolted Joints under Fatigue Testing

The real-time measurement of the surface acoustic wave (SAW) was carried out during the fatigue testing of the bolted joints of aluminum alloy plates with a frequency of 0.001 Hz. SAW distributions in the bolted region were compared with those obtained by the synchronized measurement, in which the ultrasonic wave was generated in synchronization with a loading cycle in the fatigue testing with a frequency of 10 Hz. At different numbers of fatigue cycles, the intensity of the reflection from the fatigue crack obtained by the real-time measurement was in good agreement with that obtained by the synchronized measurement. In the real-time measurement, the reflection intensity and profile changed with the stress level in a loading cycle, which were in good agreement with those obtained by the synchronized SAW measurement. From these results, it was confirmed that the SAW distributions obtained by the synchronized measurement is coincident with ones obtained by the real-time measurement in one loading cycle, and is not influenced by the measurement conditions.     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

用拉格朗日函数分析不同模量静不定桁架内力

在外载荷作用下的不同模量静不定桁架平衡问题,是任意有限多个自变量的多元函数在任意有限多个约束条件下的极值问题,对采用拉格朗日乘数法求解此类极值问题进行了数学证明．通过求解不同模量静不定桁架极限载荷的几个算例,阐述拉格朗日乘数法在计算不同模量静不定桁架极限载荷中的应用．研究结果表明：采用拉格朗日乘数法求解不同模量静不定桁架极限载荷的通用性较强,用拉格朗日乘数法求解不同模量静不定桁架极限载荷的方法不但克服了常规方法需利用几何关系建立协调方程的缺陷,且具有力学概念清晰直观、计算过程简便、便于工程设计人员在实际中掌握和应用．     

Green function on two-phase saturated medium by concentrated force in two-dimensional displacement field

The Green function on two-phase saturated medium by concentrated force has a broad and important use in seismology, seismic engineering, soil mechanics, geophysics, dynamic foundation theory and so on. According to the Green function on two-phase saturated medium by concentrated force in three-dimentional displacement field obtained by Ding Bo-yang et al. , it gives out the Green function in two-dimensional displacement field by infinite integral method along x3 -direction derived by De Hoop and Manolis. The method adopted in the thesis is simpler. The result will be simplified to the boundary element method of dynamic problem.     

A crystal plasticity theory for latent hardening by glide twinning through dislocation transmutation and twin accommodation effects

Imagine a residual glide twin interface advancing in a grain under the action of a monotonic stress. Close to the grain boundary, the shape change caused by the twin is partly accommodated by kinks and partly by slip emissions in the parent; the process is known as accommodation effects. When reached by the twin interface, slip dislocations in the parent undergo twinning shear. The twinning shear extracts from the parent dislocation a twinning disconnection, and thereby releases a transmuted dislocation in the twin. Transmutation populates the twin with dislocations of diverse modes. If the twin deforms by double twinning, double-transmutation occurs even if the twin retwins by the same mode or detwins by a stress reversal. If the twin deforms only by slip, transmutation is single. Whether single or double, dislocation transmutation is irreversible. The multiplicity of dislocation modes increases upon strain, since the twin finds more dislocations to transmute upon further slip of the parent and further growth of the twin. Thus, the process induces an increasing latent hardening rate in the twin. Under profuse twinning conditions, typical of double-lattice structures, this rate-increasing latent hardening combined with crystal rotation to hard orientations by twinning is consistent with a regime of increasing hardening rate, known as Regime II or Regime B. In this paper, we formulate governing equation of the above transmutation and accommodation effects in a crystal plasticity framework. We use the dislocation density based model originally proposed by Beyerlein and Tomé (2008) to derive the effect of latent hardening in a transmuting twin. The theory is expected to contribute to surmounting the difficulty that current models have to simultaneously predict under profuse twinning, the stress-strain curves, intermediate deformation textures, and intermediate twin volume fractions.     

乳化炸药密度对其压力减敏的影响

在相同的乳胶基质内分别添加2%、3%、4%、5%的空心玻璃微球,2%、3%、4%、5%的膨胀珍珠岩和0.10%、0.15%、0.20%、0.25%的化学发泡剂制备了3组乳化炸药。测试了各组乳化炸药未受压时和受冲击波动态压缩作用之后的水下爆炸冲击波,以波峰值计算他们的压力减敏度。比较了各组乳化炸药的压力减敏度大小,分析了密度对乳化炸药压力减敏的影响。结果表明:密度较大的乳化炸药压力减敏度较小,膨胀珍珠岩或化学发泡剂的添加量越大,乳化炸药越容易产生压力减敏作用;空心玻璃微球的添加量由2%增加到4%,乳化炸药的压力减敏度增加,但当空心玻璃微球的添加量由4%增加到5%后,乳化炸药的压力减敏度反而减小。密度对乳化炸药压力减敏影响的主要原因在于密度调整剂周围的乳胶基质破乳。     

A Numerical Study of the Distribution of Water in Partially Saturated Porous Rock

The distribution of water and air phases in small blocks of porous sandstone is examined by using a simulated annealing technique that finds the minimum interfacial energy distributions at different saturations. Simulations are based on existing sandstone microstructures that were determined by X-ray microtomography. At low saturations, some of the water is distributed in films along the walls of larger pore spaces, and connects to pendular structures in the crevices and smaller pores. As the amount of water in the pores increases the water films become thicker and pores fill from the pendular structures. The distribution of water voxels in the pore space is examined by calculating interfacial areas, by classifying water voxels as to whether they lie within films or clusters, and by determining the size and distribution of these film clusters. An exponential relationship is found between the fraction of water voxels in the films and the degree of saturation. In addition, the dependency of small-sample electrical conductivity on saturation is examined by using a random walk method.     

Shock-vortex shear-layer interaction in the transonic flow around a supercritical airfoil at high Reynolds number in buffet conditions

This paper provides a conceptual analysis and a computational model for how the unsteady ‘buffeting’ phenomenon develops in transonic, low incidence flow around a supercritical aerofoil, the OAT15A, at Reynolds number of 3.3 million. It is shown how a low-frequency buffet mode is amplified in the shock-wave region and then develops upstream and downstream interaction with the alternating von Kármán eddies in the wake past the trailing-edge as well as with the shear-layer, Kelvin–Helmholtz vortices. These interactions are tracked by wavelet analysis, autoregressive (AR) modelling and by Proper Orthogonal Decomposition. The frequency modulation of the trailing-edge instability modes is shown in the spectra and in the wall-pressure fluctuations. The amplitude modulation of the buffet and von Kármán modes has been also quantified by POD analysis. The thinning of the shear layers, both at the outer edge of the turbulent boundary layers and the wake, caused by an ‘eddy-blocking’ mechanism is modelled by stochastic forcing of the turbulent kinetic energy and dissipation, by small-scale straining of the higher-order POD modes. The benefits from thinning the shear-layers by taking into account the interfacial dynamics are clearly shown in the velocity profiles, and wall pressure distribution in comparison with the experimental data.     

Relationships of field traffic and tillage to corn yields and soil properties

Esperiements were conducted during the summer of 1979 in which field plots oon s Ste. Rosalie clay soil and a Ste. Amable sandy loam soil were subjected to different levels of compaction by machinery, and subsequently treated by moldboard plowing and discing, chiselling and subsoiling by a winged tool. A silage corn crop was grown on all plots and measurements were made of soil bulk densities, penetration resistance of soils and plant yields. The results indicated that the compaction of the soil, if not subsequently loosened by a tillage operation, caused a marked reduction in plant yields. A nnarrow range of dry bulk density produced the optimum silage corn yields in the two experimental soils. The soil densities in this range were obtained by any of the three tillage treatments, as well as by the rototiller treatment, without machinery traffic.     

Elastic Equilibrium of the Half-Space Revisited. Mindlin and Boussinesq Problems

The displacement caused in an isotropic elastic half-space by a point force localized on or beneath its surface is calculated here by a new method. These classical problems are known as Boussinesq and, respectively, Mindlin problems. The motivation for the present work resides in the fact that the original solutions involve some particular procedures, required by the complexity of the boundary conditions, which may limit their general application. The solutions presented here are obtained by including in a generalized Poisson equation the values of the function and its derivatives on the boundary, and by using in-plane Fourier transforms. This method is general and can be extended to other, similar problems.     

A self- similar solution of a shock propagation in a mixture of a non-ideal gas and small solid particles

Similarity solutions are obtained for unsteady, one-dimensional self-similar flow behind a strong shock wave, driven by a moving piston, in a dusty gas. The dusty gas is assumed to consist of a mixture of small solid particles and a non-ideal gas, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-condition is maintained and variable energy input is continuously supplied by the piston. Solutions are obtained under both the isothermal and adiabatic conditions of the flow-field. The spherical case is worked out in detail to investigate to what extent the flow-field behind the shock is influenced by the non-idealness of the gas in the mixture as well as by the mass concentration of the solid particles, by the ratio of density of the solid particles to the initial density of the mixture and by the energy input due to moving piston. A comparison is also made between isothermal and adiabatic cases.     

Diffusion of Tracer in Altered Tonalite: Experiments and Simulations with Heterogeneous Distribution of Porosity

Numerical time-domain-diffusion simulations were used for studying the diffusion behavior of tracer molecules in rock matrix with homogeneous and heterogeneous porosity. As the heterogeneous sample in these simulations, a 3D tomographic image of altered tonalite was used, in which the mineral components and the pores resolved by X-ray microtomography were represented by their respective intragranular porosities determined previously by the ¹⁴C-PMMA method. The apparent diffusion coefficient of a tracer in altered tonalite was determined experimentally, and was then used in the simulations. In the altered tonalite analyzed, inclusion of heterogeneity in the porosity increased the diffusion coefficient by 16 %. Altered and pristine feldspar was the main mineral component in the sample (72 %), and it also provided the dominant contribution to tracer diffusion, explaining alone 52 % of the diffusion coefficient. The large pores resolved by microtomography (6 %) and altered and pristine mica (22 %) gave an equal contribution to the diffusion coefficient. The simulation method applied was also validated by comparing the results to both an analytical and a numerical solution to the diffusion equation in a homogenous medium. In addition, the method was compared to discrete-time random-walk simulations in the case of randomly placed overlapping spheres.     

Linear and non-linear rheological behaviour of cement and silica suspensions. Effect of polymer addition

Oil well cement pastes and model silica suspensions demonstrate similar rheology: in oscillatory shear, beyond a critical stress, a sharp transition is ob- served between gel and liquid behaviour. In creep tests, an apparent yield stress and shear-thinning are followed by the appearance of shear thickening. The minimum viscosity measured in steady shear is close in value to the complex viscosity obtained from oscillatory measurements. The observations can be explained by the formation of liquid trapping aggregates whose compactness may be estimated by fitting the Tsenoglou model, and whose cohesion is reflected in the rigidity of the gel and in the critical strain (or stress) of gel dissolution. Substituting cement or silica particles by polymer redispersible powder causes a decrease of the storage modulus in the gel state and a lower viscosity, while leaving the general features of the flow curve unchanged. Decrease in material rigidity may be due to a weaker inter-particle attraction generated by the polymer presence. The decrease in viscosity is explained by a lessening of water entrapped within the aggregates, which now contain polymer particles which are less hydrophilic than either cement or silica.     

Analysis of liquid circulation and mixing in a partitioned electrolytic tank

The performance of an electrochemical process depends critically on the mobility of the reacting species or ions towards the electrode surface. In this work, a partitioned electrolytic cell is studied. Here the fluid flow is induced by gases which evolve at the electrode surface. The liquid circulation induced by the rising bubbles is primarily responsible for mixing. In this study, the liquid circulation in a cell where an alkaline solution of water is electrolyzed using different Nickel designs of electrodes is investigated using PIV. For each electrode, the optimum operating conditions such as voltage and concentration of electrolyte which resulted in good mixing are found. The flow-field is quantified by calculating time averaged velocity profiles along the horizontal line and by analyzing the temporal variation of liquid velocity at a point. It is found that there are differences in the circulation and hence vorticity in the two compartments, anode and cathode. The effect of gas evolution on mixing between the two chambers is studied by taking uric acid in the cathode half and NaOH in the anode half. The flow induced by the evolved gas bubbles leads to convective mixing in the two chambers. The mixing time is calculated by measuring the variation of current with time under potentiostatic conditions. This is verified by measuring the pH in anode and cathode compartments during the electrolysis.     

Interaction of an Ellipsoidal Inclusion with an Elliptic Crack in an Elastic Medium

An approach is considered to how to allow for the interaction between an ellipsoidal heterogeneity (inclusion) and an elliptic crack in an elastic medium. Using the superposition of perturbed stress states, the boundary conditions are satisfied on the ellipsoidal surface by the method of equivalent inclusion and on the crack surface by the least-squares method. A numerical analysis is carried out. Typical mechanical effects are revealed. In the calculations, the stress state near the ellipsoidal heterogeneity is approximated by a polynomial of the second degree in Cartesian coordinates, whereas the load on the crack surface is simulated by a polynomial of the fourth degree in Cartesian coordinates. In particular cases, the results are in good agreement with the data obtained by other authors     

A shape memory alloy model by a second order phase transition

Motivated by the experimental results of the paper [1] and unlike the general theories of shape memory alloys (SMAs), in this paper we suggest for such materials a phase field model by a second order phase transition. So that, with this new system we obtain a simulation of phase dynamics very convenient to describe the natural behavior of these materials. The differential system is governed by the motion equation, the heat equation and the Ginzburg–Landau (GL) equation and by a constitutive law between the phase field, the temperature, the strain and the stress. The use of this new model is characterized by new potentials of the GL equation and by a new dependence on the temperature in the constitutive equation. Using this new model, we obtain simulations in better agreement with experimental data and respect to previous work [2].