首页 > 最新文献

Physical Review Fluids最新文献

英文 中文
Investigating the effect of turbulence on hemolysis through cell-resolved fluid-structure interaction simulations of individual red blood cells 通过对单个红细胞进行细胞分辨流固耦合模拟,研究湍流对溶血的影响
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-26 DOI: 10.1103/physrevfluids.9.073102
Grant Rydquist, Mahdi Esmaily
Existing hemolysis algorithms are often constructed for laminar flows that expose red blood cells (RBCs) to a constant rate of shear. It remains an open question whether such models are applicable to turbulent flows, where there is a significant variation in shear rate along cell trajectories. To evaluate the effect of turbulence on hemolysis, we perform cell-resolved simulations of isolated RBCs in turbulent channel flow at Reτ=180 and 360 and compare them against the results obtained from laminar flow simulations at an equivalent wall shear stress. The RBCs are modeled as isolated cells in an unbounded domain with the viscosity of the bulk fluid used for the surrounding fluid. This comparison shows that, while the laminar flow generally induces greater stretch in the cell in a time-averaged sense, cells experience an overall larger deformation in turbulence. This difference is attributed to extreme events in turbulence that occasionally create bursts of high shear conditions, which, consequently, induce a large deformation in the cells. Associating damage with the most extreme deformation regimes, we observe that, in the worst case, the turbulent flow can produce deformation in the cell that is higher than the absolute maximum value in the analogous laminar case approximately 14% of the time. Additionally, the Reτ=180 universally induced greater deformation in the cells than the Reτ=360 case, suggesting that increasing the range of scales in the flow does not necessarily yield greater deformation when all other parameters are kept constant. A strong direct correlation (R>0.8) between shear rate and deformation metrics was observed in turbulence. The correlation against Q-criterion is inverse and weaker (R0.26), but once the shear contribution is subtracted, it improves in terms of areal dilatation (R0.6).
现有的溶血算法通常是针对层流设计的,层流会使红细胞(RBC)受到恒定的剪切力作用。此类模型是否适用于湍流仍是一个未决问题,因为在湍流中,细胞轨迹上的剪切率会发生显著变化。为了评估湍流对溶血的影响,我们在 Reτ=180 和 360 条件下对湍流通道流中的孤立 RBC 进行了细胞分辨模拟,并与等效壁面剪应力下的层流模拟结果进行了比较。RBC 被模拟为无界域中的孤立细胞,周围流体的粘度采用散装流体的粘度。比较结果表明,层流通常会在细胞中引起更大的时间平均意义上的拉伸,而在湍流中细胞会经历更大的整体变形。这种差异可归因于湍流中的极端事件,这些事件偶尔会产生突发的高剪切力,从而导致细胞产生较大变形。将破坏与最极端的变形状态联系起来,我们发现,在最糟糕的情况下,湍流在大约 14% 的时间内产生的单元变形高于类似层流情况下的绝对最大值。此外,与 Reτ=360 的情况相比,Reτ=180 在所有情况下都会引起细胞更大的变形,这表明在所有其他参数保持不变的情况下,增加流动的尺度范围并不一定会产生更大的变形。在湍流中观察到剪切速率与变形指标之间存在很强的直接相关性(R>0.8)。与 Q 值标准的相关性较弱(R≈-0.26),但一旦减去剪切力的贡献,则在面积扩张方面有所改善(R≈-0.6)。
{"title":"Investigating the effect of turbulence on hemolysis through cell-resolved fluid-structure interaction simulations of individual red blood cells","authors":"Grant Rydquist, Mahdi Esmaily","doi":"10.1103/physrevfluids.9.073102","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.073102","url":null,"abstract":"Existing hemolysis algorithms are often constructed for laminar flows that expose red blood cells (RBCs) to a constant rate of shear. It remains an open question whether such models are applicable to turbulent flows, where there is a significant variation in shear rate along cell trajectories. To evaluate the effect of turbulence on hemolysis, we perform cell-resolved simulations of isolated RBCs in turbulent channel flow at <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msub><mtext>Re</mtext><mi>τ</mi></msub><mo>=</mo><mn>180</mn></mrow></math> and 360 and compare them against the results obtained from laminar flow simulations at an equivalent wall shear stress. The RBCs are modeled as isolated cells in an unbounded domain with the viscosity of the bulk fluid used for the surrounding fluid. This comparison shows that, while the laminar flow generally induces greater stretch in the cell in a time-averaged sense, cells experience an overall larger deformation in turbulence. This difference is attributed to extreme events in turbulence that occasionally create bursts of high shear conditions, which, consequently, induce a large deformation in the cells. Associating damage with the most extreme deformation regimes, we observe that, in the worst case, the turbulent flow can produce deformation in the cell that is higher than the absolute maximum value in the analogous laminar case approximately <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mn>14</mn><mo>%</mo></mrow></math> of the time. Additionally, the <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msub><mi>Re</mi><mi>τ</mi></msub><mo>=</mo><mn>180</mn></mrow></math> universally induced greater deformation in the cells than the <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msub><mi>Re</mi><mi>τ</mi></msub><mo>=</mo><mn>360</mn></mrow></math> case, suggesting that increasing the range of scales in the flow does not necessarily yield greater deformation when all other parameters are kept constant. A strong direct correlation (<math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi>R</mi><mo>&gt;</mo><mn>0.8</mn></mrow></math>) between shear rate and deformation metrics was observed in turbulence. The correlation against <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>Q</mi></math>-criterion is inverse and weaker (<math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi>R</mi><mo>≈</mo><mo>−</mo><mn>0.26</mn></mrow></math>), but once the shear contribution is subtracted, it improves in terms of areal dilatation (<math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi>R</mi><mo>≈</mo><mo>−</mo><mn>0.6</mn></mrow></math>).","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"9 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141785246","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}
引用次数: 0
Evolution of turbulence using a random jet array 使用随机射流阵列的湍流演变
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-26 DOI: 10.1103/physrevfluids.9.074610
Arefe Ghazi Nezami, Blair Anne Johnson
Random jet arrays (RJAs) have been shown to be effective in generating zero mean flow homogeneous isotropic turbulence. While many laboratory studies have investigated the flow in these facilities, there are several remaining questions regarding the evolution of turbulence, from the development of turbulence to where it decays, along with understanding how input energy from the jet array transfers into different turbulent flow characteristics. To address these questions, we perform a series of laboratory experiments in which we alter the parameters of the randomized algorithm, along with the jet spacing and outlet velocity of the jets. We first determine the location where turbulence transitions to a fully developed state and show that it is a function of jet penetration length, LJ, and effective jet spacing, Se. We identify three distinct regions for the spatial decay of turbulence in RJA facilities and notably, we find different decay rates, unlike previous studies that report only one spatial decay rate using similar facilities. These regions are shown to depend on the variations of input parameters yet independent of the strength of the mean flow. We also find the strength of the mean flow does not affect the homogeneity, nor the production, transport, or advection terms of the turbulent kinetic energy budget equation. Finally, we address a longstanding question toward estimating turbulence metrics with an RJA based on the input parameters. We define an efficiency parameter that provides insight into the transfer rate of input power to the dissipation rate of the generated turbulence.
随机射流阵列(RJAs)已被证明可以有效地产生零平均流均质各向同性湍流。虽然许多实验室研究已经对这些设施中的流动进行了调查,但关于湍流的演变,从湍流的发展到湍流的衰减,以及了解射流阵列的输入能量如何转移到不同的湍流特性等,还存在一些问题。为了解决这些问题,我们进行了一系列实验室实验,在实验中我们改变了随机算法的参数以及射流间距和射流出口速度。我们首先确定了湍流过渡到充分发展状态的位置,并证明它是射流穿透长度 LJ 和有效射流间距 Se 的函数。我们确定了 RJA 设备中湍流空间衰减的三个不同区域,值得注意的是,我们发现了不同的衰减速率,这与之前使用类似设备只报告一种空间衰减速率的研究不同。研究表明,这些区域取决于输入参数的变化,但与平均流的强度无关。我们还发现,平均流的强度不会影响均质性,也不会影响湍流动能预算方程中的生成、传输或平流项。最后,我们解决了一个长期存在的问题,即根据输入参数用 RJA 估算湍流度量。我们定义了一个效率参数,该参数可帮助我们了解输入功率向所产生湍流的耗散率的转移率。
{"title":"Evolution of turbulence using a random jet array","authors":"Arefe Ghazi Nezami, Blair Anne Johnson","doi":"10.1103/physrevfluids.9.074610","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074610","url":null,"abstract":"Random jet arrays (RJAs) have been shown to be effective in generating zero mean flow homogeneous isotropic turbulence. While many laboratory studies have investigated the flow in these facilities, there are several remaining questions regarding the evolution of turbulence, from the development of turbulence to where it decays, along with understanding how input energy from the jet array transfers into different turbulent flow characteristics. To address these questions, we perform a series of laboratory experiments in which we alter the parameters of the randomized algorithm, along with the jet spacing and outlet velocity of the jets. We first determine the location where turbulence transitions to a fully developed state and show that it is a function of jet penetration length, <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi mathvariant=\"script\">L</mi><mi mathvariant=\"script\">J</mi></msub></math>, and effective jet spacing, <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi>S</mi><mi>e</mi></msub></math>. We identify three distinct regions for the spatial decay of turbulence in RJA facilities and notably, we find different decay rates, unlike previous studies that report only one spatial decay rate using similar facilities. These regions are shown to depend on the variations of input parameters yet independent of the strength of the mean flow. We also find the strength of the mean flow does not affect the homogeneity, nor the production, transport, or advection terms of the turbulent kinetic energy budget equation. Finally, we address a longstanding question toward estimating turbulence metrics with an RJA based on the input parameters. We define an efficiency parameter that provides insight into the transfer rate of input power to the dissipation rate of the generated turbulence.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"31 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141785441","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}
引用次数: 0
Dynamics of a magnetic particle in an oscillating magnetic field 振荡磁场中的磁性粒子动力学
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-26 DOI: 10.1103/physrevfluids.9.074303
I. Misra, V. Kumaran
The orientation dynamics of and the torque fluctuations due to a spheroidal magnetic particle in an oscillating magnetic field are analyzed in the Stokes flow regime. For a permanent dipole, the dynamics depends on ω, the ratio of the magnetic field frequency, and the viscous relaxation rate. For ω1, the particle executes oscillations with amplitude (ω)1 about its initial orientation. The average torque is zero because the particle does not execute complete rotations, and the root mean square of the torque fluctuations scaled by the characteristic magnetic torque tends to a constant in this limit. For ω1, the orientation is close to the magnetic field direction for most of the oscillation period, and it rapidly rotates when the field passes through extrema. The scaled root mean square of the torque fluctuations is proportional to (ω)1/2 in this limit. The particle orientation aligns along the magnetic field direction for different models of induced dipoles if the magnetization is nonhysteretic. For the hysteretic Stoner-Wohlfarth model, the dynamics also depends on the parameter h0, the ratio of the Zeeman energy, and the anisotropy energy. For h01, the magnetic moment oscillates about one pole of the orientation vector, and the orientation vector rapidly rotates when the field passes through extrema in a manner similar to that for a permanent dipole. For h01, the magnetic moment switches between the two poles of the orientation vector, and the orientation vector executes small amplitude oscillations about the field direction. There is a discontinuous transition between the oscillating and switching magnetic moment which depends on h0 and the initial orientation.
在斯托克斯流制度下,分析了振荡磁场中球形磁性粒子的取向动力学和扭矩波动。对于永久偶极子,动态取决于ω†、磁场频率比和粘性弛豫率。对于 ω†≫1,粒子以振幅 ∼(ω†)-1 围绕其初始方向摆动。由于粒子不会进行完全旋转,因此平均转矩为零,在此极限下,转矩波动的均方根与特征磁转矩的比例趋于恒定。对于ω†≪1,在振荡周期的大部分时间里,方向接近磁场方向,当磁场通过极值时,方向迅速旋转。在此极限下,转矩波动的均方根比例与 (ω†)-1/2 成正比。如果磁化是非滞后的,那么对于不同的感应偶极子模型,粒子方向会沿着磁场方向排列。对于滞磁性斯通纳-沃尔法特模型,动力学还取决于参数 h0、泽曼能和各向异性能的比值。当 h0≪1 时,磁矩围绕取向矢量的一个磁极摆动,当磁场通过极值时,取向矢量迅速旋转,其方式与永久偶极类似。当 h0≫1 时,磁矩在取向矢量的两极之间切换,取向矢量围绕磁场方向做小振幅振荡。磁矩振荡和切换之间存在不连续的过渡,这取决于 h0 和初始方位。
{"title":"Dynamics of a magnetic particle in an oscillating magnetic field","authors":"I. Misra, V. Kumaran","doi":"10.1103/physrevfluids.9.074303","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074303","url":null,"abstract":"The orientation dynamics of and the torque fluctuations due to a spheroidal magnetic particle in an oscillating magnetic field are analyzed in the Stokes flow regime. For a permanent dipole, the dynamics depends on <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msup><mi>ω</mi><mo>†</mo></msup></math>, the ratio of the magnetic field frequency, and the viscous relaxation rate. For <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msup><mi>ω</mi><mo>†</mo></msup><mo>≫</mo><mn>1</mn></mrow></math>, the particle executes oscillations with amplitude <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mo>∼</mo><msup><mrow><mo>(</mo><msup><mi>ω</mi><mo>†</mo></msup><mo>)</mo></mrow><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow></math> about its initial orientation. The average torque is zero because the particle does not execute complete rotations, and the root mean square of the torque fluctuations scaled by the characteristic magnetic torque tends to a constant in this limit. For <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msup><mi>ω</mi><mo>†</mo></msup><mo>≪</mo><mn>1</mn></mrow></math>, the orientation is close to the magnetic field direction for most of the oscillation period, and it rapidly rotates when the field passes through extrema. The scaled root mean square of the torque fluctuations is proportional to <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msup><mrow><mo>(</mo><msup><mi>ω</mi><mo>†</mo></msup><mo>)</mo></mrow><mrow><mo>−</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msup></math> in this limit. The particle orientation aligns along the magnetic field direction for different models of induced dipoles if the magnetization is nonhysteretic. For the hysteretic Stoner-Wohlfarth model, the dynamics also depends on the parameter <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi>h</mi><mn>0</mn></msub></math>, the ratio of the Zeeman energy, and the anisotropy energy. For <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msub><mi>h</mi><mn>0</mn></msub><mo>≪</mo><mn>1</mn></mrow></math>, the magnetic moment oscillates about one pole of the orientation vector, and the orientation vector rapidly rotates when the field passes through extrema in a manner similar to that for a permanent dipole. For <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><msub><mi>h</mi><mn>0</mn></msub><mo>≫</mo><mn>1</mn></mrow></math>, the magnetic moment switches between the two poles of the orientation vector, and the orientation vector executes small amplitude oscillations about the field direction. There is a discontinuous transition between the oscillating and switching magnetic moment which depends on <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi>h</mi><mn>0</mn></msub></math> and the initial orientation.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"44 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775932","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}
引用次数: 0
Shape reconfiguration for underwater propeller efficiency improvement 重新配置形状以提高水下螺旋桨效率
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-26 DOI: 10.1103/physrevfluids.9.074402
Tristan Aurégan, Sylvain Courrech du Pont, Benjamin Thiria
We investigate experimentally the propulsive efficiency of a propeller in water with chordwise flexible blades that deform under the action of fluid loading. Using a scale model experiment, we record the deformation of the blades as well as the thrust and torque generated by the rotor. The use of flexible materials can improve the resilience to changing external conditions: with optimal flexibility, the blades deform and remain efficient under off-design conditions. We derive a theoretical law for blade tip deformation and show good agreement with experiments. Our results suggest that, using only the blade flexibility alone, we are able to program the blade deformation to passively adopt an optimized shape for efficient propulsion within a given parameter range.
我们通过实验研究了带弦向柔性叶片的螺旋桨在水中的推进效率,这种叶片在流体载荷的作用下会发生变形。通过比例模型实验,我们记录了叶片的变形以及转子产生的推力和扭矩。使用柔性材料可以提高对外部条件变化的适应能力:在最佳柔性条件下,叶片可以变形,并在非设计条件下保持高效。我们推导出了叶尖变形的理论规律,并与实验结果表明了良好的一致性。我们的研究结果表明,仅利用叶片的柔韧性,我们就能对叶片变形进行编程,使其在给定的参数范围内被动地采用优化形状,从而实现高效推进。
{"title":"Shape reconfiguration for underwater propeller efficiency improvement","authors":"Tristan Aurégan, Sylvain Courrech du Pont, Benjamin Thiria","doi":"10.1103/physrevfluids.9.074402","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074402","url":null,"abstract":"We investigate experimentally the propulsive efficiency of a propeller in water with chordwise flexible blades that deform under the action of fluid loading. Using a scale model experiment, we record the deformation of the blades as well as the thrust and torque generated by the rotor. The use of flexible materials can improve the resilience to changing external conditions: with optimal flexibility, the blades deform and remain efficient under off-design conditions. We derive a theoretical law for blade tip deformation and show good agreement with experiments. Our results suggest that, using only the blade flexibility alone, we are able to program the blade deformation to passively adopt an optimized shape for efficient propulsion within a given parameter range.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"51 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775935","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}
引用次数: 0
Rupture of a surfactant-laden draining thin film 含有表面活性剂的排水薄膜的破裂
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-25 DOI: 10.1103/physrevfluids.9.074004
Atul S. Vivek, Ranabir Dey, Harish N. Dixit
Surfactant-laden thin liquid films overlaid on solid substrates are encountered in a variety of industrial and biological settings. As these films reach submicron thickness, they tend to become unstable owing to the influence of long-range dispersion forces. In the current study, we investigate how gravitational drainage affects the stability attributes of such thin liquid films. Using scaling arguments, we demonstrate that gravity and dispersion forces can exert their influence simultaneously over a wide range of film thicknesses. In the lubrication limit, we carry out linear stability analysis and nonlinear simulations to understand the evolution of draining thin films. Linear stability indicates the existence of two unstable modes and two cutoff wave numbers, as opposed to a single unstable mode and a unique cutoff wave number observed in stationary films. It is also found that surfactant-laden flowing films are more stable than stationary films with surfactants as well as draining films with clean interfaces. The origin of stabilization is identified as the enhanced surfactant perturbations generated due to drainage. We demonstrate that films exhibiting intermediate levels of surfactant activity and significant drainage exhibit the lowest rates of disturbance growth, leading to extending the time of rupture.
在各种工业和生物环境中都会遇到覆盖在固体基底上的含有表面活性剂的液体薄膜。当这些薄膜厚度达到亚微米级时,由于长程分散力的影响,它们往往会变得不稳定。在当前的研究中,我们探讨了重力排水如何影响此类液体薄膜的稳定性属性。我们利用缩放论证证明,重力和分散力可以在很宽的薄膜厚度范围内同时施加影响。在润滑极限下,我们进行了线性稳定性分析和非线性模拟,以了解排水薄膜的演变过程。线性稳定性表明存在两个不稳定模式和两个截止波数,而静止薄膜中只有一个不稳定模式和一个唯一的截止波数。研究还发现,含有表面活性剂的流动薄膜比含有表面活性剂的静止薄膜以及具有清洁界面的排水薄膜更稳定。稳定的原因是由于排水产生的表面活性剂扰动增强。我们证明,表现出中等水平表面活性剂活性和显著排水的薄膜的扰动增长率最低,从而延长了破裂时间。
{"title":"Rupture of a surfactant-laden draining thin film","authors":"Atul S. Vivek, Ranabir Dey, Harish N. Dixit","doi":"10.1103/physrevfluids.9.074004","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074004","url":null,"abstract":"Surfactant-laden thin liquid films overlaid on solid substrates are encountered in a variety of industrial and biological settings. As these films reach submicron thickness, they tend to become unstable owing to the influence of long-range dispersion forces. In the current study, we investigate how gravitational drainage affects the stability attributes of such thin liquid films. Using scaling arguments, we demonstrate that gravity and dispersion forces can exert their influence simultaneously over a wide range of film thicknesses. In the lubrication limit, we carry out linear stability analysis and nonlinear simulations to understand the evolution of draining thin films. Linear stability indicates the existence of two unstable modes and two cutoff wave numbers, as opposed to a single unstable mode and a unique cutoff wave number observed in stationary films. It is also found that surfactant-laden flowing films are more stable than stationary films with surfactants as well as draining films with clean interfaces. The origin of stabilization is identified as the enhanced surfactant perturbations generated due to drainage. We demonstrate that films exhibiting intermediate levels of surfactant activity and significant drainage exhibit the lowest rates of disturbance growth, leading to extending the time of rupture.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"55 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775936","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}
引用次数: 0
Analytical solutions for long-time steady state Boussinesq gravity currents flowing along a horizontal boundary of finite length 沿有限长度水平边界流动的长时间稳态布辛斯重力流的解析解
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-25 DOI: 10.1103/physrevfluids.9.074803
Safir Haddad, Samuel Vaux, Kevin Varrall, Olivier Vauquelin
This paper presents analytical solutions for a steady turbulent miscible gravity current flowing along a horizontal rigid boundary of finite length into a quiescent uniform environment. These solutions are obtained from the governing equations (mass, momentum, and buoyancy) originally proposed by Ellison and Turner [<span>J. Fluid Mech.</span> <b>6</b>, 423 (1959)] for a buoyant layer of fluid in the Boussinesq approximation. For a constant drag coefficient <math xmlns="http://www.w3.org/1998/Math/MathML"><msub><mi>C</mi><mi>d</mi></msub></math> and the specific entrainment law <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mi>E</mi><mo>∝</mo><msup><mtext>Ri</mtext><mrow><mo>−</mo><mn>1</mn></mrow></msup></mrow><mo>,</mo><mo> </mo><mtext>Ri</mtext></math> being the local Richardson number, we first derived a system of coupled ordinary differential equations describing the longitudinal evolution of the velocity <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>u</mi></math>, the height <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>h</mi></math>, the density deficit <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>η</mi></math>, and the Richardson number <math xmlns="http://www.w3.org/1998/Math/MathML"><mtext>Ri</mtext></math> of the current. For an initially supercritical flow <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mo>(</mo><msub><mtext>Ri</mtext><mn>0</mn></msub><mrow><mspace width="0.16em"></mspace><mo><</mo><mspace width="0.16em"></mspace></mrow><mn>1</mn><mo>)</mo></mrow></math>, explicit relations are found for <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mi>u</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>,</mo><mo> </mo><mrow><mi>h</mi><mo>(</mo><mi>x</mi><mo>)</mo><mo>,</mo></mrow></math> and <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mi>η</mi><mo>(</mo><mi>x</mi><mo>)</mo></mrow></math> solely as a function of the Richardson number <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mtext>Ri</mtext><mo>(</mo><mi>x</mi><mo>)</mo></mrow></math>. The longitudinal evolution of the Richardson number is then theoretically obtained from a universal function <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>F</mi></math> which can be tabulated and, as in the present paper, also plotted. The function <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>F</mi></math> allows us to determine (and only from the knowledge of the boundary conditions at the source) whether the flow remains supercritical over the whole length of the rigid boundary, or might transit towards a subcritical state (<math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mtext>Ri</mtext><mrow><mspace width="0.16em"></mspace><mo>></mo><mspace width="0.16em"></mspace></mrow><mn>1</mn></mrow></math>). In this latter case, the mathematical resolution is modified by including a discontinuity similar to a hydraulic jump. The location and amplitude of this discontinuity are calculated from an additional univers
本文提出了沿有限长度水平刚性边界流向静态均匀环境的稳定湍流混杂重力流的解析解。这些解法来自 Ellison 和 Turner [J. Fluid Mech. 6, 423 (1959)]最初提出的布森斯克近似浮力层流体的控制方程(质量、动量和浮力)。对于恒定的阻力系数 Cd 和特定的夹带定律 E∝Ri-1(Ri 是当地的理查德森数),我们首先导出了一个耦合常微分方程系统,描述了水流的速度 u、高度 h、密度亏损 η 和理查德森数 Ri 的纵向演变。对于初始超临界流(Ri0<1),可以发现 u(x)、h(x) 和 η(x) 完全是理查森数 Ri(x) 的函数。理查德森数的纵向演化可以通过一个通用函数 F 从理论上得到,该函数可以制表,在本文中还可以绘制成图。通过函数 F,我们可以确定(而且只能根据对源头边界条件的了解)流动是在刚性边界的整个长度上保持超临界状态,还是可能向亚临界状态(Ri>1)过渡。在后一种情况下,通过加入一个类似于水力跃迁的不连续性来修改数学分辨率。这种不连续性的位置和振幅是通过附加的通用函数 G 和注入条件计算得出的。该方法最后扩展到为其他经典夹带定律提供解析解。
{"title":"Analytical solutions for long-time steady state Boussinesq gravity currents flowing along a horizontal boundary of finite length","authors":"Safir Haddad, Samuel Vaux, Kevin Varrall, Olivier Vauquelin","doi":"10.1103/physrevfluids.9.074803","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074803","url":null,"abstract":"This paper presents analytical solutions for a steady turbulent miscible gravity current flowing along a horizontal rigid boundary of finite length into a quiescent uniform environment. These solutions are obtained from the governing equations (mass, momentum, and buoyancy) originally proposed by Ellison and Turner [&lt;span&gt;J. Fluid Mech.&lt;/span&gt; &lt;b&gt;6&lt;/b&gt;, 423 (1959)] for a buoyant layer of fluid in the Boussinesq approximation. For a constant drag coefficient &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;msub&gt;&lt;mi&gt;C&lt;/mi&gt;&lt;mi&gt;d&lt;/mi&gt;&lt;/msub&gt;&lt;/math&gt; and the specific entrainment law &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mi&gt;E&lt;/mi&gt;&lt;mo&gt;∝&lt;/mo&gt;&lt;msup&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;mrow&gt;&lt;mo&gt;−&lt;/mo&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/mrow&gt;&lt;mo&gt;,&lt;/mo&gt;&lt;mo&gt; &lt;/mo&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;/math&gt; being the local Richardson number, we first derived a system of coupled ordinary differential equations describing the longitudinal evolution of the velocity &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;u&lt;/mi&gt;&lt;/math&gt;, the height &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;h&lt;/mi&gt;&lt;/math&gt;, the density deficit &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;η&lt;/mi&gt;&lt;/math&gt;, and the Richardson number &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;/math&gt; of the current. For an initially supercritical flow &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;msub&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;mn&gt;0&lt;/mn&gt;&lt;/msub&gt;&lt;mrow&gt;&lt;mspace width=\"0.16em\"&gt;&lt;/mspace&gt;&lt;mo&gt;&lt;&lt;/mo&gt;&lt;mspace width=\"0.16em\"&gt;&lt;/mspace&gt;&lt;/mrow&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/math&gt;, explicit relations are found for &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mi&gt;u&lt;/mi&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mi&gt;x&lt;/mi&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;mo&gt;,&lt;/mo&gt;&lt;mo&gt; &lt;/mo&gt;&lt;mrow&gt;&lt;mi&gt;h&lt;/mi&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mi&gt;x&lt;/mi&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;mo&gt;,&lt;/mo&gt;&lt;/mrow&gt;&lt;/math&gt; and &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mi&gt;η&lt;/mi&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mi&gt;x&lt;/mi&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/math&gt; solely as a function of the Richardson number &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;mo&gt;(&lt;/mo&gt;&lt;mi&gt;x&lt;/mi&gt;&lt;mo&gt;)&lt;/mo&gt;&lt;/mrow&gt;&lt;/math&gt;. The longitudinal evolution of the Richardson number is then theoretically obtained from a universal function &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;F&lt;/mi&gt;&lt;/math&gt; which can be tabulated and, as in the present paper, also plotted. The function &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;F&lt;/mi&gt;&lt;/math&gt; allows us to determine (and only from the knowledge of the boundary conditions at the source) whether the flow remains supercritical over the whole length of the rigid boundary, or might transit towards a subcritical state (&lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mtext&gt;Ri&lt;/mtext&gt;&lt;mrow&gt;&lt;mspace width=\"0.16em\"&gt;&lt;/mspace&gt;&lt;mo&gt;&gt;&lt;/mo&gt;&lt;mspace width=\"0.16em\"&gt;&lt;/mspace&gt;&lt;/mrow&gt;&lt;mn&gt;1&lt;/mn&gt;&lt;/mrow&gt;&lt;/math&gt;). In this latter case, the mathematical resolution is modified by including a discontinuity similar to a hydraulic jump. The location and amplitude of this discontinuity are calculated from an additional univers","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"41 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775939","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}
引用次数: 0
Reactive Rayleigh-Taylor turbulence: Influence of mixing on the growth and displacement of the mixing zone 反应雷利-泰勒湍流:混合对混合区增长和位移的影响
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-25 DOI: 10.1103/physrevfluids.9.074609
Kevin Ley, Olivier Soulard, Jérôme Griffond, Antoine Briard, Serge Simoëns
The purpose of this paper is to investigate the effects of molecular mixing on the evolution of a reactive Rayleigh-Taylor turbulent mixing zone. In this regard, we derive algebraic relations showing that an increase in the mixing level leads to a slowing of the growth of the mixing zone width. We also show the existence of a maximum displacement velocity of the mixing zone center. These predictions are assessed using both direct numerical simulations and large eddy simulations.
本文旨在研究分子混合对反应雷利-泰勒湍流混合区演变的影响。在这方面,我们推导出代数关系,表明混合程度的增加会导致混合区宽度的增长放缓。我们还证明了混合区中心最大位移速度的存在。我们利用直接数值模拟和大涡流模拟对这些预测进行了评估。
{"title":"Reactive Rayleigh-Taylor turbulence: Influence of mixing on the growth and displacement of the mixing zone","authors":"Kevin Ley, Olivier Soulard, Jérôme Griffond, Antoine Briard, Serge Simoëns","doi":"10.1103/physrevfluids.9.074609","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074609","url":null,"abstract":"The purpose of this paper is to investigate the effects of molecular mixing on the evolution of a reactive Rayleigh-Taylor turbulent mixing zone. In this regard, we derive algebraic relations showing that an increase in the mixing level leads to a slowing of the growth of the mixing zone width. We also show the existence of a maximum displacement velocity of the mixing zone center. These predictions are assessed using both direct numerical simulations and large eddy simulations.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"351 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775937","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}
引用次数: 0
Analysis of a turbulent round jet based on direct numerical simulation data at large box and high Reynolds number 基于大方框和高雷诺数下直接数值模拟数据的湍流圆形射流分析
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-25 DOI: 10.1103/physrevfluids.9.074608
Cat Tuong Nguyen, Martin Oberlack
We have conducted a direct numerical simulation of a turbulent round jet at a previously unattained Reynolds number of Re=3500 based on the jet diameter D and jet-inlet bulk velocity Ub in a particularly long box of 75D. To achieve very fast convergence to self-similarity, we used a turbulent pipe flow at the same Reynolds number and length 5D as the upstream inflow boundary condition. This indeed results in a very rapid emergence of self-similarity already at very small axial distances z compared to all turbulent jet data published so far. Not only for the mean velocities and the Reynolds stresses as well as the budgets of the Reynolds stress tensor and the turbulent kinetic energy, a nearly perfect classical scaling based on the normalized radius η=r/z in the range z/D=2565 is shown, but also for the probability density function (PDF) of the axial velocity Uz as well as the associated skewness and kurtosis. All budget terms have been calculated directly, resulting in a marginal error in the balance. An almost completely Gaussian behavior of the PDF for the axial velocity is observed on the jet axis, while a clear deviation with increasingly heavy tails is evident with increasing distance from the axis.
我们在一个 75D 的特长箱中,根据射流直径 D 和射流入口体积速度 Ub,在 Re=3500 的雷诺数条件下对湍流圆形射流进行了直接数值模拟。为了实现自相似性的快速收敛,我们使用了相同雷诺数和长度为 5D 的湍流管道流作为上游流入边界条件。与迄今为止公布的所有湍流射流数据相比,这确实导致在非常小的轴向距离 z 时就已经非常快速地出现了自相似性。在 z/D=25-65 范围内,不仅平均速度和雷诺应力以及雷诺应力张量和湍流动能的预算显示了基于归一化半径 η=r/z 的近乎完美的经典缩放,而且轴向速度 Uz 的概率密度函数(PDF)以及相关的偏度和峰度也显示了近乎完美的经典缩放。所有预算项都是直接计算得出的,因此平衡中存在边际误差。在射流轴上可以观察到轴向速度的概率密度函数几乎完全呈高斯分布,而随着与轴线距离的增加,明显出现了尾部越来越大的偏差。
{"title":"Analysis of a turbulent round jet based on direct numerical simulation data at large box and high Reynolds number","authors":"Cat Tuong Nguyen, Martin Oberlack","doi":"10.1103/physrevfluids.9.074608","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.074608","url":null,"abstract":"We have conducted a direct numerical simulation of a turbulent round jet at a previously unattained Reynolds number of <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mtext>Re</mtext><mo>=</mo><mn>3500</mn></mrow></math> based on the jet diameter <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>D</mi></math> and jet-inlet bulk velocity <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi>U</mi><mtext>b</mtext></msub></math> in a particularly long box of <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mn>75</mn><mi>D</mi></mrow></math>. To achieve very fast convergence to self-similarity, we used a turbulent pipe flow at the same Reynolds number and length <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mn>5</mn><mi>D</mi></mrow></math> as the upstream inflow boundary condition. This indeed results in a very rapid emergence of self-similarity already at very small axial distances <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mi>z</mi></math> compared to all turbulent jet data published so far. Not only for the mean velocities and the Reynolds stresses as well as the budgets of the Reynolds stress tensor and the turbulent kinetic energy, a nearly perfect classical scaling based on the normalized radius <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi>η</mi><mo>=</mo><mi>r</mi><mo>/</mo><mi>z</mi></mrow></math> in the range <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mi>z</mi><mo>/</mo><mi>D</mi><mo>=</mo><mn>25</mn><mo>−</mo><mn>65</mn></mrow></math> is shown, but also for the probability density function (PDF) of the axial velocity <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><msub><mi>U</mi><mi>z</mi></msub></math> as well as the associated skewness and kurtosis. All budget terms have been calculated directly, resulting in a marginal error in the balance. An almost completely Gaussian behavior of the PDF for the axial velocity is observed on the jet axis, while a clear deviation with increasingly heavy tails is evident with increasing distance from the axis.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"12 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-25","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141785443","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}
引用次数: 0
Heat transport in three-layer turbulent thermal convection 三层湍流热对流中的热传输
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-24 DOI: 10.1103/physrevfluids.9.073501
Xiao-Zheng Zhao, Can Qiu, Sheng-Qi Zhou, Yi-Zhen Li, Heng-Dong Xi, Ke-Qing Xia
We report an experimental study of heat transport in a three-layer turbulent Rayleigh-Bénard convection. The experiments were conducted in a cylindrical cell (with diameter <math xmlns="http://www.w3.org/1998/Math/MathML"><mi>D</mi></math>) filled with a FC77 layer with height <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mi>H</mi><mo>=</mo><mi>D</mi></mrow></math>. A very thin layer of water and a very thin layer of mercury were introduced to the top and bottom of the FC77 layer to provide slippery boundary conditions. We performed high spatial resolution temperature measurements across the water-FC77 and FC77-mercury interfaces, determined the temperatures at the two interfaces, the Rayleigh number (Ra) and the Nusselt number (Nu) across the FC77 layer. The experiments were conducted in the Ra range of <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mn>2.81</mn><mo>×</mo><msup><mn>10</mn><mn>9</mn></msup></mrow></math> to <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mn>1.24</mn><mo>×</mo><msup><mn>10</mn><mn>11</mn></msup></mrow></math> for the FC77 layer. It is found that not only the amplitude but also the scaling exponent (with Ra) of Nu is greatly enhanced in this three-layer system compared to the canonical single-layer system, especially in the high Ra range. In particular, <math xmlns="http://www.w3.org/1998/Math/MathML"><mtext>Nu</mtext></math> first scales as <math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mtext>Ra</mtext><mrow><mn>0.31</mn></mrow></msup></math> and then <math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mtext>Ra</mtext><mrow><mn>0.38</mn></mrow></msup></math> when Ra exceeds a transitional Rayleigh number <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><msub><mtext>Ra</mtext><mi>t</mi></msub><mo>=</mo><mn>2.52</mn><mo>×</mo><msup><mn>10</mn><mn>10</mn></msup></mrow></math>, whereas in the canonical single-layer FC77 case, <math xmlns="http://www.w3.org/1998/Math/MathML"><mtext>Nu</mtext></math> is found to scale as <math xmlns="http://www.w3.org/1998/Math/MathML"><msup><mtext>Ra</mtext><mrow><mn>0.26</mn></mrow></msup></math>. Temperature measurements show that the boundary condition above and below the FC77 layer is asymmetric especially when <math xmlns="http://www.w3.org/1998/Math/MathML"><mrow><mtext>Ra</mtext><mo>></mo><msub><mtext>Ra</mtext><mi>t</mi></msub></mrow></math>: the temperature drop across the top half (in contact with the water layer) of the FC77 layer is smaller than that across the bottom half (in contact with the mercury layer), and the top thermal boundary layer (TBL) becomes thinner and follows a steeper scaling with <math xmlns="http://www.w3.org/1998/Math/MathML"><mtext>Ra</mtext></math> compared to the bottom TBL. We consider a hypothetical experiment where the top and the bottom boundary conditions are symmetric, denoted as a “water-FC77-water” three-layer system, in which the temperature drop across the bottom boundary layer <math xmlns="h
我们报告了对三层湍流雷利-贝纳德对流中热量传输的实验研究。实验在一个充满 FC77 层(高度为 H=D)的圆柱形单元(直径为 D)中进行。在 FC77 层的顶部和底部分别引入了极薄的水层和极薄的汞层,以提供滑动边界条件。我们在水-FC77 和 FC77-汞界面上进行了高空间分辨率温度测量,确定了两个界面上的温度、FC77 层上的雷利数(Ra)和努塞尔特数(Nu)。实验在 FC77 层的 Ra 范围 2.81×109 至 1.24×1011 之间进行。实验发现,与典型的单层系统相比,在这种三层系统中,Nu 的振幅和缩放指数(随 Ra 变化)都大大增强,尤其是在高 Ra 范围内。特别是,当 Ra 超过过渡瑞利数 Rat=2.52×1010 时,Nu 首先按 Ra0.31 的比例缩放,然后按 Ra0.38 的比例缩放,而在典型的单层 FC77 情况下,Nu 按 Ra0.26 的比例缩放。温度测量结果表明,FC77 层上下的边界条件是不对称的,尤其是当 Ra>Rat 时:FC77 层上半层(与水层接触)的温降小于下半层(与水银层接触)的温降,而且顶部热边界层(TBL)变得更薄,与底部 TBL 相比,随 Ra 变化的比例更陡峭。我们考虑了一个顶部和底部边界条件对称的假设实验,称为 "水-FC77-水 "三层系统,其中底部边界层的温降 ΔTb 与顶部边界层的温降 ΔTt 相同。我们发现,在这个水-FC77-水三层体系中,随着 Ra 的增加,Nu 与 Ra 的比例关系从 Nu∼Ra0.31 过渡到 Nu∼Ra0.46,过渡期间的 Ra 与之前确定的 Rat 相同。仔细观察水层、FC77 层和汞层的 Ra 演变,可以发现 Nu vs Ra 缩放的过渡是由于薄水层从传导状态过渡到对流状态,而汞层仍然处于传导状态。
{"title":"Heat transport in three-layer turbulent thermal convection","authors":"Xiao-Zheng Zhao, Can Qiu, Sheng-Qi Zhou, Yi-Zhen Li, Heng-Dong Xi, Ke-Qing Xia","doi":"10.1103/physrevfluids.9.073501","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.073501","url":null,"abstract":"We report an experimental study of heat transport in a three-layer turbulent Rayleigh-Bénard convection. The experiments were conducted in a cylindrical cell (with diameter &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mi&gt;D&lt;/mi&gt;&lt;/math&gt;) filled with a FC77 layer with height &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mi&gt;H&lt;/mi&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;mi&gt;D&lt;/mi&gt;&lt;/mrow&gt;&lt;/math&gt;. A very thin layer of water and a very thin layer of mercury were introduced to the top and bottom of the FC77 layer to provide slippery boundary conditions. We performed high spatial resolution temperature measurements across the water-FC77 and FC77-mercury interfaces, determined the temperatures at the two interfaces, the Rayleigh number (Ra) and the Nusselt number (Nu) across the FC77 layer. The experiments were conducted in the Ra range of &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mn&gt;2.81&lt;/mn&gt;&lt;mo&gt;×&lt;/mo&gt;&lt;msup&gt;&lt;mn&gt;10&lt;/mn&gt;&lt;mn&gt;9&lt;/mn&gt;&lt;/msup&gt;&lt;/mrow&gt;&lt;/math&gt; to &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mn&gt;1.24&lt;/mn&gt;&lt;mo&gt;×&lt;/mo&gt;&lt;msup&gt;&lt;mn&gt;10&lt;/mn&gt;&lt;mn&gt;11&lt;/mn&gt;&lt;/msup&gt;&lt;/mrow&gt;&lt;/math&gt; for the FC77 layer. It is found that not only the amplitude but also the scaling exponent (with Ra) of Nu is greatly enhanced in this three-layer system compared to the canonical single-layer system, especially in the high Ra range. In particular, &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mtext&gt;Nu&lt;/mtext&gt;&lt;/math&gt; first scales as &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;msup&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mrow&gt;&lt;mn&gt;0.31&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/math&gt; and then &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;msup&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mrow&gt;&lt;mn&gt;0.38&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/math&gt; when Ra exceeds a transitional Rayleigh number &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;msub&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mi&gt;t&lt;/mi&gt;&lt;/msub&gt;&lt;mo&gt;=&lt;/mo&gt;&lt;mn&gt;2.52&lt;/mn&gt;&lt;mo&gt;×&lt;/mo&gt;&lt;msup&gt;&lt;mn&gt;10&lt;/mn&gt;&lt;mn&gt;10&lt;/mn&gt;&lt;/msup&gt;&lt;/mrow&gt;&lt;/math&gt;, whereas in the canonical single-layer FC77 case, &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mtext&gt;Nu&lt;/mtext&gt;&lt;/math&gt; is found to scale as &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;msup&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mrow&gt;&lt;mn&gt;0.26&lt;/mn&gt;&lt;/mrow&gt;&lt;/msup&gt;&lt;/math&gt;. Temperature measurements show that the boundary condition above and below the FC77 layer is asymmetric especially when &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mrow&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mo&gt;&gt;&lt;/mo&gt;&lt;msub&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;mi&gt;t&lt;/mi&gt;&lt;/msub&gt;&lt;/mrow&gt;&lt;/math&gt;: the temperature drop across the top half (in contact with the water layer) of the FC77 layer is smaller than that across the bottom half (in contact with the mercury layer), and the top thermal boundary layer (TBL) becomes thinner and follows a steeper scaling with &lt;math xmlns=\"http://www.w3.org/1998/Math/MathML\"&gt;&lt;mtext&gt;Ra&lt;/mtext&gt;&lt;/math&gt; compared to the bottom TBL. We consider a hypothetical experiment where the top and the bottom boundary conditions are symmetric, denoted as a “water-FC77-water” three-layer system, in which the temperature drop across the bottom boundary layer &lt;math xmlns=\"h","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"27 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141775938","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}
引用次数: 0
Modeling a spheroidal squirmer through a complex fluid 复杂流体中的球形唧筒建模
IF 2.7 3区 物理与天体物理 Q2 PHYSICS, FLUIDS & PLASMAS Pub Date : 2024-07-23 DOI: 10.1103/physrevfluids.9.073303
Zhenyu Ouyang, Chen Liu, Zhaowu Lin, Jianzhong Lin
We simulate a spheroidal swimmer through a complex fluid, modeled by the Giesekus constitutive equation incorporating fluid inertia. We develop a spheroidal swimmer model and exert it in a direct-forcing fictitious domain method framework. This model extends the conventional spherical “squirmer,” representing a microswimmer generating self-propulsion through tangential surface waves at its boundaries. We vary the swimmer's aspect ratio (AR) and Weissenberg number (Wi; the ratio of fluid elastic force to viscous force), respectively, in the range of 1.5AR8 and 0.5Wi10. Our results show that, an inertial spheroidal puller with a small |β| (a swimming intensity parameter) swims faster than the counterpart subjected to the Stokes flow regime—a departure from the observed pattern in spherical pullers. Within the Giesekus fluid medium, an augmented mobility factor α correlates with an increased squirmer velocity, while a larger AR contributes significantly to the speed enhancement of a neutral squirmer in the presence of fluid inertia. Meanwhile, we explore the squirmer's energy expenditure and hydrodynamic efficiency, finding that a slenderer, inertial squirmer with a vigorous swimming intensity expends more energy, contrasting with the reduced energy expenditure associated with a smaller intensity. Notably, a larger AR positively correlates with squirmer efficiency, displaying an advantageous relationship with swimming speed.
我们用包含流体惯性的 Giesekus 构成方程模拟了球形游泳者在复杂流体中的运动。我们建立了一个球形游泳者模型,并在直接强迫虚构域法框架中对其进行了应用。该模型扩展了传统的球形 "squirmer",代表了通过边界切向表面波产生自推进力的微型游泳者。我们分别在 1.5≤AR≤8 和 0.5≤Wi≤10 的范围内改变游泳者的长宽比(AR)和魏森伯格数(Wi;流体弹性力与粘性力之比)。我们的研究结果表明,具有较小 |β|(游泳强度参数)的惯性球形拉力器的游泳速度快于受斯托克斯流机制影响的球形拉力器,这与在球形拉力器中观察到的模式不同。在 Giesekus 流体介质中,增大的流动因子 α 与松鼠速度的增加相关,而在存在流体惯性的情况下,较大的 AR 会显著提高中性松鼠的速度。同时,我们还探讨了松鼠的能量消耗和流体动力学效率,发现在游泳强度较大的情况下,纤细的惯性松鼠会消耗更多的能量,而在游泳强度较小的情况下,能量消耗则会减少。值得注意的是,较大的AR与松鼠的效率呈正相关,与游泳速度呈优势关系。
{"title":"Modeling a spheroidal squirmer through a complex fluid","authors":"Zhenyu Ouyang, Chen Liu, Zhaowu Lin, Jianzhong Lin","doi":"10.1103/physrevfluids.9.073303","DOIUrl":"https://doi.org/10.1103/physrevfluids.9.073303","url":null,"abstract":"We simulate a spheroidal swimmer through a complex fluid, modeled by the Giesekus constitutive equation incorporating fluid inertia. We develop a spheroidal swimmer model and exert it in a direct-forcing fictitious domain method framework. This model extends the conventional spherical “squirmer,” representing a microswimmer generating self-propulsion through tangential surface waves at its boundaries. We vary the swimmer's aspect ratio (AR) and Weissenberg number (Wi; the ratio of fluid elastic force to viscous force), respectively, in the range of <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mn>1.5</mn><mo>≤</mo><mi>AR</mi><mo>≤</mo><mn>8</mn></mrow></math> and <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mn>0.5</mn><mo>≤</mo><mi>Wi</mi><mo>≤</mo><mn>10</mn></mrow></math>. Our results show that, an inertial spheroidal puller with a small <math xmlns=\"http://www.w3.org/1998/Math/MathML\"><mrow><mo>|</mo><mi>β</mi><mo>|</mo></mrow></math> (a swimming intensity parameter) swims faster than the counterpart subjected to the Stokes flow regime—a departure from the observed pattern in spherical pullers. Within the Giesekus fluid medium, an augmented mobility factor <i>α</i> correlates with an increased squirmer velocity, while a larger AR contributes significantly to the speed enhancement of a neutral squirmer in the presence of fluid inertia. Meanwhile, we explore the squirmer's energy expenditure and hydrodynamic efficiency, finding that a slenderer, inertial squirmer with a vigorous swimming intensity expends more energy, contrasting with the reduced energy expenditure associated with a smaller intensity. Notably, a larger AR positively correlates with squirmer efficiency, displaying an advantageous relationship with swimming speed.","PeriodicalId":20160,"journal":{"name":"Physical Review Fluids","volume":"38 1","pages":""},"PeriodicalIF":2.7,"publicationDate":"2024-07-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141754059","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}
引用次数: 0
期刊
Physical Review Fluids
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1