Sliding plate rheometry of planar oriented concentrated fiber suspension

The rheology of concentrated planar fiber suspensions is investigated. A new experimental technique for fiber suspensions based on a sliding plate rheometer incorporating a shear stress transducer is developed. It is shown that this instrument works well for the tested material systems. The rheological behavior in steady shear is subsequently investigated. The results can be largely explained by a combination of frictional and hydrodynamic interaction. Despite this evidence of friction no yield stress could be detected for the investigated shear rates. It was also found that the fiber aspect ratio did not influence the steady shear viscosity.     

Fiber orientation kinetics of a concentrated short glass fiber suspension in startup of simple shear flow

The common approach for simulating the evolution of fiber orientation during flow in concentrated suspensions is to use an empirically modified form of Jeffery's equation referred to as the Folgar–Tucker (F-T) model. Direct measurements of fiber orientation were performed in the startup of shear flow for a 30 wt% short glass fiber-filled polybutylene terephthalate (PBT-30); a matrix that behaves similar to a Newtonian fluid. Comparison between predictions based on the F-T model and the experimental fiber orientation show that the model over predicts the rate of fiber reorientation. Rheological measurements of the stress growth functions show that the stress overshoot phenomenon approaches a steady state at a similar strain as the fiber microstructure, at roughly 50 units. However, fiber orientation measurements suggest that a steady state is not reached as the fiber orientation continues to slowly evolve, even up to 200 strain units. The addition of a “slip” parameter to the F-T model improved the model predictions of the fiber orientation and rheological stress growth functions.     

Comparative numerical study of single and two-phase models of nanofluid heat transfer in wavy channel

The main purpose of this study is to survey numerically comparison of two- phase and single phase of heat transfer and flow field of copper-water nanofluid in a wavy channel. The computational fluid dynamics （CFD） prediction is used for heat transfer and flow prediction of the single phase and three different two-phase models （mixture, volume of fluid （VOF）, and Eulerian）. The heat transfer coefficient, temperature, and velocity distributions are investigated. The results show that the differences between the temperature fie].d in the single phase and two-phase models are greater than those in the hydrodynamic tleld. Also, it is found that the heat transfer coefficient predicted by the single phase model is enhanced by increasing the volume fraction of nanoparticles for all Reynolds numbers; while for the two-phase models, when the Reynolds number is low, increasing the volume fraction of nanoparticles will enhance the heat transfer coefficient in the front and the middle of the wavy channel, but gradually decrease along the wavy channel.     

An experimental study of flow-induced fiber orientation and concentration distributions in a concentrated suspension flow through a slit channel containing a cylinder

Flow-induced fiber orientation and concentration distributions were measured in a concentrated fiber suspension (CFS) and a dilute one (DFS). The channel has a thin slit geometry containing a circular cylinder. In the previous work, many researchers have qualitatively studied fiber orientation and concentration distributions in injection-molded products of fiber-reinforced plastics. In the present work, however, they are quantitatively estimated by direct observation of fibers in the concentrated suspension flow. For the CFS, some fibers rotate in an expansion part between the channel wall and the circular cylinder, and the fiber orientation becomes almost random state. On the other hand, fibers are perfectly aligned along the flow direction owing to the elongational flow near the centerline downstream of the cylinder. The fiber concentration has a flat distribution except near the channel wall and the centerline. For the DFS a minimum in the fiber concentration distribution was clearly observed on the centerline, and two peaks beside the centerline and near the channel wall. This characteristic distribution is caused by the fiber-wall and fiber-cylinder interactions. It is found that the obstacle such as the circular cylinder in the channel significantly affects the fiber orientation downstream of the obstacle for the CFD, while it affects the fiber concentration distribution for the DFS.     

裂纹扩展寿命模型的对比研究

目前，已提出了许多种随机载荷谱裂纹扩展模型，但还很少有人对这些模型进行对比研究，并对工程应用提出建议。本文使用战斗机综合任务谱，对标准试件，用不同的裂纹扩展模型进行实际测算，并与试验值比较。通过研究得出了一些有益的结论这对实际工程应用十分有益。     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

A study of floc breakup and formation in flowing concentrated fiber suspensions

An experimental study was carried out to examine the effects of fiber concentration and average flow velocity on the breakup of fiber flocs and the subsequent reflocculation process downstream of a converging/diverging flow section used as a turbulence generator. Flocculation measurements were made using fiber optic probes and are presented in terms of floc-size distributions. These were complemented with pressure fluctuation measurements used to indicate the state of fluidization of the suspension.     

Nonlinear rheological behavior of a concentrated spherical silica suspension

