{"title":"Influence of surface slip on hydrodynamics and flow field around a two-dimensional hydrofoil at a moderate Reynolds number","authors":"Manfu Zhu, Weixi Huang, Liran Ma, Jianbin Luo","doi":"10.1063/5.0203389","DOIUrl":"https://doi.org/10.1063/5.0203389","url":null,"abstract":"In the present study, the effects of surface slip on the hydrodynamics and flow around a two-dimensional National Advisory Committee for Aeronautics 0012 hydrofoil are systematically investigated by numerical methods. The objective is to fully understand the effects of surface slip on the streamlined body. Three slip positions (both surfaces, the upper surface, the lower surface) and eight slip lengths (in a wide range from 1 to 500 μm) under 0°–10° angles of attack are fully investigated at a moderate Reynolds number of 1.0 × 106. Surface slip has been found to increase lift and reduce drag by postponing the flow transition, laminar separation bubble, and flow separation on the hydrofoil surface under both surfaces and the upper surface slip conditions. Slip has also been found to induce upshift of the mean velocity profile, decrease the displacement thickness, and mitigate the turbulent kinetic energy in the flow field. However, counterintuitive phenomenon occurs under the lower surface slip condition, where the total drag of the hydrofoil is increased compared to that under the no slip condition. Total drag increase is found mainly due to the increase in the pressure drag under small slip lengths and relatively large angles of attack. Flow maps demonstrating the complex interaction between different surface slip conditions and the flow field are further presented. The results suggest that surface slip can not only reduce drag, but also increase the drag of the streamlined body, which shall provide valuable insights for practical applications of slippery materials.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"4 9","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141045462","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}
Pratik Mahyawansi, Sumit R. Zanje, Abbas Sharifi, Dwayne McDaniel, Arturo S. Leon
{"title":"Experimental investigation of storm sewer geyser using a large-scale setup","authors":"Pratik Mahyawansi, Sumit R. Zanje, Abbas Sharifi, Dwayne McDaniel, Arturo S. Leon","doi":"10.1063/5.0199012","DOIUrl":"https://doi.org/10.1063/5.0199012","url":null,"abstract":"The storm sewer geyser is a process where an air–water mixture violently erupts from a manhole. Despite the low hydrostatic pressure, violent eruptions can achieve a height of tens of meters above the ground. This current study experimentally investigates large-scale violent geysers using a large air pocket inserted from a pressurized air tank. The total length of the pipe system is approximately 88 m with a 0.1572 m diameter pipe. This large-scale experiment facilitates the investigation of spontaneous geyser eruptions. This study identifies the role of air–water volume ratio and coefficient of pressure (ratio of absolute initial static pressure to initial dynamic pressure) on the geyser intensity using eruption images and pressure plots. A total of 116 cases are tested, in which the volume ratio is parametrically increased from 0 to 1.1 under various operating conditions. A geyser score is defined to quantify the geyser eruption nature based on visual observations. The key findings are as follows: first, a sharp transition in geyser intensity is observed at the critical volume ratio of 0.5, and pre-transition and post-transition intensity exhibit a linear relationship with the volume ratio; and second, the critical volume ratio linearly varies with the coefficient of pressure.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"1 8","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141047935","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}
{"title":"Controlling Couette flow by alternating axial mass flux","authors":"Sebastian A. Altmeyer","doi":"10.1063/5.0197106","DOIUrl":"https://doi.org/10.1063/5.0197106","url":null,"abstract":"This paper presents numerical simulations of the Taylor vortex flow under the influence of an externally applied alternating axial mass flux (through-flow) in a Taylor–Couette system with axial periodic boundary conditions. Such an axially modulating flow can lead to a significant variation in the onset of primary instabilities. Depending on the system parameters, the effect can be both stabilizing and destabilizing, i.e., shifting the bifurcation threshold to larger or smaller control parameters, respectively. It is found that the system response around the primary instability is sensitive to and critically influenced by an alternating mass flux, particularly the modulation frequency. We show that such an alternating axial flow represents an easily and, more importantly, precisely controllable key parameter to change the nonlinear system response from subcritical to supercritical behavior and vice versa. Furthermore, we observe different parameter regimes with regular and irregular intermittent flow dynamics.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"46 28","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141141916","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}
Chenxi Qin, L. Duan, Duoyin Wang, Bingchuang Duan, Wei Liu
{"title":"A local scour model for single pile on silty seabed considering soil cohesion (SedCohFOAM): Model and validation","authors":"Chenxi Qin, L. Duan, Duoyin Wang, Bingchuang Duan, Wei Liu","doi":"10.1063/5.0207743","DOIUrl":"https://doi.org/10.1063/5.0207743","url":null,"abstract":"In this study, the sediment transport two-phase flow model named SedFOAM is expanded to include soil cohesion, creating a new model named SedCohFOAM within OpenFOAM. The local scouring flume experiment involving a pile on silty seabed and sandy seabed is conducted in a curved flume. Due to the influence of cohesion, the scouring depth at different locations on sandy seabed is 15%–18% greater than that on silty seabed. Observations from this experiment informed the analysis of force balance, wherein the agglomerated silt particles are modeled as large singular entities and the cohesive force is treated as a downward influence that keeps the aggregated particles stationary. Meanwhile, the experimental outcomes are utilized to validate the accuracy of the SedCohFOAM model. The numerical findings demonstrated that SedCohFOAM can simulate the flow field distribution around the pile, variations in seabed shear stress, and alterations in seabed surface morphology. Compared with the SedFOAM model, the SedCohFOAM model has a significantly reduced simulation error in simulating scour on silty seabed. When comparing the cross-sectional profiles of the scour holes derived from the flume experiments with those simulated by SedCohFOAM, it was observed that the ultimate-equilibrium scour depth predicted by the model is consistently lower, but the scour radius in the numerical simulations is larger. The deviation from the experimental results is nearly within 8%, while when the flow velocity is high, the simulation error of the simulated scouring depth behind the pile and the scouring radius in front of pile is amplified.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"45 17","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141142360","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}
{"title":"A wall-modeled immersed boundary/large eddy simulation method and its application to simulating heart valve flows","authors":"Jingyang Wang, T. Pu, Chunhua Zhou","doi":"10.1063/5.0198734","DOIUrl":"https://doi.org/10.1063/5.0198734","url":null,"abstract":"In this work, a wall-modeled immersed boundary (IB)/large eddy simulation (LES) method is extended to the simulation of moving-boundary flows. The used non-equilibrium algebraic wall model is based on an assumed velocity profile, the coefficients of which are determined from physical constraints provided by the full turbulent-boundary-layer equations. To implement the wall model in an IB method named the local domain-free discretization (DFD) method, a local coordinate system fixed on the moving body is introduced. Thus, wall modeling is transformed into a local two-dimensional problem and the complexity of implementation of the wall model is reduced. In the present LES-DFD method, the tangential velocity at an exterior dependent node is determined via wall shear stress prescribed by the wall model. To reduce computational cost for simulating an internal flow with moving boundaries, the stationary boundaries are handled by the body-fitted-grid method and the moving boundaries by the local DFD method. There is no need of an auxiliary grid for solving the non-equilibrium algebraic wall model. Therefore, the inbuilt advantage of an IB method can be retained when simulating moving-boundary problems, and the economy of equilibrium wall models can also be preserved. The present method is applied to simulating the pulsatile flows through a bileaflet mechanical heart valve implanted in a model aorta. The predicted results show an acceptable agreement with the referenced experimental measurements or numerical results at much higher resolution and the applicability of the non-equilibrium wall model to LES of complex moving-boundary flows is verified.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"23 6","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141030910","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}
{"title":"Field inversion machine learning augmented turbulence modeling for time-accurate unsteady flow","authors":"Lean Fang, Ping He","doi":"10.1063/5.0207704","DOIUrl":"https://doi.org/10.1063/5.0207704","url":null,"abstract":"Field inversion machine learning (FIML) has the advantages of model consistency and low data dependency and has been used to augment imperfect turbulence models. However, the solver-intrusive field inversion has a high entry bar, and existing FIML studies focused on improving only steady-state or time-averaged periodic flow predictions. To break this limit, this paper develops an open-source FIML framework for time-accurate unsteady flow, where both spatial and temporal variations of flow are of interest. We augment a Reynolds-Averaged Navier–Stokes (RANS) turbulence model's production term with a scalar field. We then integrate a neural network (NN) model into the flow solver to compute the above augmentation scalar field based on local flow features at each time step. Finally, we optimize the weights and biases of the built-in NN model to minimize the regulated spatial-temporal prediction error between the augmented flow solver and reference data. We consider the spatial-temporal evolution of unsteady flow over a 45° ramp and use only the surface pressure as the training data. The unsteady-FIML-trained model accurately predicts the spatial-temporal variations of unsteady flow fields. In addition, the trained model exhibits reasonably good prediction accuracy for various ramp angles, Reynolds numbers, and flow variables (e.g., velocity fields) that are not used in training, highlighting its generalizability. The FIML capability has been integrated into our open-source framework DAFoam. It has the potential to train more accurate RANS turbulence models for other unsteady flow phenomena, such as wind gust response, bubbly flow, and particle dispersion in the atmosphere.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"90 5","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141032359","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}
{"title":"Numerical simulations of bio-inspired approaches to enhance underwater swimming efficiency","authors":"Rahul Kumar, S. Padhee, D. Samanta","doi":"10.1063/5.0201926","DOIUrl":"https://doi.org/10.1063/5.0201926","url":null,"abstract":"The present study discusses the numerical simulation results of swimming similar to manta rays. The complex three-dimensional kinematics of manta rays were implemented to unravel the intricacies of its propulsion mechanisms by using the discrete vortex method (DVM). The DVM replaces the requirement for a structured grid across the computational domain with a collection of vortex elements. This method simplifies grid generation, especially for intricate geometries, resulting in time and effort savings in meshing complex shapes. By modeling the pectoral fins with discrete panels and utilizing vortex rings to represent circulation and wake, the study accurately computes the pressure distribution, circulation distribution, lift coefficient, and thrust coefficient of the manta ray. This study focuses on the modulation of aerodynamic performance by altering the span length and the length change ratio during the downstroke and upstroke motion (SV). The manta ray's three-dimensional vortex configurations comprise a combination of vortex rings, vortex contrails, and horseshoe vortices. Analysis of the three-dimensional vortex structure indicates the presence of multiple vortex rings and horseshoe vortex rings at higher SV values, while adequate formation of horseshoe vortices is not observed at lower SV values. In terms of propulsive performance, both lift and thrust increase with SV, while the propulsive efficiency demonstrates its peak at SV = 1.75. The analysis reveals that at higher SV values, the net thrust generated primarily originates from the tip of the fins. Moreover, the study illustrates a significant enhancement in propulsive efficiency, particularly in association with optimal Strouhal numbers ranging between 0.3 and 0.4. The key findings of this study may be used in efficient design of agile autonomous underwater vehicles for marine exploration and surveillance applications.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"6 4","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141034626","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}
Shaofeng Du, Yue Xiao, Qiao Li, Shaowei Wang, Moli Zhao
{"title":"Stability of viscoelastic film on a slippery inclined plane","authors":"Shaofeng Du, Yue Xiao, Qiao Li, Shaowei Wang, Moli Zhao","doi":"10.1063/5.0210250","DOIUrl":"https://doi.org/10.1063/5.0210250","url":null,"abstract":"The linear and weakly nonlinear stability of viscoelastic film flowing down a slippery inclined plane is investigated analytically. Under the assumption of the long wave approximation, the first-order Benny equation of Oldroyd-B fluid thin film with slip condition is obtained. Through the normal mode analysis, the neutral stability curve and the temporal growth rates are calculated to explore the linear stability of the film. Linear results show that the critical Reynolds number decreases with the increase in slip length and viscoelastic parameter and that the liquid film may exhibit pure elastic instability. For the nonlinear stability analysis, both hydrodynamic instability and elastic instability are discussed. The primary bifurcations in the phase plane are identified by calculating the Landau coefficient, i.e., the unconditional stable region, the supercritical region, the subcritical region, and the explosive region. The dependence of primary bifurcation regions upon the slip length and Deborah number are studied, and the results indicate that the slip boundary and viscoelasticity destabilizes the flow. According to the Ginzburg–Landau equation, the threshold amplitude of the nonlinear equilibrium solution is analyzed as well.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":" 96","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141132150","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}
{"title":"Reference solutions for compressible single-phase flows in heated and cooled ducts","authors":"S. Schropff, F. Petitpas, E. Daniel","doi":"10.1063/5.0209500","DOIUrl":"https://doi.org/10.1063/5.0209500","url":null,"abstract":"Analytical/quasi-analytical solutions are proposed for a steady, compressible, single-phase flow in a rectilinear duct subjected to heating followed by cooling. The flow is driven by the pressure ratio between an upstream tank and a variable outlet pressure. The article proposes a methodology to determine the full flow behavior, as a function of pressure ratio and heat-flux distribution. Following an analogy done with the study of compressible flows in nozzles, a behavioral classification of non-adiabatic compressible flows is proposed through the definition of critical pressure ratios. It is demonstrated that a critical pressure ratio distinguishes subsonic and supersonic outlet regimes and that there cannot be a steady shock wave in such configuration. The behavior of this critical pressure ratio is studied for limit cases of heat flux, delineating physical boundaries. An abacus is also proposed for a given couple of heating and cooling powers, as both values are needed to characterize the flow. Results are studied for parameters such as pressure ratio and outlet heat power. A short validation of a numerical simulation tool is provided, yielding excellent results and very small relative errors.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"108 11","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141135484","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}
Binzhen Zhou, Yi Xiao, Kanglixi Ding, Lei Wang, Yifeng Yang, P. Jin
{"title":"Optimal strategy of the asymmetric wave energy converter survival in extreme waves","authors":"Binzhen Zhou, Yi Xiao, Kanglixi Ding, Lei Wang, Yifeng Yang, P. Jin","doi":"10.1063/5.0208825","DOIUrl":"https://doi.org/10.1063/5.0208825","url":null,"abstract":"Enhancing the survival performance of wave energy converters (WECs) in extreme wave conditions is crucial, and reducing wave loads is a key aspect of this. Placing the device underwater has been recognized as a beneficial strategy, yet the determination of the optimal submerged depth and the effects of varying wave conditions remain ambiguous. To address this, the study numerically analyzes the total forces in both horizontal and vertical directions, along with their harmonic components, across different wave configurations. A computational fluid dynamics method is employed to investigate a triangular-baffle bottom-shaped oscillating floater, which is known for its high energy conversion efficiency. The findings indicate that submerging the device to a depth equivalent to half the actual focused amplitude (1/2Ab) is the most effective strategy in the given sea state, offering superior wave force reduction vertically and robust performance horizontally. The analysis of harmonics reveals the significant contribution of high-order components to the total wave forces. Additionally, the study examines the impact of focused wave amplitudes and peak frequencies, showing that although force reductions are lessened in more extreme conditions, the optimal submerged depth of 1/2Ab still yields near 30% reduction in total vertical force and 22% in total horizontal force. This research provides theoretical insight that can guide the enhancement of WECs' survival capabilities in practical engineering applications.","PeriodicalId":509470,"journal":{"name":"Physics of Fluids","volume":"9 4p2","pages":""},"PeriodicalIF":0.0,"publicationDate":"2024-05-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141141516","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}