Coclite Alessandro, Coclite Giuseppe Maria, De Tommasi Domenico
School of Engineering, Università della Basilicata, 85100 Potenza, Italy.
Dipartimento di Meccanica, Matematica e Management, 70126 Politecnico di Bari, Italy.
Nanomaterials (Basel). 2020 Nov 3;10(11):2190. doi: 10.3390/nano10112190.
In this paper, a Multi Relaxation Time Lattice Boltzmann scheme is used to describe the evolution of a non-Newtonian fluid. Such method is coupled with an Immersed-Boundary technique for the transport of arbitrarily shaped objects navigating the flow. The no-slip boundary conditions on immersed bodies are imposed through a convenient forcing term accounting for the hydrodynamic force generated by the presence of immersed geometries added to momentum equation. Moreover, such forcing term accounts also for the force induced by the shear-dependent viscosity model characterizing the non-Newtonian behavior of the considered fluid. Firstly, the present model is validated against well-known benchmarks, namely the parabolic velocity profile obtained for the flow within two infinite laminae for five values of the viscosity model exponent, = 0.25, 0.50, 0.75, 1.0, and 1.5. Then, the flow within a squared lid-driven cavity for = 1000 and 5000 (being the Reynolds number) is computed as a function of for a shear-thinning ( < 1) fluid. Indeed, the local decrements in the viscosity field achieved in high-shear zones implies the increment in the local Reynolds number, thus moving the position of near-walls minima towards lateral walls. Moreover, the revolution under shear of neutrally buoyant plain elliptical capsules with different Aspect Ratio ( = 2 and 3) is analyzed for shear-thinning ( < 1), Newtonian ( = 1), and shear-thickening ( > 1) surrounding fluids. Interestingly, the power law by Huang et al. describing the revolution period of such capsules as a function of the Reynolds number and the existence of a critical value, c, after which the tumbling is inhibited in confirmed also for non-Newtonian fluids. Analogously, the equilibrium lateral position yeq of such neutrally buoyant capsules when transported in a plane-Couette flow is studied detailing the variation of yeq as a function of the Reynolds number as well as of the exponent .
在本文中,采用多松弛时间格子玻尔兹曼格式来描述非牛顿流体的演化。该方法与浸入边界技术相结合,用于任意形状物体在流场中的输运。通过在动量方程中添加一个方便的强迫项来施加浸入物体上的无滑移边界条件,该项考虑了浸入几何形状所产生的流体动力。此外,该强迫项还考虑了由表征所考虑流体非牛顿行为的剪切依赖粘度模型所引起的力。首先,针对著名的基准问题对当前模型进行验证,即对于粘度模型指数(\alpha)的五个值((\alpha = 0.25)、(0.50)、(0.75)、(1.0)和(1.5)),计算在两个无限平板间流动所得到的抛物线速度分布。然后,针对(\alpha = 1000)和(5000)((Re)为雷诺数),计算了方腔顶盖驱动流中剪切变稀((\alpha < 1))流体的流动,(\alpha)的函数。实际上,在高剪切区域实现的粘度场局部减小意味着局部雷诺数增加,从而使近壁最小值位置向侧壁移动。此外,针对剪切变稀((\alpha < 1))、牛顿((\alpha = 1))和剪切增稠((\alpha > 1))的周围流体,分析了不同纵横比((AR = 2)和(3))的中性浮力平面椭圆胶囊在剪切作用下的旋转。有趣的是,Huang等人描述此类胶囊旋转周期作为雷诺数函数的幂律以及存在一个临界值(Re_c)(超过该值翻滚被抑制),对于非牛顿流体也得到了证实。类似地,研究了此类中性浮力胶囊在平面库埃特流中输运时的平衡横向位置(y_{eq}),详细说明了(y_{eq})作为雷诺数以及指数(\alpha)的函数变化。