Daniel B. V. Santos, Gustavo P. Oliveira, N. Mangiavacchi, Prashant Valluri, G. Anjos
This work’s goal is to numerically investigate the interactions between two gas bubbles in a fluid flow in a circular cross-section channel, both in the presence and in the absence of gravitational forces, with several Reynolds and Weber numbers. The first bubble is placed at the center of the channel, while the second is near the wall. Their positions are set in such a way that a dynamic interaction is expected to occur due to their velocity differences. A finite element numerical tool is utilized to solve the incompressible Navier–Stokes equations and simulate two-phase flow using an unfitted mesh to represent the fluid interface, akin to the front-tracking method. The results show that the velocity gradient influences bubble shapes near the wall. Moreover, lower viscosity and surface tension force account for more significant interactions, both in the bubble shape and in the trajectory. In this scenario, it can be observed that one bubble is trapped in the other’s wake, with the proximity possibly allowing the onset of coalescence. The results obtained contribute to a deeper understanding of two-phase inner flows.
{"title":"Numerical Investigation of Gas Bubble Interaction in a Circular Cross-Section Channel in Shear Flow","authors":"Daniel B. V. Santos, Gustavo P. Oliveira, N. Mangiavacchi, Prashant Valluri, G. Anjos","doi":"10.3390/fluids9020032","DOIUrl":"https://doi.org/10.3390/fluids9020032","url":null,"abstract":"This work’s goal is to numerically investigate the interactions between two gas bubbles in a fluid flow in a circular cross-section channel, both in the presence and in the absence of gravitational forces, with several Reynolds and Weber numbers. The first bubble is placed at the center of the channel, while the second is near the wall. Their positions are set in such a way that a dynamic interaction is expected to occur due to their velocity differences. A finite element numerical tool is utilized to solve the incompressible Navier–Stokes equations and simulate two-phase flow using an unfitted mesh to represent the fluid interface, akin to the front-tracking method. The results show that the velocity gradient influences bubble shapes near the wall. Moreover, lower viscosity and surface tension force account for more significant interactions, both in the bubble shape and in the trajectory. In this scenario, it can be observed that one bubble is trapped in the other’s wake, with the proximity possibly allowing the onset of coalescence. The results obtained contribute to a deeper understanding of two-phase inner flows.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139595174","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Mosquitoes’ self-generated air movements around their antennae, especially at the wing-beat frequency, are crucial for both obstacle avoidance and mating communication. However, the characteristics of these air movements are not well clarified. In this study, the air movements induced by wing tones (sound generated by flapping wings in flight) around the antennae of a mosquito-like model (Culex quinquefasciatus, male) are investigated using the acoustic analogy method. Both the self-generated wing tone and the wing tone reflected from the ground are calculated. Given that the tiny changes in direction and magnitude of air movements can be detected by the mosquito’s antennae, a novel method is introduced to intuitively characterize the air movements induced by the wing tone. The air movements are decomposed into two basic modes (oscillation and revolution). Our results show that, without considering the scattering on the mosquito’s body, the self-generated sound wave of the wing-beat frequency around the antennae mainly induces air oscillation, with the velocity amplitude exceeding the mosquito’s hearing threshold of the male wingbeat frequency by two orders of magnitude. Moreover, when the model is positioned at a distance from the ground greater than approximately two wing lengths, the reflected sound wave at the male wingbeat frequency attenuates below the hearing threshold. That is, the role of reflected wing tone in the mosquito’s obstacle avoidance mechanism appears negligible. Our findings and method may provide insight into how mosquitoes avoid obstacles when their vision is unavailable and inspire the development of collision avoidance systems in micro-aerial vehicles.
{"title":"Characterization of the Wing Tone around the Antennae of a Mosquito-like Model","authors":"Yongtao Wang, Zhiteng Zhou, Zhuoyu Xie","doi":"10.3390/fluids9020031","DOIUrl":"https://doi.org/10.3390/fluids9020031","url":null,"abstract":"Mosquitoes’ self-generated air movements around their antennae, especially at the wing-beat frequency, are crucial for both obstacle avoidance and mating communication. However, the characteristics of these air movements are not well clarified. In this study, the air movements induced by wing tones (sound generated by flapping wings in flight) around the antennae of a mosquito-like model (Culex quinquefasciatus, male) are investigated using the acoustic analogy method. Both the self-generated wing tone and the wing tone reflected from the ground are calculated. Given that the tiny changes in direction and magnitude of air movements can be detected by the mosquito’s antennae, a novel method is introduced to intuitively characterize the air movements induced by the wing tone. The air movements are decomposed into two basic modes (oscillation and revolution). Our results show that, without considering the scattering on the mosquito’s body, the self-generated sound wave of the wing-beat frequency around the antennae mainly induces air oscillation, with the velocity amplitude exceeding the mosquito’s hearing threshold of the male wingbeat frequency by two orders of magnitude. Moreover, when the model is positioned at a distance from the ground greater than approximately two wing lengths, the reflected sound wave at the male wingbeat frequency attenuates below the hearing threshold. That is, the role of reflected wing tone in the mosquito’s obstacle avoidance mechanism appears negligible. Our findings and method may provide insight into how mosquitoes avoid obstacles when their vision is unavailable and inspire the development of collision avoidance systems in micro-aerial vehicles.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139599596","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
J. González-Trejo, R. Miranda-Tello, Ruslan Gabbasov, C. A. Real-Ramírez, F. Cervantes-De-La-Torre
This work studies how the sliding-gate valve (SGV) modifies the features and the dynamic behavior of the outlet jets for flat-bottom and well-bottom bifurcated submerged entry nozzles (SENs) used in continuous casting machines. Three conditions for the SGV were studied: no obstruction, moderate obstruction, and severe obstruction. The experimental study used a scaled model, employing cold water as the working fluid. A high-frequency analysis of the flow inside the SEN’s bore arriving at the outlet ports was performed by employing the particle image velocimetry (PIV) technique. Low-frequency measurements of the volumetric flow at the exit port were obtained by splitting the exit jet into four quadrants and employing digital flowmeters. It was observed that reducing the SGV clearance increases the turbulence of the flow inside the SEN bore, but the flow displays ordered rather than erratic fluctuations. Flowmeter measurements showed that, regardless of the level of obstruction in the SGV, the outlet jets on flat-bottom and the well-bottom SENs have dynamic behaviors and features with significant differences. This finding is relevant because the flow distribution inside the outlet ports is directly related to the jet’s wideness, affecting the recirculation pattern inside the mold and, therefore, the quality of the finished steel slab.
这项工作研究了滑动闸阀（SGV）如何改变连铸机中使用的平底和井底分叉浸没式入口喷嘴（SEN）的出口喷流特征和动态行为。研究了 SGV 的三种情况：无阻塞、中等阻塞和严重阻塞。实验研究使用了一个缩放模型，以冷水作为工作流体。利用粒子图像测速仪（PIV）技术对到达出口端口的 SEN 内孔流动情况进行了高频分析。通过将出口射流分成四个象限并使用数字流量计，对出口端口的体积流量进行了低频测量。测量结果表明，减小 SGV 间隙会增加 SEN 孔内的湍流，但湍流是有序的，而不是无规律的波动。流量计测量结果表明，无论 SGV 的阻塞程度如何，平底和井底 SEN 的出口喷流都具有显著差异的动态行为和特征。这一发现很有意义，因为出口端口内的流量分布与喷流的宽窄直接相关，会影响模具内的再循环模式，从而影响成品钢坯的质量。
{"title":"Experimental Analysis of the Influence of the Sliding-Gate Valve on Submerged Entry Nozzle Outlet Jets","authors":"J. González-Trejo, R. Miranda-Tello, Ruslan Gabbasov, C. A. Real-Ramírez, F. Cervantes-De-La-Torre","doi":"10.3390/fluids9010030","DOIUrl":"https://doi.org/10.3390/fluids9010030","url":null,"abstract":"This work studies how the sliding-gate valve (SGV) modifies the features and the dynamic behavior of the outlet jets for flat-bottom and well-bottom bifurcated submerged entry nozzles (SENs) used in continuous casting machines. Three conditions for the SGV were studied: no obstruction, moderate obstruction, and severe obstruction. The experimental study used a scaled model, employing cold water as the working fluid. A high-frequency analysis of the flow inside the SEN’s bore arriving at the outlet ports was performed by employing the particle image velocimetry (PIV) technique. Low-frequency measurements of the volumetric flow at the exit port were obtained by splitting the exit jet into four quadrants and employing digital flowmeters. It was observed that reducing the SGV clearance increases the turbulence of the flow inside the SEN bore, but the flow displays ordered rather than erratic fluctuations. Flowmeter measurements showed that, regardless of the level of obstruction in the SGV, the outlet jets on flat-bottom and the well-bottom SENs have dynamic behaviors and features with significant differences. This finding is relevant because the flow distribution inside the outlet ports is directly related to the jet’s wideness, affecting the recirculation pattern inside the mold and, therefore, the quality of the finished steel slab.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139610182","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The understanding of the boundary layer flame flashback (BLF) has considerably improved in recent decades, driven by the increasing focus on clean energy and the need to address the operational issues associated with flashback. This study investigates the influence of the Lewis number (Le) on symmetric flame shapes under the critical conditions for a laminar boundary layer flashback in cylindrical tubes. It has been found that the transformation of the flame shape from a mushroom to a tulip happens in a tube of a given radius, as the thermal expansion coefficient and Le are modified. A smaller Lewis number results in a local increase in the burning rate at the flame tip, with the flame being able to propagate closer to the wall, which significantly increases the flashback propensity, in line with previous findings. In cases with a Lewis number smaller than unity, a higher thermal expansion results in a flame propagation happening closer to the wall, thus facing a weaker oncoming flow and, consequently, becoming more prone to flashback. For Le > 1, the effect of the increase in the thermal expansion coefficient on the flashback tendency is much less pronounced.
近几十年来，人们对边界层火焰回火（BLF）的认识有了很大提高，其原因是对清洁能源的日益关注以及需要解决与回火相关的运行问题。本研究探讨了刘易斯数（Le）对圆柱形管内层流边界层回火临界条件下对称火焰形状的影响。研究发现，在给定半径的管子中，随着热膨胀系数和 Le 的变化，火焰形状会从蘑菇形转变为郁金香形。路易斯数越小，火焰顶端的局部燃烧率就越高，火焰可以传播到更靠近管壁的地方，这就大大增加了回火的可能性，这与之前的研究结果是一致的。在路易斯数小于 1 的情况下，较高的热膨胀会导致火焰传播更靠近墙壁，从而面对较弱的迎气流，因此更容易发生回火。当 Le > 1 时，热膨胀系数的增加对回火倾向的影响要小得多。
{"title":"Combined Impact of the Lewis Number and Thermal Expansion on Laminar Flame Flashback in Tubes","authors":"Kai Huang, Louis Benteux, Wenhu Han, D. Valiev","doi":"10.3390/fluids9010028","DOIUrl":"https://doi.org/10.3390/fluids9010028","url":null,"abstract":"The understanding of the boundary layer flame flashback (BLF) has considerably improved in recent decades, driven by the increasing focus on clean energy and the need to address the operational issues associated with flashback. This study investigates the influence of the Lewis number (Le) on symmetric flame shapes under the critical conditions for a laminar boundary layer flashback in cylindrical tubes. It has been found that the transformation of the flame shape from a mushroom to a tulip happens in a tube of a given radius, as the thermal expansion coefficient and Le are modified. A smaller Lewis number results in a local increase in the burning rate at the flame tip, with the flame being able to propagate closer to the wall, which significantly increases the flashback propensity, in line with previous findings. In cases with a Lewis number smaller than unity, a higher thermal expansion results in a flame propagation happening closer to the wall, thus facing a weaker oncoming flow and, consequently, becoming more prone to flashback. For Le > 1, the effect of the increase in the thermal expansion coefficient on the flashback tendency is much less pronounced.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139612475","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
We present an analysis of the effect of particle inertia and thermal inertia on the heat transfer in a turbulent shearless flow, where an inhomogeneous passive temperature field is advected along with inertial point particles by a homogeneous isotropic velocity field. Eulerian–Lagrangian direct numerical simulations are carried out in both one- and two-way coupling regimes and analyzed through single-point statistics. The role of particle inertia and thermal inertia is discussed by introducing a new decomposition of particle second-order moments in terms of correlations involving Lagrangian acceleration and time derivative of particles. We present how particle relaxation times mediate the level of particle velocity–temperature correlation, which gives particle contribution to the overall heat transfer. For each thermal Stokes number, a critical Stokes number is individuated. The effect of particle feedback on the attenuation or enhancement of fluid temperature variance is presented. We show that particle feedback enhances fluid temperature variance for Stokes numbers less than one and damps is for larger than one Stokes number, regardless of the thermal Stokes number, even if this effect is amplified by an increasing thermal inertia.
{"title":"The Role of Particle Inertia and Thermal Inertia in Heat Transfer in a Non-Isothermal Particle-Laden Turbulent Flow","authors":"Hamid Reza Zandi Pour, Michele Iovieno","doi":"10.3390/fluids9010029","DOIUrl":"https://doi.org/10.3390/fluids9010029","url":null,"abstract":"We present an analysis of the effect of particle inertia and thermal inertia on the heat transfer in a turbulent shearless flow, where an inhomogeneous passive temperature field is advected along with inertial point particles by a homogeneous isotropic velocity field. Eulerian–Lagrangian direct numerical simulations are carried out in both one- and two-way coupling regimes and analyzed through single-point statistics. The role of particle inertia and thermal inertia is discussed by introducing a new decomposition of particle second-order moments in terms of correlations involving Lagrangian acceleration and time derivative of particles. We present how particle relaxation times mediate the level of particle velocity–temperature correlation, which gives particle contribution to the overall heat transfer. For each thermal Stokes number, a critical Stokes number is individuated. The effect of particle feedback on the attenuation or enhancement of fluid temperature variance is presented. We show that particle feedback enhances fluid temperature variance for Stokes numbers less than one and damps is for larger than one Stokes number, regardless of the thermal Stokes number, even if this effect is amplified by an increasing thermal inertia.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139525746","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Leonardo Geronzi, B. Fanni, Bart De Jong, Gerben Roest, Sasa Kenjeres, S. Celi, M. E. Biancolini
The treatment for asthma and chronic obstructive pulmonary disease relies on forced inhalation of drug particles. Their distribution is essential for maximizing the outcomes. Patient-specific computational fluid dynamics (CFD) simulations can be used to optimize these therapies. In this regard, this study focuses on creating a parametric model of the human respiratory tract from which synthetic anatomies for particle deposition analysis through CFD simulation could be derived. A baseline geometry up to the fourth generation of bronchioles was extracted from a CT dataset. Radial basis function (RBF) mesh morphing acting on a dedicated tree structure was used to modify this baseline mesh, extracting 1000 synthetic anatomies. A total of 26 geometrical parameters affecting branch lengths, angles, and diameters were controlled. Morphed models underwent CFD simulations to analyze airflow and particle dynamics. Mesh morphing was crucial in generating high-quality computational grids, with 96% of the synthetic database being immediately suitable for accurate CFD simulations. Variations in wall shear stress, particle accretion rate, and turbulent kinetic energy across different anatomies highlighted the impact of the anatomical shape on drug delivery and deposition. The study successfully demonstrates the potential of tree-structure-based RBF mesh morphing in generating parametric airways for drug delivery studies.
{"title":"A Parametric 3D Model of Human Airways for Particle Drug Delivery and Deposition","authors":"Leonardo Geronzi, B. Fanni, Bart De Jong, Gerben Roest, Sasa Kenjeres, S. Celi, M. E. Biancolini","doi":"10.3390/fluids9010027","DOIUrl":"https://doi.org/10.3390/fluids9010027","url":null,"abstract":"The treatment for asthma and chronic obstructive pulmonary disease relies on forced inhalation of drug particles. Their distribution is essential for maximizing the outcomes. Patient-specific computational fluid dynamics (CFD) simulations can be used to optimize these therapies. In this regard, this study focuses on creating a parametric model of the human respiratory tract from which synthetic anatomies for particle deposition analysis through CFD simulation could be derived. A baseline geometry up to the fourth generation of bronchioles was extracted from a CT dataset. Radial basis function (RBF) mesh morphing acting on a dedicated tree structure was used to modify this baseline mesh, extracting 1000 synthetic anatomies. A total of 26 geometrical parameters affecting branch lengths, angles, and diameters were controlled. Morphed models underwent CFD simulations to analyze airflow and particle dynamics. Mesh morphing was crucial in generating high-quality computational grids, with 96% of the synthetic database being immediately suitable for accurate CFD simulations. Variations in wall shear stress, particle accretion rate, and turbulent kinetic energy across different anatomies highlighted the impact of the anatomical shape on drug delivery and deposition. The study successfully demonstrates the potential of tree-structure-based RBF mesh morphing in generating parametric airways for drug delivery studies.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139525929","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pavel Bulat, Pavel Chernyshov, Nikolay Prodan, Konstantin Volkov
The article explores flow behavior around thick airfoils at low Reynolds numbers and the potential application of energy methods to manipulate the flow field for increased lift and reduced drag. The study relies on a set of propulsion airfoils calculated using a combined approach of solving the inverse problem of aerodynamics and applying stochastic global optimization methods. The calculations consider the transition from laminar to turbulent flow regimes, which significantly affects lift and airfoil drag. The suitability of different turbulence models for airfoil modeling in low Reynolds numbers is discussed, and numerical simulation results determine the lift coefficient dependence on angle of attack and the optimal air flow rate taken from the airfoil surface for each angle of attack. The accuracy of different turbulence models is analyzed by comparing numerical simulation results to physical experiment data.
{"title":"Control of Aerodynamic Characteristics of Thick Airfoils at Low Reynolds Numbers Using Methods of Boundary Layer Control","authors":"Pavel Bulat, Pavel Chernyshov, Nikolay Prodan, Konstantin Volkov","doi":"10.3390/fluids9010026","DOIUrl":"https://doi.org/10.3390/fluids9010026","url":null,"abstract":"The article explores flow behavior around thick airfoils at low Reynolds numbers and the potential application of energy methods to manipulate the flow field for increased lift and reduced drag. The study relies on a set of propulsion airfoils calculated using a combined approach of solving the inverse problem of aerodynamics and applying stochastic global optimization methods. The calculations consider the transition from laminar to turbulent flow regimes, which significantly affects lift and airfoil drag. The suitability of different turbulence models for airfoil modeling in low Reynolds numbers is discussed, and numerical simulation results determine the lift coefficient dependence on angle of attack and the optimal air flow rate taken from the airfoil surface for each angle of attack. The accuracy of different turbulence models is analyzed by comparing numerical simulation results to physical experiment data.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139527366","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Vincent Proulx-Cabana, Guilhem Michon, Eric Laurendeau
The aim of this article is to investigate the parameter sensitivity of the (Non-Linear) Unsteady Vortex Lattice Method-Vortex Particle Method [(NL-)UVLM-VPM] with Particle Strength Exchange-Large Eddy Simulations (PSE-LES) method on a lower Reynolds number rotor. The previous work detailed the method, but introduced parameters whose influence were not investigated. Most importantly, the Vreman model coefficient was chosen arbitrarily and was not suitable to ensure stability for this lower Reynolds number rotor simulation. In addition, the previous work presented a consistency study where geometry and time discretization were refined simultaneously. The present article starts with a comparative literature review of potential methods used to solve the aerodynamics of an isolated hovering rotor. This review highlights the differences in modeling, discretizations, sensitivity analysis, validation cases, and the results chosen by the different studies. Then, a transparent and thorough parametric study of the method is presented alongside discussions of the observed results and their physical interpretation regarding the flow. The sensitivity analysis is performed for the three free parameters of UVLM, namely Vatistas core size, the geometry and the temporal discretizations, and then for the three additional parameters introduced by UVLM-VPM, which are the Vreman model coefficient, the particle spacing, and the conversion time. The effect of different databases in the non-linear coupling is also shown. The method is shown to be consistent with both geometry and temporal refinements. It is also consistent with the expected behavior of the different parameters change, including the numerical stability that depends on the strength of the LES diffusion controlled by the Vreman model coefficient. The effect of discretization refinement presented here not only shows the integrated coefficients where different errors can cancel each other, but also looks at their convergence and where relevant, the distributed loads and tip singularity position. Finally, the aerodynamics results of the method are compared for different databases and with higher fidelity Unsteady Reynolds Averaged Navier–Stokes (URANS) 3D results on a lower Reynolds number rotor.
{"title":"Parametrization Effects of the Non-Linear Unsteady Vortex Method with Vortex Particle Method for Small Rotor Aerodynamics","authors":"Vincent Proulx-Cabana, Guilhem Michon, Eric Laurendeau","doi":"10.3390/fluids9010024","DOIUrl":"https://doi.org/10.3390/fluids9010024","url":null,"abstract":"The aim of this article is to investigate the parameter sensitivity of the (Non-Linear) Unsteady Vortex Lattice Method-Vortex Particle Method [(NL-)UVLM-VPM] with Particle Strength Exchange-Large Eddy Simulations (PSE-LES) method on a lower Reynolds number rotor. The previous work detailed the method, but introduced parameters whose influence were not investigated. Most importantly, the Vreman model coefficient was chosen arbitrarily and was not suitable to ensure stability for this lower Reynolds number rotor simulation. In addition, the previous work presented a consistency study where geometry and time discretization were refined simultaneously. The present article starts with a comparative literature review of potential methods used to solve the aerodynamics of an isolated hovering rotor. This review highlights the differences in modeling, discretizations, sensitivity analysis, validation cases, and the results chosen by the different studies. Then, a transparent and thorough parametric study of the method is presented alongside discussions of the observed results and their physical interpretation regarding the flow. The sensitivity analysis is performed for the three free parameters of UVLM, namely Vatistas core size, the geometry and the temporal discretizations, and then for the three additional parameters introduced by UVLM-VPM, which are the Vreman model coefficient, the particle spacing, and the conversion time. The effect of different databases in the non-linear coupling is also shown. The method is shown to be consistent with both geometry and temporal refinements. It is also consistent with the expected behavior of the different parameters change, including the numerical stability that depends on the strength of the LES diffusion controlled by the Vreman model coefficient. The effect of discretization refinement presented here not only shows the integrated coefficients where different errors can cancel each other, but also looks at their convergence and where relevant, the distributed loads and tip singularity position. Finally, the aerodynamics results of the method are compared for different databases and with higher fidelity Unsteady Reynolds Averaged Navier–Stokes (URANS) 3D results on a lower Reynolds number rotor.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139621960","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The usual concept of simulation methods for turbulent flows is to impose a certain (partial) flow resolution. This concept becomes problematic away from limit regimes of no or an almost complete flow resolution: discrepancies between the imposed and actual flow resolution may imply an unreliable model behavior and high computational cost to compensate for simulation deficiencies. An exact mathematical approach based on variational analysis provides a solution to these problems. Minimal error continuous eddy simulation (CES) designed in this way enables simulations in which the model actively responds to variations in flow resolution by increasing or decreasing its contribution to the simulation as required. This paper presents the first application of CES methods to a moderately complex, relatively high Reynolds number turbulent flow simulation: the NASA wall-mounted hump flow. It is shown that CES performs equally well or better than almost resolving simulation methods at a little fraction of computational cost. Significant computational cost and performance advantages are reported in comparison to popular partially resolving simulation methods including detached eddy simulation and wall-modeled large eddy simulation. Characteristic features of the asymptotic flow structure are identified on the basis of CES simulations.
湍流模拟方法的通常概念是施加一定的（部分）流动分辨率。在没有或几乎完全没有流动分辨率的极限状态下，这一概念就会出现问题：强加的流动分辨率与实际流动分辨率之间的差异可能意味着不可靠的模型行为和弥补模拟缺陷的高计算成本。基于变分分析的精确数学方法为这些问题提供了解决方案。采用这种方法设计的最小误差连续涡模拟（CES）可以在模拟过程中使模型主动响应流动分辨率的变化，根据需要增加或减少其对模拟的贡献。本文首次将 CES 方法应用于中度复杂、雷诺数相对较高的湍流模拟：NASA 壁挂式驼峰流。结果表明，CES 的性能与解析模拟方法相当，甚至更好，而计算成本仅为后者的一小部分。与流行的部分解析模拟方法（包括分离涡模拟和壁面建模大涡模拟）相比，CES 在计算成本和性能方面具有显著优势。在 CES 模拟的基础上，确定了渐近流动结构的特征。
{"title":"Continuous Eddy Simulation vs. Resolution-Imposing Simulation Methods for Turbulent Flows","authors":"A. Fagbade, Stefan Heinz","doi":"10.3390/fluids9010022","DOIUrl":"https://doi.org/10.3390/fluids9010022","url":null,"abstract":"The usual concept of simulation methods for turbulent flows is to impose a certain (partial) flow resolution. This concept becomes problematic away from limit regimes of no or an almost complete flow resolution: discrepancies between the imposed and actual flow resolution may imply an unreliable model behavior and high computational cost to compensate for simulation deficiencies. An exact mathematical approach based on variational analysis provides a solution to these problems. Minimal error continuous eddy simulation (CES) designed in this way enables simulations in which the model actively responds to variations in flow resolution by increasing or decreasing its contribution to the simulation as required. This paper presents the first application of CES methods to a moderately complex, relatively high Reynolds number turbulent flow simulation: the NASA wall-mounted hump flow. It is shown that CES performs equally well or better than almost resolving simulation methods at a little fraction of computational cost. Significant computational cost and performance advantages are reported in comparison to popular partially resolving simulation methods including detached eddy simulation and wall-modeled large eddy simulation. Characteristic features of the asymptotic flow structure are identified on the basis of CES simulations.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139441110","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Some non-Newtonian fluids have time-dependent rheological properties like a shear stress that depends on the shear history or a stress overshoot that is a function of the resting time, when fluid movement is started. The rheological properties of such complex fluids may not stay constant while they are used in an industrial process, and it is therefore desirable to measure these properties frequently and with a simple and robust device like a pipe rheometer. This paper investigated how the time-dependent rheological properties of a thixotropic and viscoelastic shear-thinning fluid can be extracted from differential pressure measurements obtained at different flowrates along a circular pipe section. The method consists in modeling the flow of a thixotropic version of a Quemada fluid and solving the inverse problem in order to find the model parameters using the measurements made in steady-state conditions. Also, a Maxwell linear viscoelastic model was used to reproduce the stress overshoot observed when starting circulation after a resting period. The pipe rheometer was designed to have the proper features necessary to exhibit the thixotropic and viscoelastic effects that were needed to calibrate the rheological model parameters. The accuracy of rheological measurements depends on understanding the effects that can influence the observations and on a proper design that takes advantage of these side effects instead of attempting to eliminate them.
{"title":"Calibration of Thixotropic and Viscoelastic Shear-Thinning Fluids Using Pipe Rheometer Measurements","authors":"E. Cayeux","doi":"10.3390/fluids9010023","DOIUrl":"https://doi.org/10.3390/fluids9010023","url":null,"abstract":"Some non-Newtonian fluids have time-dependent rheological properties like a shear stress that depends on the shear history or a stress overshoot that is a function of the resting time, when fluid movement is started. The rheological properties of such complex fluids may not stay constant while they are used in an industrial process, and it is therefore desirable to measure these properties frequently and with a simple and robust device like a pipe rheometer. This paper investigated how the time-dependent rheological properties of a thixotropic and viscoelastic shear-thinning fluid can be extracted from differential pressure measurements obtained at different flowrates along a circular pipe section. The method consists in modeling the flow of a thixotropic version of a Quemada fluid and solving the inverse problem in order to find the model parameters using the measurements made in steady-state conditions. Also, a Maxwell linear viscoelastic model was used to reproduce the stress overshoot observed when starting circulation after a resting period. The pipe rheometer was designed to have the proper features necessary to exhibit the thixotropic and viscoelastic effects that were needed to calibrate the rheological model parameters. The accuracy of rheological measurements depends on understanding the effects that can influence the observations and on a proper design that takes advantage of these side effects instead of attempting to eliminate them.","PeriodicalId":12397,"journal":{"name":"Fluids","volume":null,"pages":null},"PeriodicalIF":1.9,"publicationDate":"2024-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139534673","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：