Linear and nonlinear viscoelastic properties were examined for a 50 wt% suspension of spherical silica particles (with radius of 40 nm) in a viscous medium, 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The effective volume fraction of the particles evaluated from zero-shear viscosities of the suspension and medium was 0.53. At a quiescent state the particles had a liquid-like, isotropic spatial distribution in the medium. Dynamic moduli G^* obtained for small oscillatory strain (in the linear viscoelastic regime) exhibited a relaxation process that reflected the equilibrium Brownian motion of those particles. In the stress relaxation experiments, the linear relaxation modulus G(t) was obtained for small step strain (0.2) while the nonlinear relaxation modulus G(t, ) characterizing strong stress damping behavior was obtained for large (>0.2). G(t, ) obeyed the time-strain separability at long time scales, and the damping function h() (–G(t, )/G(t)) was determined. Steady flow measurements revealed shear-thinning of the steady state viscosity () for small shear rates (< ^–1; = linear viscoelastic relaxation time) and shear-thickening for larger (>^–1). Corresponding changes were observed also for the viscosity growth and decay functions on start up and cessation of flow, + (t, ) and _– (t, ). In the shear-thinning regime, the and dependence of ₊(t,) and _–(t,) as well as the dependence of () were well described by a BKZ-type constitutive equation using the G(t) and h() data. On the other hand, this equation completely failed in describing the behavior in the shear-thickening regime. These applicabilities of the BKZ equation were utilized to discuss the shearthinning and shear-thickening mechanisms in relation to shear effects on the structure (spatial distribution) and motion of the suspended particles.Dedicated to the memory of Prof. Dale S. Parson     

Nonlinear rheology of a concentrated spherical silica suspension:

Time-dependent nonlinear flow behavior was investigated for a model hard-sphere suspension, a 50 wt% suspension of spherical silica particles (radius = 40 nm; effective volume fraction = 0.53) in a 2.27/1 (wt/wt) ethylene glycol/glycerol mixture. The suspension had two stress components, the Brownian stress ^B and the hydrodynamic stress ^H After start-up of flow at various shear rates , the viscosity growth function ₊ (t, ) was measured with time t until it reached the steady state. The viscosity decay function _– (t, ) was measured after cessation of flow from the steady as well as transient states. At low where the steady state viscosity ( ) exhibited the shear-thinning, the _– (t, ) and ⁺ (t, ) data were quantitatively described with a BKZ constitutive equation utilizing data for nonlinear relaxation moduli G (t, ). This result enabled us to attribute the thinning behavior to the decrease of the Brownian contribution ^B = ^B / (considered in the BKZ equation through damping of G (t, )). On the other hand, at high where ( ) exhibited the thickening, the BKZ prediction largely deviated from the ₊ (t, ) and ₊ (t, ) data, the latter obtained after cessation of steady flow. This result suggested that the thickening was due to an enhancement of the hydrodynamic contribution ^H = ^H / (not considered in the BKZ equation). However, when the flow was stopped at the transient state and only a small strain (<0.2) was applied, ^H was hardly enhanced and the _– (t, ) data agreed with the BKZ prediction. Correspondingly, the onset of thickening of ₊ (t, ) was characterized with a -insensitive strain ( 0.2). On the basis of these results, the enhancement of ^H (thickening mechanism) was related to dynamic clustering of the particles that took place only when the strain applied through the fast flow was larger than a characteristic strain necessary for close approach/collision of the particles.     

A numerical simulation on diffuser‐nozzle based piezoelectric micropumps with two different numerical models

This study has been conducted to compare two different numerical models for the evaluation of the performance characteristics of a diffuser‐nozzle based piezoelectric micropump. Here, the transient displacements of the membrane and the flow characteristics in the piezoelectric micropump have been closely investigated with the FSI model and the prescribed deformation model for two different frequencies. It has been found that the behaviour of the membrane computed with the FSI model is not in accordance with that with the prescribed deformation model, and that the net flow rate with the FSI model is larger than that with the prescribed deformation model. Therefore, the choice of numerical model is very important in conducting numerical analysis for piezoelectric micropumps. The results of this study can be utilized as basic data for the design and analysis of piezoelectric micropumps. Copyright © 2006 John Wiley & Sons, Ltd.     

两种颗粒湍流扩散模型数值模拟气液两相流泄漏扩散的比较

采用相间耦合的欧拉-拉格朗日方法,模拟了装有液化气（丙烷）的容器出现小孔或裂缝时,发生泄漏后的气液两相扩散过程。分别应用随机轨道模型及颗粒群模型来考察湍流对液滴扩散的影响,并与确定性轨道进行了比较。给出了数学物理模型,计算结果与实验数据进行了对比。结果表明：采用随机轨道模型能较好地描述液滴的湍流扩散,适用于有液相蒸发的两相流扩散问题。     

On the viscosity of a concentrated suspension of solid spheres

Explanations of the very high viscosities of concentrated suspensions of spheres based on the dissipation in squeezing flow between particles pairs are shown to be in error. The dissipation in pair interactions is always of the order of that generated by shearing motions, and this dissipation is of too low an order in concentration for single pair interactions to explain the observed viscosities.     

Physical characterisation and yield stress of a concentrated Miscanthus suspension

