Dynamic fluid–structure interaction of an elastic capsule in a viscous shear flow at moderate Reynolds number

The unsteady behavior of a 2-D circular elastic capsule was investigated in three viscous shear flows. An immersed boundary method (IBM) has been used to solve the dynamic fluid-structure interaction of the capsule. Computations were carried out in finite parameter ranges where the Reynolds number is Re=1-40 and the capillary number is Ca=0.0005-0.05, which is the ratio of the external viscous shear stress to the resistant elastic tensions of the membrane. For the simple shear flow, the effect of inertia on the transient behavior of the capsule was studied. For the pulsatile shear flow, two values of the peak fluid strain, T_f=1 and 5, were considered for the quasi-steady capsule mechanics. The capsule shows a cyclic structural response that includes subharmonics as the Reynolds number is elevated to 10 and 40. The capsule dynamic response includes a phase lag, which is a function of the capillary number, the Reynolds number, and the peak fluid strain. Finally, the capsule flowing in the Couette flow shows lateral migration due to the transient lift force, which is higher for lower Ca and higher Re. When capsules with diverse elasticity are dispersed along the velocity gradient, the capsule with a hard membrane experienced greater lift than the one with a soft membrane.     

In situ observations of unusual drop deformation and wobbling in simple shear flow

Deformation and wobbling of a liquid drop immersed in a liquid matrix were studied under mild shear conditions for various viscosity ratios. In situ visualization experiments were conducted on a homemade transparent Couette cell incorporated to the Paar Physica MCR500 shear rheometer. The effect of drop or matrix elasticity was examined and was found to play a major role in both deformation and wobbling processes. Experimental results were compared to Jackson and Tucker (J Rheol 47:659–682, 2003), Maffettone and Minale (J Non-Newton Fluid Mech 78:227–241, 1998) and Yu and Bousmina (J Rheol 47:1011–1039, 2003) ellipsoidal models. It was found that the agreement between the Newtonian models and the experimental results required an increase in the drop viscosity. Such increment in viscosity was found to scale with the first normal stress difference.     

Numerical study of single bubble motion in liquid metal exposed to a longitudinal magnetic field

The paper presents numerical simulations modeling the ascent of an argon bubble in liquid metal with and without an external magnetic field. The governing equations for the fluid and the electric potential are discretized on a uniform Cartesian grid and the bubble is represented with a highly efficient immersed boundary method. The simulations performed were conducted matching experiments under the same conditions so that sound validation is possible. The three-dimensional trajectory of the bubble is analyzed quantitatively and related to the flow structures in the wake. Indeed, the substantial impact of the magnetic field in the bubble trajectory results from its influence on the wake. Quantitative data describing the selective damping of vortex structures are provided and discussed. As a result of applying a longitudinal field, the time-averaged bubble rise velocity increases for large bubbles, it reaches a maximum and then decreases when further increasing the magnetic interaction parameter. For small bubbles, the time-averaged bubble rise velocity decreases when increasing the magnetic field. The bubble Strouhal number as a dimensionless frequency is reduced with the application of a magnetic field for all bubbles considered and the zig–zag trajectory of the bubble becomes more rectilinear.     

The orientational behavior of multiwall carbon nanotubes in polycarbonate in simple shear flow

This paper describes the changes in the orientation of multiwall carbon nanotubes (MWCNT) in polycarbonate as determined by transient and oscillatory shear rheology. It is well known from rheological studies on composites with macroscopic fibers that the overshoot in transient shear viscosity is caused by the change in orientation distribution of these fibers. This study shows that although an overshoot in transient shear viscosity of MWCNT/polycarbonate is measured at shear rates as low as 0.1 s^− 1, the MWCNT network is disturbed only at considerably higher shear rates. Scanning electron microscopy micrographs and oscillatory shear show that MWCNT in thermoplastic composites will only be oriented at high shear rates. Simultaneous measurements of the electrical conductivity during rheological start-up shear and oscillatory measurements show large differences between electrical and mechanical relaxation behaviors. The viscosity of the composite seems to depend strongly on the MWCNT network density, whereas the proximity of the tubes at the network points seems to determine the electrical properties of the MWCNT composite.     

