Journal of Non-Newtonian Fluid Mechanics最新文献

{"title":"CFD-DEM simulation of dynamic filtration in heterogeneous porous media with viscoplastic and shear-thinning fluids","authors":"Ayrton Cavallini Zotelle ,&nbsp;Vinicius Gustavo Poletto ,&nbsp;Felipe Barboza Pereira ,&nbsp;Fernando Cesar De Lai ,&nbsp;Renato do Nascimento Siqueira ,&nbsp;Silvio Luiz de Mello Junqueira","doi":"10.1016/j.jnnfm.2025.105558","DOIUrl":"10.1016/j.jnnfm.2025.105558","url":null,"abstract":"<div><div>This study employs a coupled Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) to investigate the dynamic filtration behavior of non-Newtonian fluids in highly permeable porous media. The research explores the influence of rheological parameters, specifically the plasticity number (<span><math><mrow><mi>P</mi><mi>l</mi></mrow></math></span>) and power-law index (<span><math><mi>n</mi></math></span>), over filter cake shape, filtered mass, and flow rate through the porous medium. Simulations consider Newtonian, shear-thinning, and viscoplastic fluids flowing through a porous matrix composed of a staggered array of cylindrical obstacles. The study also explores the impact of mesh resolution and regularization parameters on simulation stability. Results indicate that higher plasticity enhances sealing performance by avoiding settling and building evenly thickened particulate layers. Meanwhile, shear-thinning behavior increases local viscosity in low-shear zones, reducing the local flow rate and, therefore, mass retention. The combination of low <span><math><mi>n</mi></math></span> and high <span><math><mrow><mi>P</mi><mi>l</mi></mrow></math></span> yields the most effective filtration, minimizing fluid loss and retained mass. Findings highlight the critical role of fluid rheology in optimizing dynamic filtration and suggest that tailoring <span><math><mrow><mi>P</mi><mi>l</mi></mrow></math></span> and <span><math><mi>n</mi></math></span> can significantly improve the process.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105558"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145925954","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}

{"title":"Flow-induced migration of aggregating suspensions in pipe flow","authors":"Soham Jariwala ,&nbsp;Norman J. Wagner,&nbsp;Antony N. Beris","doi":"10.1016/j.jnnfm.2025.105533","DOIUrl":"10.1016/j.jnnfm.2025.105533","url":null,"abstract":"<div><div>Aggregating suspensions can be found in many materials, such as food products, biological fluids, printer inks, paints, and slurries. These suspensions display unique viscoelastic and thixotropic behavior due to the way the agglomerates formed by interparticle attraction undergo elastic deformation, aggregation, and breakage. In this work we show how it is possible to describe inhomogeneities induced by stress-induced migration by using a population balance-based constitutive model [Mwasame et al. AIChE J. 63 (2017) 517-531]. An important advantage of this model over phenomenological constitutive models, such as structure kinetics models, is that it can describe both the kinetics of aggregation and breakage as well as the migration fluxes in terms of the local microscopic structural descriptors, such as the volume fraction of primary particles, the fractal dimension, and the moments of the agglomerate size distribution.</div><div>In the present work, we adapt this mesoscale structural description to resolve the flow-induced migration behavior exhibited in pipe flows by coupling with a modified diffusive flux model proposed by Phillips et al. [Phys. Fluids 4 (1992) 30-40]. Our study focuses on fully-developed, pressure-driven (Poiseuille) flows in a tube, both steady and transient. We use numerical simulations based on Chebyshev orthogonal polynomial approximations of the variables along the radial directions within an efficient Galerkin weighted residuals methodology. The development of concentration inhomogeneities is investigated along with their effects on the stress that arises from thixotropy and viscoelasticity in both steady state and transient flows. Additionally, we use wall slip to appropriately model the particle-free layer than forms near the tube wall. Of potential significance to applications is the observation that the superposition of an oscillatory pressure gradient to a steady one can lead to a reduction in the power dissipation for a given average flow rate.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105533"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145791261","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}