We present yield stress measurements of a concentrated suspension of Miscanthus x giganteus, a biomass substrate for the production of biofuels and chemicals. The particle porosity and the apparent solid packing fraction are quantified, and a simple relationship between volume and mass fractions proposed. The relationship between yield stress and solid concentration was measured with the maximum torque method of Dzuy and Boger (J Rheol 27: 321–349, 1983) for concentrations between 11 and 15 %, just below maximum packing. In the range of concentrations considered, the measured dependence of the yield stress on concentration appears to be significantly stronger than that reported for corn stover, a plant for which reliable data is available: the apparent power law exponent is found to be larger than 10, while the corresponding exponent reported for corn is in the range 4–6. We argue that differences in particle rigidity between corn and Miscanthus may play a role in explaining this result.     

Effective viscosity of a concentrated suspension of composite porous spherical particles

Meccanica - We present an analytical study of the effective viscosity of concentrated suspension of porous spherical particles with a rigid core, under the creeping flow conditions. It is assumed...     

Comparative study of two corrective gradient models in the simulation of multiphase flows using moving particle semi-implicit method

In this study, two corrective gradient models (CGMs) are compared in the simulation of multiphase flows. Linear consistency of the gradient model of moving particle semi-implicit (MPS) method has been recovered by introducing corrective matrix. However, it is found that particles tend to disperse along the streamline while using the CGM proposed in a previous study. Particle shifting (PS) schemes are necessary to reduce the irregularity of particle distribution to stabilize the calculation. To enhance the accuracy and stability, another CGM with dummy particle (CGMD) was proposed in our previous study. This enhanced CGM is featured by linear consistency and purely repulsive pressure gradient force. In this study, this enhanced CGM is modified and applied to multiphase flow simulation. Comparative study suggests that the modified CGM with PS scheme is capable of calculating various multiphase flows and predicting the interface evolution both clearly and accurately.     

The merging of two co-rotating vortices: a numerical study

The merging of two-dimensional co-rotating vortices is analysed through direct numerical simulations at large Reynolds numbers. It is shown how the Reynolds number affects each of the three phases that characterise this phenomenon. In the first phase, we examine the merging onset and focus on its definition. During the second rapid phase, the contributions of various flow regions upon the dynamics of a vortex are quantitatively studied. These regions are respectively the companion vortex, the filaments and an intermediate zone between vortices and filaments. The third phase is interpreted in terms of an advection diffusion process. Finally the final profile and circulation of the merged vortex is determined: the two thirds of the total circulation of the two initial vortices is contained in the newly formed vortex.     

A numerical study of constitutive models endowed with Pom-Pom molecular attributes

The branched polymer melts are modeled respectively in this investigation by the existing XPP and PTT–XPP models, along with the proposed S-MDCPP (Single/Simplified Modified Double Convected Pom-Pom) model developed on the basis of the existing MDCPP model. A pressure stabilized mass equation is formulated with the finite increment calculus (FIC) process to restrain and further eliminate spurious oscillations of pressure field due to the incompressibility of fluids. The discrete elastic viscous stress splitting (DEVSS) technique is employed, in order to retain an elliptic contribution in the weak form of the momentum equation. An inconsistent streamline-upwind (SU) method is applied to spatially discretize the constitutive equations. The mass, momentum conservation and constitutive equations are discretized and solved by the iterative stabilized fractional step algorithm along with the Crank–Nicolson implicit difference scheme. Thus the finite elements with equal low-order interpolation approximations for velocity–pressure–stress variables can be devised to numerically simulate the viscoelastic contraction flows for branched LDPE melts. The influences of the three viscoelastic constitutive models and the branched arms at the end of the Pom-Pom molecule on the rheological behaviors occurring in this complex flow are discussed. The numerical results demonstrate that the proposed S-MDCPP model is capable of reproducing some properties similar to those predicted by the XPP model in high shear flow and, on the other hand, reproducing some properties similar to those predicted by the PTT–XPP model in high elongational flow. Furthermore, the proposed S-MDCPP model is capable of well identifying the macromolecule topological structures of branched polymer melts.     

金属靶板侵彻数值模拟对比研究

数值模拟是研究靶板侵彻问题的重要手段。为对相关软件的选取提供参考,采用当前广泛使用的3种有限元软件（LS-DYNA、ABAQUS、PAM-CRASH）对同一个靶板侵彻实验进行了数值模拟,并比较了各模拟软件的优缺点。研究结果表明:3种软件的模拟结果与实验结果均基本一致;对平头形弹体侵彻的模拟效果普遍优于对半球形弹体的模拟效果;在弹体剩余速度、冲塞块速度2个方面,3种有限元软件的数值模拟结果与实验结果大致吻合。其中,ABAQUS和PAM-CRASH对上述特征量的模拟结果更接近实验结果,平均相对误差普遍小于15%,且PAM-CRASH多数模拟结果误差都小于10%;而在弹体变形量上,3种软件的模拟结果均与实验结果相差较大。此外,LS-DYNA在模拟时的鲁棒性较好,模拟结果较稳定,PAM-CRASH对参数设置最敏感,ABAQUS在计算时间和效果上较平衡。