Rolling/sliding of a particle on a flat wall in a linear shear flow at finite Re

Recently Lee and Balachandar proposed analytically-based expressions for drag and lift coefficients for a spherical particle moving on a flat wall in a linear shear flow at finite Reynolds number. In order to evaluate the accuracy of these expressions, we have conducted direct numerical simulations of a rolling particle for shear Reynolds number up to 100. We assume that the particle rolls on a horizontal flat wall with a small gap separating the particle from the wall (L = 0.505) and thus avoiding the logarithmic singularity. The influence of the shear Reynolds number and the translational velocity of the particle on the hydrodynamic forces of the particle was investigated under both transient and the final drag-free and torque-free steady state. It is observed that the quasi-steady drag and lift expressions of Lee and Balachandar provide good approximation for the terminal state of the particle motion ranging from perfect sliding to perfect rolling. With regards to transient particle motion in a wall-bounded shear flow it is observed that the above validated quasi-steady drag and lift forces must be supplemented with appropriate wall-corrected added-mass and history forces in order to accurately predict the time-dependent approach to the terminal steady state. Quantitative comparison with the actual particle motion computed in the numerical simulations shows that the theoretical models quite effective in predicting rolling/sliding motion of a particle in a wall-bounded shear flow at moderate Re.     

Evaluation of subgrid-scale models in large-eddy simulation of flow past a two-dimensional block

Large-eddy simulations of flow past a two-dimensional (2D) block were performed to evaluate four subgrid-scale (SGS) models: (i) the traditional Smagorinsky model, (ii) the Lagrangian dynamic model, (iii) the Lagrangian scale-dependent dynamic model, and (iv) the modulated gradient model. An immersed boundary method was employed to simulate the 2D block boundaries on a uniform Cartesian grid. The sensitivity of the simulation results to grid refinement was investigated by using four different grid resolutions. The velocity streamlines and the vertical profiles of the mean velocities and variances were compared with experimental results. The modulated gradient model shows the best overall agreement with the experimental results among the four SGS models. In particular, the flow recirculation, the reattachment position and the vertical profiles are accurately reproduced with a relative coarse grid resolution of (N_x × N_y × N_z=) 160 × 40 × 160 (n_x × n_z = 13 × 16 covering the block). Besides the modulated gradient model, the Lagrangian scale-dependent dynamic model is also able to give reasonable prediction of the flow statistics with some discrepancies compared with the experimental results. Relatively poor performance by the Lagrangian dynamic model and the Smagorinsky model is observed, with simulated recirculating patterns that differ from the measured ones. Analysis of the turbulence kinetic energy (TKE) budget in this flow shows evidence of a strong production of TKE in the shear layer that forms as the flow is deflected around the block.     

Effect of confinement on droplet deformation in shear flow

The dynamics of a single droplet under shear flow between two parallel plates is investigated by using the immersed boundary method. The immersed boundary method is appropriate for simulating the drop-ambient fluid interface. We apply a volume-conserving method using the normal vector of the surface to prevent mass loss of the droplet. In addition, we present a surface remeshing algorithm to cope with the distortion of droplet interface points caused by the shear flow. This mesh quality improvement in conjunction with the volume-conserving algorithm is particularly essential and critical for long time evolutions. We study the effect of wall confinement on the droplet dynamics. Numerical simulations show good agreement with previous experimental results and theoretical models.     

Rheological behavior of a dilute suspension of relatively coarse deformable particles in a simple shear flow

Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 55–61, January–February, 1990.     

双柱状颗粒在线性剪切流场中的动力学分析

为研究柱状颗粒在线性剪切流场中的运动状态和受力情况,本文以颗粒长径比为2,颗粒之间的初始距离ΔS_Py=4D为例,基于直接力浸入边界法数值模拟了双柱状颗粒在三维线性剪切流场中的运动过程。根据模拟结果分析了柱状颗粒周围流场参数分布,在考虑壁面对颗粒的影响和颗粒之间相互影响的条件下,研究了颗粒的受力和运动的变化,探索了流体曳力导致柱状颗粒迁移和转动的规律。研究结果表明,双柱状颗粒在线性剪切流场中易向速度大的流体区域运动;前后两颗粒运动状态和轨迹不同,颗粒之间距离较近时,曳力会产生较大的波动;只有当颗粒在壁面附近时,滞后颗粒才能追上领先颗粒,两颗粒发生牵引、翻滚和分离过程。     

