Computers & FluidsPub Date : 2025-04-08DOI: 10.1016/j.compfluid.2025.106638
G. Fourtakas , R. Vacondio , B.D. Rogers
{"title":"Divergence cleaning for weakly compressible smoothed particle hydrodynamics","authors":"G. Fourtakas , R. Vacondio , B.D. Rogers","doi":"10.1016/j.compfluid.2025.106638","DOIUrl":"10.1016/j.compfluid.2025.106638","url":null,"abstract":"<div><div>This paper presents a divergence cleaning formulation for the velocity in the weakly compressible smoothed particle hydrodynamics (SPH) scheme. The proposed hyperbolic/parabolic divergence cleaning, ensures that the velocity divergence, <em>div</em>(<strong>u</strong>), is minimised throughout the simulation. The divergence equation is coupled with the momentum conservation equation through a scalar field <em>ψ</em>. A parabolic term is added to the time-evolving divergence equation, resulting in a hyperbolic/parabolic form, dissipating acoustic waves with a speed of sound proportional to the local Mach number in order to maximise dissipation of the velocity divergence, preventing unwanted diffusion of the pressure field. The <em>div</em>(<strong>u</strong>)-SPH algorithm is implemented in the open-source weakly compressible SPH solver DualSPHysics. The new formulation is validated against a range of challenging 2-D test cases including the Taylor-Green vortices, patch impact test, jet impinging on a surface, and wave impact in a sloshing tank. The results show that the new formulation reduces the divergence in the velocity field by at least one order of magnitude which prevents spurious numerical noise and the formation of unphysical voids. The temporal evolution of the impact pressures shows that the <em>div</em>(<strong>u</strong>)-SPH formulation virtually eliminates unwanted acoustic pressure oscillations. Investigation of particle resolution confirms that the new <em>div</em>(<strong>u</strong>)-SPH formulation does not reduce the spatial convergence rate.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"295 ","pages":"Article 106638"},"PeriodicalIF":2.5,"publicationDate":"2025-04-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143870853","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}
Computers & FluidsPub Date : 2025-04-05DOI: 10.1016/j.compfluid.2025.106620
Juan Sheng, Feng He, Yitao Wang, Pengfei Hao, Xiwen Zhang, Xiangru Li
{"title":"Large Eddy simulation of a supersonic starting impinging jet","authors":"Juan Sheng, Feng He, Yitao Wang, Pengfei Hao, Xiwen Zhang, Xiangru Li","doi":"10.1016/j.compfluid.2025.106620","DOIUrl":"10.1016/j.compfluid.2025.106620","url":null,"abstract":"<div><div>Supersonic impinging jets have wide applications in many directions, such as rocket exhaust, supersonic combustor mixing, and accidental leakage of pressurized fluids. The steady processes of impinging jets have been extensively studied, but there is little research on the starting impinging jets. Considering the important applications of starting jet in rocket launches, car's airbags, and pulse jets, the under-expanded round jet impinging on a flat plate normally with a nozzle pressure ratio of 2.7 has been investigated using large eddy simulations. The distance between the impinging plate and the nozzle outlet is set to 5 times the nozzle diameter. The simulation results agree well with the experimental data. The variations of flow structures and dominant frequencies of the acoustic resonance at the starting process are presented. The radiation directions of the dominant frequencies exhibit notable variations at different stages of the starting process, and the explanation for this distinction is explored by the two-dimensional correlation analysis. The wavenumber spectra and dispersion relations are employed, showing that at different stages, the upstream-propagating guided jet modes are excited by the interaction between the Kelvin–Helmholtz wavepacket and shock cells of distinct wavenumbers. Finally, based on the instantaneous pressure distribution, the dynamic evolution characteristics of the flow structures at different stages are analyzed.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"294 ","pages":"Article 106620"},"PeriodicalIF":2.5,"publicationDate":"2025-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143817814","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}
Computers & FluidsPub Date : 2025-04-05DOI: 10.1016/j.compfluid.2025.106630
Mohammadreza Bagheri
{"title":"Modeling two-phase equilibrium using smoothed particle hydrodynamics","authors":"Mohammadreza Bagheri","doi":"10.1016/j.compfluid.2025.106630","DOIUrl":"10.1016/j.compfluid.2025.106630","url":null,"abstract":"<div><div>Smoothed Particle Hydrodynamics (SPH) is an emerging particle-based methodology that can also be used for modeling two-phase flows, currently in the early stages of development. This Lagrangian, mesh-free approach utilizes macro-scale formulations at the meso-scale, achieving computational performance comparable to micro-scale methods. This integration allows for efficient computations at larger scales than micro while facilitating detailed analysis at smaller scales than macro. This paper focuses on the study of two-phase equilibrium and droplet formation, employing the Equation of State (EOS) alongside the careful selection of an appropriate smoothing length. The majority of existing SPH literature utilizes the van der Waals (vdW) EOS for two-phase simulations. While the vdW EOS has provided foundational insights, newer models have been developed to accommodate a broader range of fluids. In this study, the Peng-Robinson EOS is employed, which separates the EOS into attractive and repulsive components, thereby enhancing modeling capabilities. This work critically examines the limitations of SPH in simulating two-phase equilibrium, deriving the smoothing length for attractive forces based on surface tension. Furthermore, it contends that employing an updated smoothing length does not accurately reflect physical realities. To the best of the author's knowledge, this research is among the few that directly integrates the Peng-Robinson Equation of State (PR EOS) and a viscosity equation of state within the SPH framework for the simulation of two-phase equilibrium.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"294 ","pages":"Article 106630"},"PeriodicalIF":2.5,"publicationDate":"2025-04-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143824335","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}
Computers & FluidsPub Date : 2025-04-03DOI: 10.1016/j.compfluid.2025.106621
M. Breuer, K. Boulbrachene, G. De Nayer
{"title":"Combination of two FSI methods and their validation based on artificial wind gusts impacting a flexible T-structure","authors":"M. Breuer, K. Boulbrachene, G. De Nayer","doi":"10.1016/j.compfluid.2025.106621","DOIUrl":"10.1016/j.compfluid.2025.106621","url":null,"abstract":"<div><div>The study focuses on the combination of two numerical approaches that are typically not used together in this manner. The first is a well-established partitioned fluid–structure interaction (FSI) simulation methodology relying on a finite-volume fluid solver for curvilinear, block-structured, body-fitted grids written in the Arbitrary Lagrangian–Eulerian (ALE) formulation, and a finite-element solver for the structural analysis. The second approach is an immersed boundary (IB) method employing a continuous and direct forcing strategy. The IB method, often applied to Cartesian grids, is also referred to as an approach to simulate fluid–structure interactions. In this study, both methods are combined to exploit their respective advantages in simulating a complex flow problem. The coupled FSI problem involves the interaction of a thin, flexible structure deforming under the dynamic load of a wind gust (task 1). The gust itself is generated by an artificial wind gust generator, which includes a <em>paddle</em> that partially obstructs the wind tunnel’s outlet, thereby defining an FSI problem of its own (task 2). For task 1, the classical partitioned ALE approach is employed, while the IB method is more appropriate for task 2. Using available experimental measurement data for both the fluid flow and the structural deformation, the combined simulation framework is first validated for the case without gust. In a second step, the more challenging FSI problem of discrete gusts impacting the T-structure is thoroughly analyzed and the predicted data are compared with the available measurement data. For both cases without and with gusts, a very good agreement between simulation and experiment is achieved, which justifies the chosen approach.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"294 ","pages":"Article 106621"},"PeriodicalIF":2.5,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143799504","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Computers & FluidsPub Date : 2025-04-03DOI: 10.1016/j.compfluid.2025.106604
Haowei Liu, Zhiyong Liu, Qiuyan Xu, Jiye Yang
{"title":"Hierarchical radial basis functions method for solving the unsteady Navier–Stokes equations","authors":"Haowei Liu, Zhiyong Liu, Qiuyan Xu, Jiye Yang","doi":"10.1016/j.compfluid.2025.106604","DOIUrl":"10.1016/j.compfluid.2025.106604","url":null,"abstract":"<div><div>The Navier–Stokes equations (NSE) are essential equations in fluid dynamics that describe the motion of viscous fluids, accurately reflecting changes in fluid velocity and pressure. It is widely used in the fields of aerodynamics analysis in aerospace, fluid flow simulation and others. We use the hierarchical radial basis functions (H-RBFs) collocation method to simulate NSE in this paper, which is essentially a meshfree method that only requires knowledge of scattered data node information without the need to build the meshgrid. Then, the trial space of H-RBFs is constructed with the help of nested sets of points and scaling the support radii of compactly supported radial basis functions. Meanwhile, the numerical solution is found within the trial space to approximate the model’s solution. Numerical tests demonstrate that the method proposed in this paper achieves high accuracy in both regular and irregular domains. Compared to compactly supported radial basis functions collocation method, H-RBFs collocation method exhibits smaller errors, particularly when the collocation points become dense. Finally, a classical experiment above flow around the cylinder is presented, demonstrating that H-RBFs collocation method can effectively simulate the formation of vortex streets.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106604"},"PeriodicalIF":2.5,"publicationDate":"2025-04-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143785978","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}
Computers & FluidsPub Date : 2025-03-26DOI: 10.1016/j.compfluid.2025.106609
Ksenia Kozhanova , Yannick Hoarau , Eric Goncalves da Silva
{"title":"A 3D numerical strategy for the computations of shock-induced bubble collapse near a wall","authors":"Ksenia Kozhanova , Yannick Hoarau , Eric Goncalves da Silva","doi":"10.1016/j.compfluid.2025.106609","DOIUrl":"10.1016/j.compfluid.2025.106609","url":null,"abstract":"<div><div>The importance of modelling two-phase flows involving shock waves arises from many engineering and medical applications. The presence of strong shock waves and their interactions with bubble interfaces, the high density ratio between phases and the large variation of material properties makes the resolution of such problems a complicated task for the numerical methods. While the variety of numerical techniques to solve these problems exist, e.g. the sharp interface or the diffuse interface methods, these strategies can lead to spurious oscillations of the solution near the interface. It is well known that it is difficult to achieve both a high order accuracy of the scheme and the monotonicity of the solution. In this paper a four-equation two-phase model is employed and integrated in an explicit fully parallelised finite-volume solver with HLLC numerical scheme coupled with WENO reconstruction methods and Hancock predictor–corrector scheme and non-uniform mesh based on stretching function in order to compute a 3D shock-induced bubble collapse near a wall. The novelty of our work is improved accuracy of computations of such a problem with optimised computational cost thanks to the non-uniform mesh introduction in 3D computations.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106609"},"PeriodicalIF":2.5,"publicationDate":"2025-03-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143746651","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}
Computers & FluidsPub Date : 2025-03-24DOI: 10.1016/j.compfluid.2025.106608
P. Olivucci , X. Shao , M. Albers , W. Schröder , R. Semaan
{"title":"Multi-fidelity modeling of net power savings of an actuated turbulent boundary layer","authors":"P. Olivucci , X. Shao , M. Albers , W. Schröder , R. Semaan","doi":"10.1016/j.compfluid.2025.106608","DOIUrl":"10.1016/j.compfluid.2025.106608","url":null,"abstract":"<div><div>We simulate and model the net power savings of an actuated turbulent boundary layer flow. The actuation is performed by spanwise traveling transverse surface waves parametrized by wavelength, amplitude, and period. The data is provided by 81 large-eddy simulations (LES) over a range of conditions. Since the numerical resolution of the skin-friction physics requires expensive large-eddy simulations, additional input power data is provided by low-cost, low-fidelity, two-dimensional simulations. An ad-hoc Gaussian Process (GP) framework is used to construct a single and a multi-fidelity surrogate model of the net power savings response to a range of actuation settings. The multi-fidelity model is shown to be able to leverage the two databases and combine the two independent constitutive models for the drag reduction and the input power. The predictive performance of the model is evaluated and compared to the single-fidelity baseline through cross-validated accuracy scores, including its probabilistic predictions. The models are queried to infer the detailed dependence of the flow response on the control parameters, to explore the existence of maxima, and to discuss the physical underpinnings of the flow.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106608"},"PeriodicalIF":2.5,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143725253","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}
Computers & FluidsPub Date : 2025-03-24DOI: 10.1016/j.compfluid.2025.106607
Nathaniel Hildebrand , Meelan M. Choudhari , Preethi V. Mysore , Balaji S. Venkatachari , Pedro Paredes
{"title":"Zonal and adaptive mesh refinement for RANS-based transition models in FUN3D","authors":"Nathaniel Hildebrand , Meelan M. Choudhari , Preethi V. Mysore , Balaji S. Venkatachari , Pedro Paredes","doi":"10.1016/j.compfluid.2025.106607","DOIUrl":"10.1016/j.compfluid.2025.106607","url":null,"abstract":"<div><div>We determine how the grid resolution and topology influence the accuracy and convergence of RANS-based transition model results by analyzing the NLF-0416 and NLR-7301 airfoils at subsonic and transonic Mach numbers, respectively. Natural and separation-induced transition scenarios are analyzed using the Langtry-Menter <span><math><mi>γ</mi></math></span>-<span><math><msub><mrow><mi>Re</mi></mrow><mrow><msub><mrow><mi>θ</mi></mrow><mrow><mi>t</mi></mrow></msub></mrow></msub></math></span> model. Multiple grid refinement techniques are investigated. First, we determine the relative effectiveness of zonal streamwise refinement near transition to turbulence for structured grids as an alternative to costly global uniform refinement. We also complement this zonal technique by globally varying the wall-normal resolution keeping the streamwise resolution fixed. The zonal streamwise refinement can accurately model separation-induced transition, but significant wall-normal resolution is needed to model natural transition that occurs on airfoils. A series of unstructured prismatic grids that have similar node counts and viscous wall spacings as the structured hexahedral grids result in solutions that are only about 1% different in terms of the lift and drag coefficients at infinite resolution according to Richardson extrapolation. We employ Mach-Hessian-based unstructured grid adaptation to natural and separation-induced transition on the NLF-0416 airfoil, which leads to both the lift and drag coefficients plateauing on coarse grids, but the converged transition locations can be inaccurate due to poor near-wall resolution. Adjoint-based grid adaptation is explored briefly, and for imposed transition with the Spalart–Allmaras model, it yields accurate solutions even for coarse grid resolutions.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106607"},"PeriodicalIF":2.5,"publicationDate":"2025-03-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143705011","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}
Computers & FluidsPub Date : 2025-03-21DOI: 10.1016/j.compfluid.2025.106605
Jiaxian Qin , Yaming Chen , Xiaogang Deng
{"title":"Unified algorithm to efficiently implement high-order schemes of differential form on unstructured grids","authors":"Jiaxian Qin , Yaming Chen , Xiaogang Deng","doi":"10.1016/j.compfluid.2025.106605","DOIUrl":"10.1016/j.compfluid.2025.106605","url":null,"abstract":"<div><div>It is well recognized that structured grid methods have the advantage of significantly lower computational cost needed for achieving higher accuracy, and unstructured grid methods are welcomed for their convenience in grid generation. Nevertheless, the major obstacle hindering the practical application of structured grid methods to complex engineering problems is the generation of high-quality grids, whereas the bottleneck for high-order unstructured grid methods lies in their huge demand for computational resource and memory. In this work, we propose a novel algorithm which enables the efficient implementation of high-order schemes of differential form on unstructured grids. By dividing the initial simplex cells locally into quadrilateral or hexahedral sub-cells, unique line-structures called Hamiltonian path can be identified. Subsequently, high-order spacial discretization can be done in a dimension-by-dimension manner along the lines, inheriting the accuracy and efficiency of structured grid methods. Meanwhile, the line-structures make the application of robust line-implicit time-marching schemes on fully unstructured grid possible, leading to excellent efficiency and robustness.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106605"},"PeriodicalIF":2.5,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143681860","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}
Computers & FluidsPub Date : 2025-03-21DOI: 10.1016/j.compfluid.2025.106602
Fandi D. Suprianto , Ming-Jyh Chern , Chin-Cheng Wang
{"title":"An enhanced analytical-based geometry processor for the volume of solid (VoS) method in CPU and GPU computations","authors":"Fandi D. Suprianto , Ming-Jyh Chern , Chin-Cheng Wang","doi":"10.1016/j.compfluid.2025.106602","DOIUrl":"10.1016/j.compfluid.2025.106602","url":null,"abstract":"<div><div>This paper introduces an analytical–based geometry processor which improves the efficiency of volume-of-solid (VoS) CFD solvers for fluid–structure interaction (FSI) scenarios with rigid body motion. The combination of the VoS method and the direct forcing immersed boundary (DFIB) method speeds up geometry creation while significantly reducing computational workload during unsteady simulations of solid structure movement or deformation. The proposed geometry processor’s effectiveness is demonstrated in the vortex-induced vibrations (VIV) scenario of a single circular cylinder, where the geometry was dynamically modified at each timestep to correspond to the cylinder’s passive motion. Comparisson studies show that the enhanced VoS function accelerates geometry construction process significantly faster than the ray-casting method while maintaining the same level of accuracy. Furthermore, its GPU implementation consistently achieves significant speedup across various computational loads, indicating superior scalability. This proposed method not only works well for simple geometries, but it also supports a wide range of multi-body curvilinear forms, including extruded spans and fully defined 3D geometries, making it suitable for a variety of applications.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"293 ","pages":"Article 106602"},"PeriodicalIF":2.5,"publicationDate":"2025-03-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143705012","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}