Pub Date : 2025-12-30DOI: 10.1016/j.euromechflu.2025.204441
M.S. Faltas , E.A. Ashmawy , Samar A. Mahrous , M. Magdy El Sayed , Kareem E. Ragab
This study investigates the quasi-steady axisymmetric thermophoretic motion of a spherical particle partially submerged at the flat interface of a semi-infinite Brinkman medium. The analysis is conducted under the assumptions of small Reynolds and Péclet numbers, while the capillary number is considered sufficiently small to preserve the flatness of the interface. The specific case of a contact angle with the flat surface is examined. To avoid singularities at the contact line, the Knudsen number is assumed to lie within the slip-flow regime. Analytical expressions are derived for the thermophoretic velocity and force acting on the half-submerged particle. Graphical results illustrate the influence of parameters such as Fourier thermal conductivity ratio, Knudsen number, medium permeability, frictional slip, and thermal stress slip. Furthermore, the limiting behavior corresponding to thermophoresis in a classical viscous fluid is discussed. Since the present solution is exact, the case of a contact angle also serves as a benchmark for validating numerical solutions at other contact angles. The findings are relevant to applications involving particle manipulation at fluid–porous interfaces, such as targeted drug delivery across biological membranes, pollutant transport at soil–air boundaries, and the design of microfluidic systems for controlled colloidal assembly.
{"title":"Thermophoresis of a spherical particle straddling a flat interface in a Brinkman medium at a 90° contact angle","authors":"M.S. Faltas , E.A. Ashmawy , Samar A. Mahrous , M. Magdy El Sayed , Kareem E. Ragab","doi":"10.1016/j.euromechflu.2025.204441","DOIUrl":"10.1016/j.euromechflu.2025.204441","url":null,"abstract":"<div><div>This study investigates the quasi-steady axisymmetric thermophoretic motion of a spherical particle partially submerged at the flat interface of a semi-infinite Brinkman medium. The analysis is conducted under the assumptions of small Reynolds and Péclet numbers, while the capillary number is considered sufficiently small to preserve the flatness of the interface. The specific case of a <span><math><mrow><mn>90</mn><mo>°</mo></mrow></math></span> contact angle with the flat surface is examined. To avoid singularities at the contact line, the Knudsen number is assumed to lie within the slip-flow regime. Analytical expressions are derived for the thermophoretic velocity and force acting on the half-submerged particle. Graphical results illustrate the influence of parameters such as Fourier thermal conductivity ratio, Knudsen number, medium permeability, frictional slip, and thermal stress slip. Furthermore, the limiting behavior corresponding to thermophoresis in a classical viscous fluid is discussed. Since the present solution is exact, the case of a <span><math><mrow><mn>90</mn><mo>°</mo></mrow></math></span> contact angle also serves as a benchmark for validating numerical solutions at other contact angles. The findings are relevant to applications involving particle manipulation at fluid–porous interfaces, such as targeted drug delivery across biological membranes, pollutant transport at soil–air boundaries, and the design of microfluidic systems for controlled colloidal assembly.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204441"},"PeriodicalIF":2.5,"publicationDate":"2025-12-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880201","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-29DOI: 10.1016/j.euromechflu.2025.204430
N. Nouaime , B. Després , M.A. Puscas , C. Fiorini
This paper uses the intrusive polynomial chaos method (IPCM) to analyze sensitivity in heat transfer problems governed by the Navier–Stokes equations with heat transfer. The intrusive polynomial chaos method incorporates uncertain variables as combinations of orthogonal polynomials, known as polynomial chaos expansions (PCEs), directly into the governing equations. This transformation turns the original deterministic PDEs into coupled deterministic equations for the PCE coefficients. We apply first-order IPCM and propose a decoupling approach for state and sensitivity systems. We discretize the state equations and their sensitivity using the Finite Element-Volume (FEV) method. We establish a stability estimate for the continuous and discrete state and sensibility equations.
{"title":"Sensitivity analysis of Navier–Stokes equations with heat transfer using the first-order polynomial chaos method and FEV discretization","authors":"N. Nouaime , B. Després , M.A. Puscas , C. Fiorini","doi":"10.1016/j.euromechflu.2025.204430","DOIUrl":"10.1016/j.euromechflu.2025.204430","url":null,"abstract":"<div><div>This paper uses the intrusive polynomial chaos method (IPCM) to analyze sensitivity in heat transfer problems governed by the Navier–Stokes equations with heat transfer. The intrusive polynomial chaos method incorporates uncertain variables as combinations of orthogonal polynomials, known as polynomial chaos expansions (PCEs), directly into the governing equations. This transformation turns the original deterministic PDEs into coupled deterministic equations for the PCE coefficients. We apply first-order IPCM and propose a decoupling approach for state and sensitivity systems. We discretize the state equations and their sensitivity using the Finite Element-Volume (FEV) method. We establish a stability estimate for the continuous and discrete state and sensibility equations.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"118 ","pages":"Article 204430"},"PeriodicalIF":2.5,"publicationDate":"2025-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145923772","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-27DOI: 10.1016/j.euromechflu.2025.204451
Sanjay Kumar Pandey, Kriti Yadav
This study presents an analytical model that extends previous formulations by incorporating axial dependence in all velocity components, thus offering a more realistic representation of dust devil dynamics. We address the limitations of the fundamental model of Vyas and Majdalani (2006), such as unbounded velocity and infinite vortex width, by modifying the stream function, which captures radial and axial variation of velocity. The velocity components are physically bounded and consistent with observed vortex structures. Additionally, the radial pressure distribution is derived, and the flow field is analysed using normalized radial, axial, and azimuthal velocities. The azimuthal velocity is computed from the radial and axial flows using the angular momentum equation. It increases from the centre, peaks near the periphery, and gradually decays to zero toward the axis. This model represents the compact features of dust devils by producing a smooth, confined structure with consistent vertical dependence, unlike alternative models that introduce abrupt transitions or unbounded behaviour. The profile maintains central stability in accordance with the Rayleigh criterion, while peripheral centrifugal instability reflects the transient and dissipative nature of dust devils. The model captures dust lifting and charge separation in dust devils, linking strong updrafts and confined flow to enhanced electrification and particle transport.
{"title":"An analytical model for whirlwinds of finite width rising with bounded velocity and decaying outflow: Application to dust devils","authors":"Sanjay Kumar Pandey, Kriti Yadav","doi":"10.1016/j.euromechflu.2025.204451","DOIUrl":"10.1016/j.euromechflu.2025.204451","url":null,"abstract":"<div><div>This study presents an analytical model that extends previous formulations by incorporating axial dependence in all velocity components, thus offering a more realistic representation of dust devil dynamics. We address the limitations of the fundamental model of Vyas and Majdalani (2006), such as unbounded velocity and infinite vortex width, by modifying the stream function, which captures radial and axial variation of velocity. The velocity components are physically bounded and consistent with observed vortex structures. Additionally, the radial pressure distribution is derived, and the flow field is analysed using normalized radial, axial, and azimuthal velocities. The azimuthal velocity is computed from the radial and axial flows using the angular momentum equation. It increases from the centre, peaks near the periphery, and gradually decays to zero toward the axis. This model represents the compact features of dust devils by producing a smooth, confined structure with consistent vertical dependence, unlike alternative models that introduce abrupt transitions or unbounded behaviour. The profile maintains central stability in accordance with the Rayleigh criterion, while peripheral centrifugal instability reflects the transient and dissipative nature of dust devils. The model captures dust lifting and charge separation in dust devils, linking strong updrafts and confined flow to enhanced electrification and particle transport.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204451"},"PeriodicalIF":2.5,"publicationDate":"2025-12-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880202","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-26DOI: 10.1016/j.euromechflu.2025.204450
Manuel Rubio , Francisco Castro-Ruiz , José Sierra-Pallares , César Barrios-Collado , Joaquín Anatol
This study examines asymmetric pumping in Liebau pumps, a type of valveless pump that generates unidirectional flow through the periodic compression of a flexible (compliant) tube. This pumping has applications in biomedical devices, microfluidics, and organ support. We investigate how the properties of the compliant tube and the operating conditions affect the pump performance, using dimensionless parameters. Experiments were performed with different configurations, varying the tube material (latex and rubber) and the fluid (water and water–glycerine mixture). The results indicate that the flow rate and resonant period depend on the stiffness of the tube and the viscous effects. It was observed that for small values of the Womersley number (), viscous effects significantly reduce the flow rate. In contrast, for large Womersley values, the semiempirical models previously proposed adequately predict the experimental behaviour. Effects such as the compliant tube depression, which were not accounted for in previous models, were also found to influence performance. This work extends the analysis of this type of pump to unexplored conditions, with the aim of expanding knowledge and enabling the use of Liebau pumps in real-world applications.
{"title":"Analysis of compliant tube properties and operating conditions in a Liebau pump","authors":"Manuel Rubio , Francisco Castro-Ruiz , José Sierra-Pallares , César Barrios-Collado , Joaquín Anatol","doi":"10.1016/j.euromechflu.2025.204450","DOIUrl":"10.1016/j.euromechflu.2025.204450","url":null,"abstract":"<div><div>This study examines asymmetric pumping in Liebau pumps, a type of valveless pump that generates unidirectional flow through the periodic compression of a flexible (compliant) tube. This pumping has applications in biomedical devices, microfluidics, and organ support. We investigate how the properties of the compliant tube and the operating conditions affect the pump performance, using dimensionless parameters. Experiments were performed with different configurations, varying the tube material (latex and rubber) and the fluid (water and water–glycerine mixture). The results indicate that the flow rate and resonant period depend on the stiffness of the tube and the viscous effects. It was observed that for small values of the Womersley number (<span><math><msup><mrow><mi>W</mi><mi>o</mi></mrow><mrow><mn>2</mn></mrow></msup></math></span>), viscous effects significantly reduce the flow rate. In contrast, for large Womersley values, the semiempirical models previously proposed adequately predict the experimental behaviour. Effects such as the compliant tube depression, which were not accounted for in previous models, were also found to influence performance. This work extends the analysis of this type of pump to unexplored conditions, with the aim of expanding knowledge and enabling the use of Liebau pumps in real-world applications.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204450"},"PeriodicalIF":2.5,"publicationDate":"2025-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145837005","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-23DOI: 10.1016/j.euromechflu.2025.204449
Aakanksha Singh, Punit Kumar
Electron-acoustic solitary waves (EASWs) in quantum plasma comprising stationary ions, cold electrons, hot electrons, and kappa-distributed electrons have been investigated. The generalized Kappa-Fermi distribution has been modified to include electrostatic energy contribution and the number density of Kappa electrons has been obtained using this modified distribution. Utilizing the quantum hydrodynamic (QHD) model, a dispersion relation has been derived for linear EAWs. Employing the standard reductive perturbation technique, a Korteweg–de Vries (KdV) equation governing the dynamics of EAWs have been derived. The quantum mechanical effects of different parameters like kappa index, Mach number and equilibrium kappa electron density have been examined on the profiles of EASWs. It is found that the presence of kappa electrons in quantum plasma leads to new results, including steeper dispersion curves, sharper and more localized solitary waves with lower value of kappa index and stronger plasma interactions with increased kappa electrons density in dense astrophysical environment.
{"title":"Electron acoustic solitary wave in quantum plasmas with Kappa electrons","authors":"Aakanksha Singh, Punit Kumar","doi":"10.1016/j.euromechflu.2025.204449","DOIUrl":"10.1016/j.euromechflu.2025.204449","url":null,"abstract":"<div><div>Electron-acoustic solitary waves (EASWs) in quantum plasma comprising stationary ions, cold electrons, hot electrons, and kappa-distributed electrons have been investigated. The generalized Kappa-Fermi distribution has been modified to include electrostatic energy contribution and the number density of Kappa electrons has been obtained using this modified distribution. Utilizing the quantum hydrodynamic (QHD) model, a dispersion relation has been derived for linear EAWs. Employing the standard reductive perturbation technique, a Korteweg–de Vries (KdV) equation governing the dynamics of EAWs have been derived. The quantum mechanical effects of different parameters like kappa index, Mach number and equilibrium kappa electron density have been examined on the profiles of EASWs. It is found that the presence of kappa electrons in quantum plasma leads to new results, including steeper dispersion curves, sharper and more localized solitary waves with lower value of kappa index and stronger plasma interactions with increased kappa electrons density in dense astrophysical environment.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204449"},"PeriodicalIF":2.5,"publicationDate":"2025-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145837006","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-23DOI: 10.1016/j.euromechflu.2025.204448
Bruno Thierry Nyatchouba Nsangue , Tang Hao , Mezoue Adiang Cyrille , Liuxiong Xu , Fuxiang Hu , Ruben Mouangue
The unsteady turbulent flow, the trawl codend, and the Bycatch reduction device (BRD) are characterized by a complex interaction. A thorough understanding of this interaction is crucial for minimizing bycatch, enhancing the escape probability of non-target species, and improving trawl selectivity. This study analyzes unsteady turbulent flow fields inside and around the bottom trawl codends equipped with BRDs, systematically exploring the reconstruction capabilities of three deep learning algorithms, involving Long Short-Term Memory (LSTM), Echo State Network (ESN), and Convolutional Neural Network-Long Short-Term Memory (CNN-LSTM). Extensive assessment encompasses the predictive reliability of the instantaneous flow fields, velocity ratio profiles, turbulent intensity, turbulent kinetic energy, Reynold stress, time series streamwise flow velocities, and training efficiency. Results indicate a full development of unsteady turbulent flow inside and around the codend without BRD, characterized by instantaneous shear layer instabilities and vortex shedding. In contrast, the codend-BRD system and water flow interaction generates two distinct regions with very low flow velocity fields behind it, characterized by vortex-shedding structures on the unsteady turbulent wake. The flow velocity deficit was greater inside the codend compared to that observed inside the combined codend-BRD system due to free water flow passage through the combined codend-BRD system due to the presence of BRD windows. A higher turbulent kinetic energy, greater momentum flux, and stronger turbulence intensities are observed inside and around the codend without BRD compared to the codend-BRD system. Additionally, the results indicate that ESN and CNN-LSTM exhibit a significant advantage in reconstructing unsteady turbulent flow parameters. The prediction of unsteady turbulent flow parameters indicates that all the three models are substantially consistent with the experimental data. However, LSTM is a judicious choice when solely the time series variables require prediction. These insights significantly advance the development of smart trawl nets, offering a pathway to enhance gear selectivity through data-driven design.
{"title":"Deep learning prediction of unsteady turbulent transformation pattern inside and around the interactions between a sorting grid and a trawl system","authors":"Bruno Thierry Nyatchouba Nsangue , Tang Hao , Mezoue Adiang Cyrille , Liuxiong Xu , Fuxiang Hu , Ruben Mouangue","doi":"10.1016/j.euromechflu.2025.204448","DOIUrl":"10.1016/j.euromechflu.2025.204448","url":null,"abstract":"<div><div>The unsteady turbulent flow, the trawl codend, and the Bycatch reduction device (BRD) are characterized by a complex interaction. A thorough understanding of this interaction is crucial for minimizing bycatch, enhancing the escape probability of non-target species, and improving trawl selectivity. This study analyzes unsteady turbulent flow fields inside and around the bottom trawl codends equipped with BRDs, systematically exploring the reconstruction capabilities of three deep learning algorithms, involving Long Short-Term Memory (LSTM), Echo State Network (ESN), and Convolutional Neural Network-Long Short-Term Memory (CNN-LSTM). Extensive assessment encompasses the predictive reliability of the instantaneous flow fields, velocity ratio profiles, turbulent intensity, turbulent kinetic energy, Reynold stress, time series streamwise flow velocities, and training efficiency. Results indicate a full development of unsteady turbulent flow inside and around the codend without BRD, characterized by instantaneous shear layer instabilities and vortex shedding. In contrast, the codend-BRD system and water flow interaction generates two distinct regions with very low flow velocity fields behind it, characterized by vortex-shedding structures on the unsteady turbulent wake. The flow velocity deficit was greater inside the codend compared to that observed inside the combined codend-BRD system due to free water flow passage through the combined codend-BRD system due to the presence of BRD windows. A higher turbulent kinetic energy, greater momentum flux, and stronger turbulence intensities are observed inside and around the codend without BRD compared to the codend-BRD system. Additionally, the results indicate that ESN and CNN-LSTM exhibit a significant advantage in reconstructing unsteady turbulent flow parameters. The prediction of unsteady turbulent flow parameters indicates that all the three models are substantially consistent with the experimental data. However, LSTM is a judicious choice when solely the time series variables require prediction. These insights significantly advance the development of smart trawl nets, offering a pathway to enhance gear selectivity through data-driven design.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204448"},"PeriodicalIF":2.5,"publicationDate":"2025-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880203","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-22DOI: 10.1016/j.euromechflu.2025.204446
Anupam M. Hiremath , Harunori N. Yoshikawa , Innocent Mutabazi
Liquid Metal Battery (LMB) is a recent technology proposed to store the energy produced by intermittent sources such as solar cells, wind turbines or tidal power plants. The liquid metal batteries operate with large current densities but with very low open circuit voltages (OCV). A LMB consists of superimposed three layers of immiscible liquids in a stable stratification: the top layer made of light liquid metal and the bottom one made of dense liquid metal alloy sandwich the electrolyte layer (often a molten salt). Each layer has its own specific fluid properties and due to large difference in thermal and electrical conductivities between the liquid metals and the electrolyte, the application of a large current density can generate thermoconvective instabilities. The present works presents, in the framework of the linear stability analysis, the effects of the physical properties and thicknesses of liquid layers on the critical parameters on thermoconvective flows induced by Joule heating in a LMB of rectangular cross-section. The thermal coupling of the LMB with its environment described by the mixed thermal boundary conditions through the Biot number modifies the temperature profiles of the conduction state and reduces the threshold of all thermoconvective flows.
{"title":"Threshold of thermoconvective flows in a liquid metal battery with thermal coupling with its environment","authors":"Anupam M. Hiremath , Harunori N. Yoshikawa , Innocent Mutabazi","doi":"10.1016/j.euromechflu.2025.204446","DOIUrl":"10.1016/j.euromechflu.2025.204446","url":null,"abstract":"<div><div>Liquid Metal Battery (LMB) is a recent technology proposed to store the energy produced by intermittent sources such as solar cells, wind turbines or tidal power plants. The liquid metal batteries operate with large current densities but with very low open circuit voltages (OCV). A LMB consists of superimposed three layers of immiscible liquids in a stable stratification: the top layer made of light liquid metal and the bottom one made of dense liquid metal alloy sandwich the electrolyte layer (often a molten salt). Each layer has its own specific fluid properties and due to large difference in thermal and electrical conductivities between the liquid metals and the electrolyte, the application of a large current density can generate thermoconvective instabilities. The present works presents, in the framework of the linear stability analysis, the effects of the physical properties and thicknesses of liquid layers on the critical parameters on thermoconvective flows induced by Joule heating in a LMB of rectangular cross-section. The thermal coupling of the LMB with its environment described by the mixed thermal boundary conditions through the Biot number modifies the temperature profiles of the conduction state and reduces the threshold of all thermoconvective flows.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204446"},"PeriodicalIF":2.5,"publicationDate":"2025-12-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145921070","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-18DOI: 10.1016/j.euromechflu.2025.204447
Yupeng Liao , Changjun Li , Wenlong Jia , Juncheng Mu , Jie He , Fan Yang , Qiaojing Huang
This paper proposes a new segregated algorithm named SIMPLE-Revised with Residual Coupling (SIMPLER*) for simulating the two-phase flow of gas-condensate in pipelines. The full coupling of velocity and pressure is ensured by introducing the momentum equation residual and the double inner iteration, which improves the computational efficiency while solving the issues of velocity-pressure decoupling and pressure overcorrection in the traditional Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. The adaptive time-stepping method is employed to reduce computational costs in transient simulations. The performance of SIMPLE, SIMPLE-Revised-Revised (SIMPLERR), and SIMPLER* algorithms is comprehensively compared. Validation across four classical benchmark tests confirms the superior robustness and numerical stability of SIMPLER*, which exhibits lower numerical diffusion and remains oscillation-free under strong pressure discontinuities. Subsequent steady-state and transient engineering cases demonstrate that the computational efficiency of the SIMPLER* algorithm is significantly better than that of the SIMPLE and SIMPLERR algorithms, both in coarse and fine grids. The proposed "adaptive time-stepping + SIMPLER* algorithm" combination method greatly reduces the computational overhead while ensuring the accuracy of the solution, providing an efficient and reliable solution for the transient simulation of complex two-phase flow with phase change.
{"title":"Insight into the numerical efficiency and stability of a residual-coupled segregated solver for gas-condensate two-phase flow in pipelines: An algorithmic study","authors":"Yupeng Liao , Changjun Li , Wenlong Jia , Juncheng Mu , Jie He , Fan Yang , Qiaojing Huang","doi":"10.1016/j.euromechflu.2025.204447","DOIUrl":"10.1016/j.euromechflu.2025.204447","url":null,"abstract":"<div><div>This paper proposes a new segregated algorithm named SIMPLE-Revised with Residual Coupling (SIMPLER*) for simulating the two-phase flow of gas-condensate in pipelines. The full coupling of velocity and pressure is ensured by introducing the momentum equation residual and the double inner iteration, which improves the computational efficiency while solving the issues of velocity-pressure decoupling and pressure overcorrection in the traditional Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) algorithm. The adaptive time-stepping method is employed to reduce computational costs in transient simulations. The performance of SIMPLE, SIMPLE-Revised-Revised (SIMPLERR), and SIMPLER* algorithms is comprehensively compared. Validation across four classical benchmark tests confirms the superior robustness and numerical stability of SIMPLER*, which exhibits lower numerical diffusion and remains oscillation-free under strong pressure discontinuities. Subsequent steady-state and transient engineering cases demonstrate that the computational efficiency of the SIMPLER* algorithm is significantly better than that of the SIMPLE and SIMPLERR algorithms, both in coarse and fine grids. The proposed \"adaptive time-stepping + SIMPLER* algorithm\" combination method greatly reduces the computational overhead while ensuring the accuracy of the solution, providing an efficient and reliable solution for the transient simulation of complex two-phase flow with phase change.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204447"},"PeriodicalIF":2.5,"publicationDate":"2025-12-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145837007","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-16DOI: 10.1016/j.euromechflu.2025.204440
Siyi Li , Zihao Zhu , Lei Xie , Yaguang Xie , Ruonan Wang , Qiang Du , Junqiang Zhu
Instabilities in rotor–stator cavities significantly influence flow dynamics and heat transfer processes within aerospace propulsion systems. Among these instabilities, the circular waves manifests within the stator boundary layer and exhibits transient behavior highly sensitive to disturbances in the basic state. To elucidate the underlying mechanisms driving this transient phenomenon, through direct numerical simulation (DNS), we systematically imposed impulsive changes and harmonic modulations on the rotation speed of the rotor, thereby inducing various types of disturbances. Our findings reveal that the emergence of circular waves is triggered by these disturbances, with the waves’ characteristics displaying marked sensitivity to the nature of the disturbances. Specifically, increasing the disturbance frequency leads to an upward migration of the circular waves’ radial position. The energy and radial extent of the circular waves exhibit a trend of first increasing and then decreasing as the disturbance frequency increases. Moreover, as the disturbance amplitude increases, the radial extent occupied by the circular waves expands, while the midpoint of their radial position remains unaltered. We further identified that when a hub rotating with the rotor, circular waves can become self-sustaining under certain conditions. Specifically, when the gap ratio (, where is the radius of the hub, is the radius of the shroud, is the half of the gap between the rotor and stator) and rotational Reynolds number (, where refers to the rotational speed of the rotor, and refers to the kinematic viscosity) are sufficiently large, disturbances on the stator side can migrate through the hub, amplify, and form disturbances on the rotor side, subsequently re-exciting circular waves on the stator. Through linear stability analysis, we determined the boundary in the parameter domain that delineates conditions for self-sustaining circular waves. This study provides a comprehensive investigation into the behavior of circular waves, shedding new light on their complex dynamics within rotor–stator cavities.
{"title":"Direct numerical simulation and linear stability analysis of circular waves in the stator boundary layer of rotor–stator cavity","authors":"Siyi Li , Zihao Zhu , Lei Xie , Yaguang Xie , Ruonan Wang , Qiang Du , Junqiang Zhu","doi":"10.1016/j.euromechflu.2025.204440","DOIUrl":"10.1016/j.euromechflu.2025.204440","url":null,"abstract":"<div><div>Instabilities in rotor–stator cavities significantly influence flow dynamics and heat transfer processes within aerospace propulsion systems. Among these instabilities, the circular waves manifests within the stator boundary layer and exhibits transient behavior highly sensitive to disturbances in the basic state. To elucidate the underlying mechanisms driving this transient phenomenon, through direct numerical simulation (DNS), we systematically imposed impulsive changes and harmonic modulations on the rotation speed of the rotor, thereby inducing various types of disturbances. Our findings reveal that the emergence of circular waves is triggered by these disturbances, with the waves’ characteristics displaying marked sensitivity to the nature of the disturbances. Specifically, increasing the disturbance frequency leads to an upward migration of the circular waves’ radial position. The energy and radial extent of the circular waves exhibit a trend of first increasing and then decreasing as the disturbance frequency increases. Moreover, as the disturbance amplitude increases, the radial extent occupied by the circular waves expands, while the midpoint of their radial position remains unaltered. We further identified that when a hub rotating with the rotor, circular waves can become self-sustaining under certain conditions. Specifically, when the gap ratio (<span><math><mrow><mi>γ</mi><mo>=</mo><mrow><mo>(</mo><mi>b</mi><mo>−</mo><mi>a</mi><mo>)</mo></mrow><mo>/</mo><mi>H</mi></mrow></math></span>, where <span><math><mi>a</mi></math></span> is the radius of the hub, <span><math><mi>b</mi></math></span> is the radius of the shroud, <span><math><mi>H</mi></math></span> is the half of the gap between the rotor and stator) and rotational Reynolds number (<span><math><mrow><mi>R</mi><mi>e</mi><mo>=</mo><mi>Ω</mi><msup><mrow><mi>b</mi></mrow><mrow><mn>2</mn></mrow></msup><mo>/</mo><mi>ν</mi></mrow></math></span>, where <span><math><mi>Ω</mi></math></span> refers to the rotational speed of the rotor, and <span><math><mi>ν</mi></math></span> refers to the kinematic viscosity) are sufficiently large, disturbances on the stator side can migrate through the hub, amplify, and form disturbances on the rotor side, subsequently re-exciting circular waves on the stator. Through linear stability analysis, we determined the boundary in the <span><math><mrow><mo>(</mo><mi>R</mi><mi>e</mi><mo>,</mo><mi>γ</mi><mo>)</mo></mrow></math></span> parameter domain that delineates conditions for self-sustaining circular waves. This study provides a comprehensive investigation into the behavior of circular waves, shedding new light on their complex dynamics within rotor–stator cavities.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204440"},"PeriodicalIF":2.5,"publicationDate":"2025-12-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145837011","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2025-12-15DOI: 10.1016/j.euromechflu.2025.204445
Valerio Lupi , Filippo Caruso Lombardi , Martino Andrea Scarpolini , Roberto Verzicco , Francesco Viola
Although blood exhibits non-Newtonian properties, it is often modelled as a Newtonian fluid in numerical studies. However, the impact of this assumption on the accuracy of the predictions of computational models is not thoroughly assessed. For this reason, in this work, we compare the effect of Newtonian and shear-thinning blood rheological models on computational haemodynamics in the left heart. The analysis is based on the fluid–structure-electrophysiology interaction (FSEI) within a patient-specific anatomy of an individual with low ejection fraction to better highlight viscosity changes at low strain rates. Owing to the pulsatile nature of the flow, considerable spatio-temporal variations of the kinematic viscosity are observed for the non-Newtonian case. Integral quantities as well as blood pressure within the cardiac chambers are found to be weakly affected by the rheological model. Substantial differences are reported, instead, between the Newtonian and non-Newtonian cases, for the wall shear stress and local haemolysis index since they are sensitive to the local kinematic viscosity, which decreases within the shear layers of the systolic and diastolic jets, as well as in the near-wall regions.
{"title":"Impact of blood rheology on left heart haemodynamics: Newtonian vs. non-Newtonian modelling","authors":"Valerio Lupi , Filippo Caruso Lombardi , Martino Andrea Scarpolini , Roberto Verzicco , Francesco Viola","doi":"10.1016/j.euromechflu.2025.204445","DOIUrl":"10.1016/j.euromechflu.2025.204445","url":null,"abstract":"<div><div>Although blood exhibits non-Newtonian properties, it is often modelled as a Newtonian fluid in numerical studies. However, the impact of this assumption on the accuracy of the predictions of computational models is not thoroughly assessed. For this reason, in this work, we compare the effect of Newtonian and shear-thinning blood rheological models on computational haemodynamics in the left heart. The analysis is based on the fluid–structure-electrophysiology interaction (FSEI) within a patient-specific anatomy of an individual with low ejection fraction to better highlight viscosity changes at low strain rates. Owing to the pulsatile nature of the flow, considerable spatio-temporal variations of the kinematic viscosity are observed for the non-Newtonian case. Integral quantities as well as blood pressure within the cardiac chambers are found to be weakly affected by the rheological model. Substantial differences are reported, instead, between the Newtonian and non-Newtonian cases, for the wall shear stress and local haemolysis index since they are sensitive to the local kinematic viscosity, which decreases within the shear layers of the systolic and diastolic jets, as well as in the near-wall regions.</div></div>","PeriodicalId":11985,"journal":{"name":"European Journal of Mechanics B-fluids","volume":"117 ","pages":"Article 204445"},"PeriodicalIF":2.5,"publicationDate":"2025-12-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145837008","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
扫码关注我们
求助内容：
应助结果提醒方式：
京公网安备 11010802042870号