Journal of Non-Newtonian Fluid Mechanics最新文献

A lattice Boltzmann flux solver with log-conformation representation for the simulations of viscoelastic flows at high Weissenberg numbers

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-26 DOI: 10.1016/j.jnnfm.2024.105351

Hua Zhang , Chang Shu , Lian-Ping Wang , Yaguang Liu

{"title":"A lattice Boltzmann flux solver with log-conformation representation for the simulations of viscoelastic flows at high Weissenberg numbers","authors":"Hua Zhang , Chang Shu , Lian-Ping Wang , Yaguang Liu","doi":"10.1016/j.jnnfm.2024.105351","DOIUrl":"10.1016/j.jnnfm.2024.105351","url":null,"abstract":"<div><div>In this work, a viscoelastic lattice Boltzmann flux solver (VLBFS) with log-conformation representation is proposed for simulating the incompressible flows of a viscoelastic fluid at high Weissenberg number conditions. Compared with the original lattice Boltzmann flux solver (LBFS), the present method has two main new features. First, the method solves the polymer constitutive equations with log-conformation representation. Second, an upwind-biased scheme is incorporated in the interpolation when performing flux reconstructions at the cell interface. With the aid of these two treatments, the numerical stability of VLBFS is significantly improved, making it capable of solving high Weissenberg number problems (HWNP). Compared with using the lattice Boltzmann method (LBM) to solve the viscoelastic fluid flow, VLBFS inherits the advantages of LBFS, such as flexible mesh generation, decoupling of the grid spacing and time interval, and low memory requirement. VLBFS can also precisely recover the macroscopic constitutive equation. The present method has been critically validated using three benchmark cases, namely, the plane Poiseuille flow, lid-driven cavity flow, and 4:1 abrupt planar contraction flow. The numerical results fully demonstrate the solver’s powerful ability in simulating HWNP.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105351"},"PeriodicalIF":2.7,"publicationDate":"2024-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142744985","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Analysis of the shear thickening behavior of a fumed silica suspension using QL-LAOS approach

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-26 DOI: 10.1016/j.jnnfm.2024.105355

E. Fernández-Díaz , F.J. Rubio-Hernández , J.F. Velázquez-Navarro

{"title":"Analysis of the shear thickening behavior of a fumed silica suspension using QL-LAOS approach","authors":"E. Fernández-Díaz , F.J. Rubio-Hernández , J.F. Velázquez-Navarro","doi":"10.1016/j.jnnfm.2024.105355","DOIUrl":"10.1016/j.jnnfm.2024.105355","url":null,"abstract":"<div><div>The feasibility of applying quasi-linear large amplitude oscillatory shear (QL-LAOS) approach to a shear thickening (ST) fumed silica suspension was tested. While the characteristic time has been used as the parameter for the original QL-LAOS analysis of shear thinning fluids, we obtained that a description based upon increasing stiffness is more appropriate for ST fumed silica suspensions. Very low <span><math><mrow><mo>(</mo><mrow><mo>≤</mo><mn>1.5</mn></mrow><mo>)</mo></mrow></math></span> third to first harmonics ratio were obtained indicating the need of alternative criteria to identify QL-LAOS behavior in ST suspensions. Consequently, a method based upon the best fit of an ellipse to the experimental Lissajous-Bowditch curves was proposed. Compliances were obtained from areas of viscous and elastic fitted ellipses. The dependence of the material functions obtained by using a Jeffrey´s mechanical viscoelastic framework with angular frequency supports the idea of ST microstructure evolves by increasing with shear the number of small hydroclusters.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105355"},"PeriodicalIF":2.7,"publicationDate":"2024-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142744986","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

The influence of thixotropy on bubble growth in thixotropic yield stress fluids: Insights from numerical simulations

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-19 DOI: 10.1016/j.jnnfm.2024.105353

Yixuan Hou , Zhao Jin , Xinzhe Que , Yongchao Zhou , Yiping Zhang