{"title":"The effects of viscoelasticity and shear thinning on the inertial instability and mixing performance in a T-channel geometry","authors":"Rebecca J. Hill ,&nbsp;Mahdi Davoodi ,&nbsp;Simon J. Haward ,&nbsp;Cláudio P. Fonte ,&nbsp;Robert J. Poole","doi":"10.1016/j.jnnfm.2026.105564","DOIUrl":"10.1016/j.jnnfm.2026.105564","url":null,"abstract":"<div><div>The T-channel, in which two opposing square inlet streams join and turn through <span><math><mrow><mtext>90°</mtext></mrow></math></span> into a rectangular outlet of equal cross-sectional area, exhibits a steady symmetry-breaking bifurcation above a critical Reynolds number. For Newtonian fluids it is well-established that this transition produces engulfment flow and consequent enhanced mixing. However, the transition behaviour and mixing of viscoelastic and shear thinning fluids in the T-channel is less well studied or understood. In this work, finite-volume flow simulations are used to investigate how non-Newtonian rheology influences both the onset of instability and the resulting mixing performance. For constant-viscosity viscoelastic models, the Oldroyd-B fluid destabilises the flow, leading to transition at lower Reynolds numbers than in the Newtonian case, while the Oldroyd-A fluid stabilises the flow, delaying the onset of the engulfment regime. This contrast highlights the influence of normal stress differences on critical conditions. For shear-thinning Giesekus and Carreau models, the critical Reynolds number depends strongly on how the Reynolds number is defined; using the zero-shear viscosity both predict destabilisation relative to the Newtonian baseline. Analysis of the outlet channel shows that the instability mechanism is governed by the strength of secondary Dean-type vortices generated by the <span><math><mrow><mtext>90°</mtext></mrow></math></span> turn, which are amplified or suppressed depending on the balance between inertia, curvature and rheology. Finally, the quantification of the mixing index reveals that, despite shifting the onset of instability, all non-Newtonian models studied reduce mixing efficiency of the mixer relative to the Newtonian case.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105564"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395605","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}

{"title":"Experimental study of shear thinning effects on solitary wave propagation: A Newtonian fluid comparison","authors":"C. Calvo ,&nbsp;A. Tamburrino ,&nbsp;C. Falcón","doi":"10.1016/j.jnnfm.2026.105560","DOIUrl":"10.1016/j.jnnfm.2026.105560","url":null,"abstract":"<div><div>We present an experimental study of viscous effects in transient non-linear long wave propagation in two Newtonian fluids and one shear thinning fluid in the laminar flow regime. Using optical measuring techniques (Fourier Transform Profilometry), we show that the wave phase speed decreases in both glycerin and carboxymethylcellulose (CMC) solutions with respect to that in water. A decrease in wave phase speed is observed, and a dispersion relation is obtained for surface waves through dimensional analysis from five dimensionless groups: the dimensionless wave celerity, the shallowness parameter, dimensionless amplitude, Reynolds number and the flow index. To complete the picture on wave propagation, an empirical dependence between the wave attenuation and the last four dimensionless groups mentioned above is found for non-linear long surface waves in our working fluids. We conclude quantitatively about viscosity effects in non-linear long wave propagation.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105560"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146077775","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}

{"title":"Turbulent wall-bounded flow of weakly shear-thinning fluid: Phenomenological modeling of eddy viscosity and law of the wall","authors":"Ernest Simon,&nbsp;Alain Liné","doi":"10.1016/j.jnnfm.2026.105563","DOIUrl":"10.1016/j.jnnfm.2026.105563","url":null,"abstract":"<div><div>This work proposes a new analytical law of the wall for turbulent, wall-bounded flows of generalized Newtonian fluids. Direct numerical simulation data from the literature are processed to analyze and model the total (molecular + turbulent) viscosity profiles in pipe flows for power-law fluids with flow indices in the range <span><math><mrow><mn>0</mn><mo>.</mo><mn>4</mn><mo>&lt;</mo><mi>n</mi><mo>&lt;</mo><mn>0</mn><mo>.</mo><mn>8</mn></mrow></math></span> and generalized Reynolds numbers between <span><math><mrow><mn>11</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mn>3</mn></mrow></msup></mrow></math></span> and <span><math><mrow><mn>20</mn><mo>×</mo><mn>1</mn><msup><mrow><mn>0</mn></mrow><mrow><mn>3</mn></mrow></msup></mrow></math></span>. The total viscosity is expressed as a nonlinear function of the wall-normal coordinate <span><math><msup><mrow><mi>y</mi></mrow><mrow><mo>+</mo></mrow></msup></math></span>, revealing a deviation from the classical linear dependence valid for Newtonian fluids. A simplified empirical formulation of the turbulent viscosity is proposed, leading to a non-logarithmic mean velocity profile that captures the effects of shear thinning. The resulting analytical velocity profiles show good agreement with literature data and reproduce the observed upward shift of the mean velocity in the near-wall region. From this law, expressions for the Fanning friction factor and drag reduction coefficient are derived. The analysis captures the drag reduction increase with decreasing flow index and decreasing generalized Reynolds number, consistent with the low drag reduction (LDR) regime reported in previous numerical studies. In order to generalize our results, it will be necessary to explore a wider range of flow indices and Reynolds numbers.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105563"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395717","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}