The flexibility effect of a plate according to various angles of attack in a free-stream

This paper presents a computational fluid–structure interaction analysis for a flexible plate in a free-stream to investigate the effects of flexibility and angle of attack on force generation. A Lattice Boltzmann Method with an immersed boundary technique using a direct forcing scheme model of the fluid is coupled to a finite element model with rectangular bending elements. We investigated the effects of various angles of attack of a flexible plate fixed at one of the end edges in a free-stream at a Reynolds number of 5000, which represents the wing flapping condition of insects and small birds in nature. The lift of the flexible plate is maintained at the large angle of attack, whereas the rigid plate shows the largest lift at angles of attack around 30–40° and then drastic reductions in the lift at the large angle of attack. If we consider the efficiency as the lift divided by the drag, the flexible plate shows better efficiency at angles of attack greater than 30° compared to the rigid plate. The better performance of the flexible plate at large angles of attack comes from the deformation of the plate, which produces an interaction between the trailing edge vortex and the short edge vortex. The horseshoe-shaped vortex produced by a large vortex interaction at the trailing edge side has an important role in increasing the lift, and the small projection area due to the deformation reduces the drag. Furthermore, we investigate the role of flexibility on the lift and the drag force of the rectangular plate in a free-stream as the Reynolds number increases. Whenever a large vortex interaction at the trailing edge side is shown, the efficiency of the rectangular plate is improved. Especially, the flexible plate shows better efficiency as the Reynolds number increases regardless of the angle of attack.     

Numerical investigation of the influence of shear band localization on the resonant behavior in the VHCF regime

The monitored resonant behavior of fatigue specimens of metastable austenitic stainless steel (AISI304) is correlated with its damage accumulation in the very high cycle fatigue (VHCF) regime. The resonant behavior is studied experimentally and shows a distinct transient characteristic. Microscopic examinations indicate that during VHCF a localized plastic deformation in shear bands arises on the specimen surface. Hence, this work focuses on the effect of damage accumulation in shear bands on the resonant behavior of AISI304 in the VHCF regime. A microstructural simulation model is proposed that takes into account specific mechanisms in shear bands proven by experimental results. The simulation model is solved numerically using the two-dimensional boundary element method and the resonant behavior is characterized by evaluating the force-displacement hysteresis loop. Simulation of shear bands agrees well with microscopic examinations and plastic deformation in shear bands influences the transient characteristic of the resonant behavior.     

Lattice Boltzmann modeling of transport phenomena in fuel cells and flow batteries

Fuel cells and flow batteries are promising technologies to address climate change and air pollution problems. An understanding of the complex multiscale and multiphysics transport phenomena occurring in these electrochemical systems requires powerful numerical tools. Over the past decades, the lattice Boltzmann(LB) method has attracted broad interest in the computational fluid dynamics and the numerical heat transfer communities, primarily due to its kinetic nature making it appropriate for modeling complex multiphase transport phenomena. More importantly, the LB method fits well with parallel computing due to its locality feature, which is required for large-scale engineering applications. In this article, we review the LB method for gas–liquid two-phase flows, coupled fluid flow and mass transport in porous media, and particulate flows. Examples of applications are provided in fuel cells and flow batteries. Further developments of the LB method are also outlined.     

Chaotic dynamics of periodically forced spheroids in simple shear flow with potential application to particle separation

In this paper, we consider the technologically important problem of periodically forced spheroids in simple shear flow and demonstrate the existence of chaotic parametric regimes. The approach used by Strand (1989) (for the Strong Brownian limit) is inappropriate in the chaotic regimes corresponding to the weak Brownian limit. Our results also indicate a strong dependence of the solutions obtained on the aspect ratio of the spheroids. This strong dependence on the aspect ratio may be utilized to separate particles from a suspension of particles having different shapes but similar sizes.     