{"title":"The influence of thixotropy on bubble growth in thixotropic yield stress fluids: Insights from numerical simulations","authors":"Yixuan Hou , Zhao Jin , Xinzhe Que , Yongchao Zhou , Yiping Zhang","doi":"10.1016/j.jnnfm.2024.105353","DOIUrl":"10.1016/j.jnnfm.2024.105353","url":null,"abstract":"<div><div>Bubble behaviors in structured fluids are of great interests in industrial applications, while there is currently a lack of understanding regarding the effect of thixotropic microstructure on the bubble formation process. To this end, this study explores the influence of thixotropy on bubble growth in thixotropic yield stress fluids by numerical simulations using the Arbitrary Lagrangian-Eulerian (ALE) method. The numerical results reveal that, with the increase in the thixotropy number, the bubbles at detachment transform from inverted conical to spherical shapes at lower gas flow rates, and from spindle to conical shapes at higher gas flow rates, along with the decreased detachment volume and time. It is also found that the effect of gas flow rate varies with different thixotropy numbers. The flow field of the structured fluid reveals that the increases in gas flow rate primarily promote the structural destruction near the bubble tip, while the increase in thixotropy number facilitate the fluid flow around the bubble, with the significant reduction of the low-shear zones and expansion of the yielded zones near the equatorial plane. As a result, modulating the fluid flow with thixotropy number mainly influences the hydrodynamic pressure on the bubble. Based on a force balance model, the forces acting on the bubble are then calculated by integrating the stress on the interface, and it is found that thixotropy number controls the bubble detachment state with the drag effect. Accordingly, the mechanisms governing the influence of thixotropy on drag effect are discussed considering the flow field characteristics and the correlations of drag coefficients. This work helps to deepen the understanding of the bubble behaviors in structured fluids.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105353"},"PeriodicalIF":2.7,"publicationDate":"2024-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142744983","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Suppression and augmentation in vortex shedding frequency due to fluid elasticity

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-19 DOI: 10.1016/j.jnnfm.2024.105348

A. Chauhan, Sana Raffi, C. Sasmal

{"title":"Suppression and augmentation in vortex shedding frequency due to fluid elasticity","authors":"A. Chauhan, Sana Raffi, C. Sasmal","doi":"10.1016/j.jnnfm.2024.105348","DOIUrl":"10.1016/j.jnnfm.2024.105348","url":null,"abstract":"<div><div>Several previous experimental and numerical studies have demonstrated the suppression of vortex shedding frequency from bluff bodies, such as circular cylinders, due to fluid elasticity induced by adding solid polymer additives to a solvent like water, even in parts per million (ppm) quantities. However, this study reveals a more complex relationship between the two using extensive two-dimensional (2D) direct numerical simulations (DNS) of flows past a circular cylinder at a fixed Reynolds number of 100. Our findings show that the vortex shedding frequency initially decreases with increasing Weissenberg number (a measure of fluid elasticity), reaches a minimum at a critical Weissenberg number, and then increases with further increments in the Weissenberg number. The same non-monotonic trend is also observed in the temporal variation of the spanwise velocity component fluctuation within the flow domain. This study aims to elucidate the reasons behind these non-monotonic trends in vortex shedding frequency and velocity component fluctuations as functions of the Weissenberg number. Our detailed analysis attributes these trends to significant alterations in the vortex-shedding mechanism as fluid elasticity increases due to the appearance of inertio-elastic instability at higher Weissenberg numbers. Our findings also align with limited experimental observations of similar unexpected behaviors in viscoelastic fluids, providing new insights into the underlying mechanisms. Moreover, the study highlights that shear-thinning behavior in viscoelastic fluids counteracts these non-monotonic trends, instead promoting a monotonic increase in vortex shedding frequency with the Weissenberg number. Finally, the 2D simulation results show both qualitative and quantitative agreement with limited three-dimensional (3D) simulations conducted at higher Weissenberg numbers where the flow may transit from 2D to 3D due to the appearance of inertio-elastic instability.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105348"},"PeriodicalIF":2.7,"publicationDate":"2024-11-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142744982","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Viscoelastic model hierarchy for fiber melt spinning of semi-crystalline polymers

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-16 DOI: 10.1016/j.jnnfm.2024.105349

Manuel Ettmüller , Walter Arne , Nicole Marheineke , Raimund Wegener

{"title":"Viscoelastic model hierarchy for fiber melt spinning of semi-crystalline polymers","authors":"Manuel Ettmüller , Walter Arne , Nicole Marheineke , Raimund Wegener","doi":"10.1016/j.jnnfm.2024.105349","DOIUrl":"10.1016/j.jnnfm.2024.105349","url":null,"abstract":"<div><div>In the fiber melt spinning of semi-crystalline polymers, the degree of crystallization can be non-homogeneous over the cross-section of the fiber, affecting the properties of the end product. For simulation-based process design, the question arises as to which fiber quantities and hence model equations must be resolved in radial direction to capture all practically relevant effects and at the same time imply a model that can be computed with reasonable effort. In this paper, we present a hierarchy of viscoelastic two-phase fiber models ranging from a complex, fully resolved and highly expensive three-dimensional description to a cross-sectionally averaged, cheap-to-evaluate one-dimensional model. In particular, we propose a novel stress-averaged one-two-dimensional fiber model, which circumvents additional assumptions on the inlet profiles needed in the established stress-resolved fiber model by Doufas et al. (2001). Simulation results demonstrate the performance and application regime of the dimensionally reduced models. The novel stress-averaged variant provides fast and reliable results, especially in the regime of low flow-enhanced crystallization.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105349"},"PeriodicalIF":2.7,"publicationDate":"2024-11-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142744981","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Enhancing polymer fiber orientation with 3D-printed shell–core structures

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-15 DOI: 10.1016/j.jnnfm.2024.105352