{"title":"Special issue - ICTAM - non-Newtonian and complex fluids. Editorial","authors":"Anke Lindner ,&nbsp;Prabhu R Nott","doi":"10.1016/j.jnnfm.2026.105561","DOIUrl":"10.1016/j.jnnfm.2026.105561","url":null,"abstract":"","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105561"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147448993","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}

{"title":"Experimental investigation of movement and deposition of woody-debris suspensions in inclined channel tests","authors":"Chyan-Deng Jan,&nbsp;Le-Trang Nguyen","doi":"10.1016/j.jnnfm.2025.105547","DOIUrl":"10.1016/j.jnnfm.2025.105547","url":null,"abstract":"&lt;div&gt;&lt;div&gt;Debris flows, which mobilize large volumes of water, sediment, and woody debris, pose significant risks to human communities and infrastructure. In wildfire-affected forested areas, the accumulation of woody debris in drainage channels is exacerbated, thereby increasing the potential for more hazardous debris flows. To examine the influence of woody debris on debris flow dynamics, an inclined channel test, allowing the observation of woody debris flow in an inclined channel and its deposition in a horizontal tank, was conducted with highly concentrated woody-debris suspensions composed of clay, silt, woody debris, and water. This study explores how variations in fine-sediment fraction (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;f&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;), woody debris proportion (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;g&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;), and woody debris size (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;S&lt;/mi&gt;&lt;mi&gt;w&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;) impact flow behavior, including the entry speed (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;V&lt;/mi&gt;&lt;mn&gt;0&lt;/mn&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;) into the horizontal tank, and deposition characteristics such as runout distance (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;), deposit width (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;W&lt;/mi&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;), deposit thickness (&lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;), and final profiles on a &lt;span&gt;&lt;math&gt;&lt;mrow&gt;&lt;mn&gt;2&lt;/mn&gt;&lt;msup&gt;&lt;mrow&gt;&lt;mn&gt;0&lt;/mn&gt;&lt;/mrow&gt;&lt;mi&gt;o&lt;/mi&gt;&lt;/msup&gt;&lt;/mrow&gt;&lt;/math&gt;&lt;/span&gt; channel slope. To examine the influence of proportions and sizes of woody debris on the entry speed, an empirical equation is presented relating &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;V&lt;/mi&gt;&lt;mn&gt;0&lt;/mn&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; to &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;f&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;, &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;g&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt;, and &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;S&lt;/mi&gt;&lt;mi&gt;w&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; using multiple linear regression analysis. The results indicate that a higher &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;f&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; and &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mrow&gt;&lt;mi&gt;v&lt;/mi&gt;&lt;mi&gt;g&lt;/mi&gt;&lt;/mrow&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; yields smaller entry speeds, leading to shorter runout distances, thicker deposits, and wider deposit extents. The tests of larger &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;S&lt;/mi&gt;&lt;mi&gt;w&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; generate larger entry speeds, resulting in longer runout distances while producing thinner and narrower deposits. Empirical equations relating &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;V&lt;/mi&gt;&lt;mn&gt;0&lt;/mn&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; to &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;L&lt;/mi&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; and &lt;span&gt;&lt;math&gt;&lt;msub&gt;&lt;mi&gt;W&lt;/mi&gt;&lt;mi&gt;R&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt;&lt;/span&gt; are also provided to further demonstrate the influence of entry speeds on the deposit characteristics. Additionally, a strong correlation was found between inclined channel test parameters (e.g., entry speed, runout distance, and maximum deposit width) and rheological parame","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105547"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145791262","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}