Experimental and numerical study of the rotation and the erosion of fillers suspended in viscoelastic fluids under simple shear flow

When a porous agglomerate immersed in a fluid is submitted to a shear flow, hydrodynamic stresses acting on its surface may cause a size reduction if they exceed the cohesive stress of the agglomerate. The aggregates forming the agglomerate are slowly removed from the agglomerate surface. Such a behaviour is known when the suspending fluid is Newtonian but unknown if the fluid is viscoelastic. By using rheo-optical tools, model fluids, carbon black agglomerates and particles of various shapes, we found that the particles had a rotational motion around the vorticity axis with a period which is independent on shape (flat particles not considered), but which is exponentially increasing with the elasticity of the medium expressed by the Weissenberg number (We). Spherical particles are always rotating for We up to 2.6 (largest investigated We in this study) but elongated particles stop rotating for We>0.9 while orienting along the flow direction. Erosion is strongly reduced by elasticity. Since finite element numerical simulation shows that elasticity increases the local stress around a particle, the origin of the erosion reduction is interpreted as an increase of cohesiveness of the porous agglomerate due to the infiltration of a viscoelastic fluid.     

A calculating method of shock wave oscillating frequency due to turbulent shear layer fluctuations in supersonic flow

One of the more severe fluctuating pressure environments encountered in supersonic orhypersonic flows is the shock wave oscillation driven by interaction of a shock wave withboundary layer.The high intensity oscillating shock wave may induce structure resonanceof a high speed vehicle.The research for the shock oscillation used to adopt empirical orsemiempirical methods because the phenomenon is very complex.In this paper atheoretical solution on shock oscillating frequency due to turbulent shear layer fluctuationshas been obtained from basic conservation equations.Moreover,we have attained theregularity of the frequency of oscillating shock varying with incoming flow Mach numbersM_∞and turning angleθ.The calculating results indicate excellent agreement withmeasurements.This paper has supplied a valuable analytical method to study aeroelasticproblems produced by shock wave oscillation.     

Dynamic behavior of two collinear anti-plane shear cracks in a functionally graded layer bonded to dissimilar half planes

In recent years, the functionally graded materials (FGMs) have been widely applied in extremely high temperate environment. In this paper, the dynamic behavior of two collinear cracks in FGM layer bonded to dissimilar half planes under anti-plane shear waves is studied by the Schmidt method. By using the Fourier transform technique, the present problem can be solved with a dual integral equation. These equations are solved using the Schmidt method. The present method is used to illustrate the fundamental behavior of the interacting cracks in FGMs under dynamic loading. Furthermore, the effects of the geometry of the interacting cracks, the shear stress wave velocity of the materials and the frequency of the incident wave on the Dynamic Stress Intensity Factor are investigated.     

Simulations of chemical transport and reaction in a suspension of cells I: an augmented forcing point method for the stationary case

A novel augmented forcing point method is presented to solve the problem of chemical transport in the fluid outside of a collection of suspended cells coupled with chemical reactions on the surfaces of the cells. In this method, the chemical concentrations and the forcing function values are determined simultaneously from an augmented system of equations. The method is more stable and accurate than predictor‐corrector‐type forcing point methods, yields the same solution as a corresponding ghost cell method with much less computational cost, and provides solutions that are pointwise second‐order accurate in space and time, even for closely‐spaced cells. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical analysis of the developing fluid flow in a circular duct rotating steadily about a parallel axis

A numerical analysis of the flow pattern in the inlet region of a circular pipe rotating steadily about an axis parallel to its own is presented. Both finite cell and finite element methods are used to analyse the problem and they give qualitatively similar results which show that a swirling fluid motion is induced in the pipe inlet region. The analyses show that the direction of swirl is opposite to that of the pipe rotation when viewed along the flow axis and that its magnitude depends on the speed of pipe rotation and throughflow Reynolds number. Neither numerical analysis predicts the marked upturn in friction factor (or pressure drop) which has been observed experimentally. However, a dependence on the pipe inlet boundary conditions is demonstrated.