Mengfan Lou, Zhenyu Ouyang

{"title":"Enhancing polymer fiber orientation with 3D-printed shell–core structures","authors":"Mengfan Lou, Zhenyu Ouyang","doi":"10.1016/j.jnnfm.2024.105352","DOIUrl":"10.1016/j.jnnfm.2024.105352","url":null,"abstract":"<div><div>This study numerically investigates the flow dynamics and fiber orientation in enhanced polymer shell–core structures during 3D printing using the smoothed particle hydrodynamics method. A microstructure-based fiber suspension model, coupled with a shear-thinning viscosity model, is employed to assess the effects of varying material viscosities, fiber aspect ratios and volume fractions, and substrate speeds on fiber orientation distribution. Our results demonstrate that fiber orientation in deposited layers is particularly sensitive to variations in the viscosity of the shell material, while changes in the core material viscosity have a smaller and more localized impact on fiber orientation. Furthermore, fiber orientation is strongly influenced by the product of the fiber aspect ratio <em>α<sub>r</sub></em> and volume fraction <em>ϕ</em>; as <em>α<sub>r</sub>ϕ</em> increases, fibers tend to align with the flow direction of materials. However, adjusting the fiber aspect ratio and volume fraction while maintaining the same <em>α<sub>r</sub>ϕ</em> results in slight changes to fiber alignment. Additionally, lower substrate speeds cause upstream material accumulation, increasing deposition layer height and creating differences in fiber orientation between the core and shell regions. An appropriate increase in substrate speed can mitigate these effects.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"335 ","pages":"Article 105352"},"PeriodicalIF":2.7,"publicationDate":"2024-11-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142745105","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Pressure drop reduction due to coupling between shear-thinning fluid flow and a weakly deformable channel wall: A reciprocal theorem approach

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-11-14 DOI: 10.1016/j.jnnfm.2024.105347

Shrihari D. Pande, Ivan C. Christov

引用次数: 0

Flow enhancement produced by a pulsatile flow of shear-thinning fluids in circular and concentric annular tubes 剪切稀化流体在圆形和同心环形管中的脉动流动所产生的流动增强效应

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-10-28 DOI: 10.1016/j.jnnfm.2024.105346

J. Báez-Amador , R. Baños , J. Arcos , F. Méndez , O. Bautista

{"title":"Flow enhancement produced by a pulsatile flow of shear-thinning fluids in circular and concentric annular tubes","authors":"J. Báez-Amador , R. Baños , J. Arcos , F. Méndez , O. Bautista","doi":"10.1016/j.jnnfm.2024.105346","DOIUrl":"10.1016/j.jnnfm.2024.105346","url":null,"abstract":"<div><div>Although the analysis of the flow enhancement of non-Newtonian fluids produced by pulsatile flows through tubes is common in the literature, the case of Carreau fluid has not been analyzed, which is the aim of this work. This study determines the flow enhancement caused by the pulsatile fluid flow through <strong>(a)</strong> a circular tube and <strong>(b)</strong> a concentric annular tube. We show that the flow rate enhancement of the shear-thinning fluid is controlled by the Carreau number <span><math><mrow><mi>C</mi><mi>u</mi></mrow></math></span>, the Womersley number <span><math><mrow><mi>W</mi><mi>o</mi></mrow></math></span>, the fluid power-law index <span><math><mi>n</mi></math></span>, the ratio between the outer and inner radii <span><math><mi>κ</mi></math></span>, a parameter <span><math><mi>β</mi></math></span> that represents the ratio between the infinite and zero-shear viscosities, and the amplitude of the oscillatory signal <span><math><mi>ɛ</mi></math></span>. In both cases <strong>(a)</strong> and <strong>(b)</strong>, a numerical solution of the start-up of the hydrodynamic is evaluated. With the aid of the velocity solution, the volumetric flow rate is determined under periodic conditions after the initial transient has vanished. Then, the fractional increase in the mean flow rate <span><math><mi>I</mi></math></span> due to the pulsatile pressure gradient is calculated. Furthermore, an asymptotic solution for small, intermediate, and very large values of the Carreau number is performed to provide physical insight into flow enhancement.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"334 ","pages":"Article 105346"},"PeriodicalIF":2.7,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142571410","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

A fully-coupled implementation of the contravariant deformation tensor formulation for viscoelastic flows 粘弹性流动的反常变形张量公式的全耦合实现

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-10-28 DOI: 10.1016/j.jnnfm.2024.105345

Nick O. Jaensson, Martien A. Hulsen

{"title":"A fully-coupled implementation of the contravariant deformation tensor formulation for viscoelastic flows","authors":"Nick O. Jaensson, Martien A. Hulsen","doi":"10.1016/j.jnnfm.2024.105345","DOIUrl":"10.1016/j.jnnfm.2024.105345","url":null,"abstract":"<div><div>We present a fully-coupled numerical scheme for computing steady and time-dependent viscoelastic flows. The scheme relies on the contravariant deformation tensor formulation and uses a Newton–Raphson iteration to solve the non-linear system of equations. The contravariant reformulation allows for the computation and implementation of the analytical Jacobian relatively easily, especially compared to other reformulations such as the log-conformation. The contravariant deformation tensor rotates in steady state shearing flows, which is solved here by “resetting” it as a pre-processing step in the numerical scheme, rather than a post-processing step. We use the finite element method with standard stabilization techniques (SUPG and DEVSS-G) for the spatial discretization. The numerical scheme is tested in three viscoelastic flow problems which are studied in terms of stability and accuracy: planar Couette flow, 2D flow around a cylinder and 3D flow around a sphere. For all problems, quadratic convergence is observed in both the difference between iterations and the residuals during the Newton–Raphson procedure. Moreover, we observe that the residuals are several orders smaller than the difference between iterations. A distinct advantage of the numerical scheme presented here, is that it significantly relaxes the requirement on the time-step size in time-dependent problems, as compared to explicit or semi-implicit methods. Moreover, steady states can be efficiently computed if the initial guess in the Newton–Raphson iteration is close enough to the solution.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"334 ","pages":"Article 105345"},"PeriodicalIF":2.7,"publicationDate":"2024-10-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142571409","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0

Slot coating of time-dependent structured materials: Effects of thixotropy on the flow dynamics and operating limits 随时间变化的结构材料槽涂层：触变性对流动动力学和工作极限的影响

IF 2.7 2区工程技术

Journal of Non-Newtonian Fluid Mechanics Pub Date : 2024-10-22 DOI: 10.1016/j.jnnfm.2024.105344

C. Sanchez-Perez, D. Maza, I.R. Siqueira, P.R. de Souza Mendes, M.S. Carvalho

{"title":"Slot coating of time-dependent structured materials: Effects of thixotropy on the flow dynamics and operating limits","authors":"C. Sanchez-Perez, D. Maza, I.R. Siqueira, P.R. de Souza Mendes, M.S. Carvalho","doi":"10.1016/j.jnnfm.2024.105344","DOIUrl":"10.1016/j.jnnfm.2024.105344","url":null,"abstract":"<div><div>Some of the most common liquid-like formulations rooted in the coating industry consist of rheologically complex structured materials such as inks, paints, and slurries. Yet, the impact of time-dependent structuring effects due to thixotropy on thin film coating applications remains elusive and still unclear. Here, we present a computational study of the effects of thixotropy on slot coating of time-dependent structured materials. By coupling a recent fluidity-based constitutive model for thixotropic materials with a well-established finite element/elliptic mesh generation method for free surface flows, we assess the role of thixotropy by comparing the predictions from the thixotropic model with those from a simple Generalized Newtonian Fluid model that uses the same flow curve for the material steady-state equilibrium viscosity. We find that thixotropy can indeed have a major impact on slot coating applications, potentially bringing strong implications not only to the flow dynamics but also to the operating limits of the process. In conclusion, the results and discussions we present in this study underscore the importance of accounting for thixotropy to genuinely model and fundamentally understand the behavior of time-dependent structured materials with complex rheology in processing flows like those ubiquitous across the broad fields of coating science and engineering.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"334 ","pages":"Article 105344"},"PeriodicalIF":2.7,"publicationDate":"2024-10-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142551993","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

引用次数: 0