{"title":"Non-monotonic die swell of liquid foams","authors":"Mehdi Habibi ,&nbsp;Etienne Jambon-Puillet ,&nbsp;Daniel Bonn","doi":"10.1016/j.jnnfm.2026.105570","DOIUrl":"10.1016/j.jnnfm.2026.105570","url":null,"abstract":"<div><div>Liquid jets extruded through narrow orifices can expand significantly upon exiting, a phenomenon known as die swell. While extensively studied for polymers, die swell in other complex fluids remains poorly understood. We present an experimental investigation of die swell in foams extruded from a capillary tube at the outlet of a syringe (whipped cream and shaving foam). The experiments show a very large die swell that is, in addition, non-monotonic: contrary to what happens for polymers, with increasing flow rate, the die swell of the foam first decreases, passes through a minimum and then increases again. We show that these two regimes can be explained by a balance of capillary and viscous forces at low pressure and the compressibility of the foam at high pressure, respectively. These results will facilitate the prediction of extruded filaments' diameter in the additive manufacturing of aerated liquids, with applications in food, medicine, or construction. We show that die swell scales as 1/Ca at low capillary numbers and as D/D₀ ∼ (1+ΔP/P₀)^1/3 at high pressures. Die swell ratios reach ∼2 for small capillaries (D₀ &lt; 1 mm) at low pressures, decrease to ∼1 at intermediate pressures (0.15–0.2 bar), then increase again at higher pressures<strong>.</strong></div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105570"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147395607","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}

{"title":"The Giesekus model revisited","authors":"Fatemeh Karami ,&nbsp;Pavlos S. Stephanou","doi":"10.1016/j.jnnfm.2025.105546","DOIUrl":"10.1016/j.jnnfm.2025.105546","url":null,"abstract":"<div><div>Since its introduction, the Giesekus model has received increased attention, particularly due to its ability to provide a non-vanishing second normal stress difference in simple shear. However, its derivation was based on a postulate regarding the linearity between the mobility tensor and the conformation tensor. In this work, we elaborate on the implications of this linearity and examine how its predictions are altered when the second-order and third-order corrections are considered. The predictions of the non-linear versions of the Giesekus model are found to be partially in better agreement with experimental and simulation data than those of the original Giesekus model in the case of uniaxial elongational flow.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105546"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145791260","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}

{"title":"A sharp computational method for simulating multiphase viscoelastic flows","authors":"Joseph V. Giliberto,&nbsp;Olivier Desjardins","doi":"10.1016/j.jnnfm.2025.105559","DOIUrl":"10.1016/j.jnnfm.2025.105559","url":null,"abstract":"<div><div>Viscoelastic constitutive equations often model the elastic stress field through the use of an elastic dumbbell model that utilizes a conformation tensor to represent the average polymer configuration in the flow field. In a liquid–gas flow environment, the conformation tensor is a discontinuous quantity that only exists in the liquid phase. This discontinuity often presents numerical challenges that can be tackled through the use of very fine meshes at the interface to ensure the stress profile is accurately captured. In contrast, this work presents a hybrid advection scheme for the discontinuous conformation tensor field that uses a semi-Lagrangian geometric flux-based scheme in the direct vicinity of the liquid–gas interface and a MUSCL scheme in the bulk of the liquid, away from the interface. This hybrid method is found to be exactly conservative and bounded, and prevents any leakage of data across the liquid–gas interface. Verification and validation of this approach is done using the case of a gas bubble rising in a viscoelastic liquid. Results of the convergence study show that the hybrid scheme is able to converge to experimental results with 32 cells across the initial diameter of the bubble, which is one-third the resolution used in other computational studies comparing against experiments. The hybrid advection scheme is then applied to the case of a viscoelastic droplet deforming in homogeneous isotropic turbulence to investigate the influence of elastic stresses on droplet morphology. Results indicate that increasing viscoelastic stresses within the droplet significantly alters its deformation dynamics. At the moderate elastic stress levels tested, the droplet forms elongated liquid filaments delaying break-up for a longer duration. As viscoelasticity is further increased, deformation is progressively suppressed, ultimately stabilizing the droplet’s shape and preventing fragmentation.</div></div>","PeriodicalId":54782,"journal":{"name":"Journal of Non-Newtonian Fluid Mechanics","volume":"348 ","pages":"Article 105559"},"PeriodicalIF":2.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145925955","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}