Computers & Fluids最新文献_第10页

An Euler–Lagrange coupling strategy for pressurized liquid jet atomization in Large Eddy Simulation 大涡模拟中高压液体射流雾化的欧拉-拉格朗日耦合策略

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-11 DOI: 10.1016/j.compfluid.2025.106625

Thibault Gioud , Thomas Laroche , Thomas Schmitt , Bénédicte Cuenot , Odier Nicolas

{"title":"An Euler–Lagrange coupling strategy for pressurized liquid jet atomization in Large Eddy Simulation","authors":"Thibault Gioud , Thomas Laroche , Thomas Schmitt , Bénédicte Cuenot , Odier Nicolas","doi":"10.1016/j.compfluid.2025.106625","DOIUrl":"10.1016/j.compfluid.2025.106625","url":null,"abstract":"<div><div>The atomization process is crucial for numerous industrial and engineering applications, ranging from fuel injection systems to spray coating technologies and remains a challenge for numerical prediction. This study presents a coupled Euler/Lagrange methodology to model the atomization processes, at affordable numerical costs. An Eulerian framework is used to resolve the liquid core and the largest liquid structures, while a Lagrangian approach tracks the smallest droplets and their interactions, addressing transport, secondary breakup, and evaporation phenomena. The Eulerian two-phase flow method considered here is a multi-fluid diffuse interface approach that assumes equilibrium of temperature, pressure, velocity, and Gibbs potentials at the liquid–gas interface. An originality of this paper is the transfer of liquid from the under resolved Eulerian liquid structures to the Lagrangian formalism when user-defined geometrical criteria are met. A validation of the coupling strategy is performed on a jet in cross flow configuration. Although the spatial distribution of the droplet volumetric flux is not exactly retrieved, other key characteristics of the spray, such as the Sauter Mean Diameter (SMD) and the outer trajectory, show very encouraging results. Finally, this paper shows the importance of taking into account the liquid core in the simulation, rather than directly injecting already atomized droplets.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"294 ","pages":"Article 106625"},"PeriodicalIF":2.5,"publicationDate":"2025-04-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143835227","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

Application of the quantum Fourier transform in a harmonic balance solver for Burgers’ equation 量子傅立叶变换在Burgers方程谐波平衡解算器中的应用

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-09 DOI: 10.1016/j.compfluid.2025.106619

Loïc Dewitte , Jérémie Roland , Frank Eulitz

{"title":"Application of the quantum Fourier transform in a harmonic balance solver for Burgers’ equation","authors":"Loïc Dewitte , Jérémie Roland , Frank Eulitz","doi":"10.1016/j.compfluid.2025.106619","DOIUrl":"10.1016/j.compfluid.2025.106619","url":null,"abstract":"<div><div>This paper explores the potential of quantum computing for computational fluid dynamics (CFD) applications, with a focus on turbomachinery CFD. A harmonic balance solver is developed for Burgers’ equation, in which discrete Fourier transforms are approximated using a hybrid quantum-classical algorithm based on the quantum Fourier transform (QFT). Three novel algorithms are presented to estimate Fourier coefficients, providing complete knowledge of their amplitudes and phases, which is not possible with the standard QFT. Their behaviour is studied theoretically and numerically, under noiseless and noisy conditions. The algorithms, whose performance is limited by the extensive sampling required to achieve a small error on the Fourier coefficients, tolerate low levels of depolarising noise. The behaviour of the hybrid solver is investigated for a baseline case with a Reynolds number of 1000, 7 harmonics and 100 grid cells, using the best performing algorithm in a noiseless setting with up to 10<sup>8</sup> samples per QFT. Residuals decrease until errors introduced by statistical uncertainty dominate the total error on the solution. Nevertheless, the hybrid solutions match the classical one closely, with RMS residuals as low as 10<sup>–4</sup>. The impact of several solver parameters on convergence and solution quality is also assessed, including the effect of noise. Although the latter rapidly degrades the solution, the solver achieves satisfactory approximations to the classical solution with 0.001% of depolarising noise. Without seeking to demonstrate a quantum advantage, this work offers valuable insights into the opportunities and challenges of quantum computing, helping readers understand how to design, implement and study quantum algorithms, as well as evaluate their impact in CFD applications.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"295 ","pages":"Article 106619"},"PeriodicalIF":2.5,"publicationDate":"2025-04-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"143860235","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

Locally energy-stable finite element schemes for incompressible flow problems: Design and analysis for equal-order interpolations 不可压缩流动问题的局部能量稳定有限元方案：等阶插值的设计与分析

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-09 DOI: 10.1016/j.compfluid.2025.106622

Hennes Hajduk , Dmitri Kuzmin , Gert Lube , Philipp Öffner

引用次数: 0

Divergence cleaning for weakly compressible smoothed particle hydrodynamics 弱可压缩光滑粒子流体力学的散度清理

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-08 DOI: 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}

引用次数: 0

Large Eddy simulation of a supersonic starting impinging jet 超音速起动撞击射流的大涡模拟

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-05 DOI: 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}

引用次数: 0

Modeling two-phase equilibrium using smoothed particle hydrodynamics 用光滑粒子流体力学建模两相平衡

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-05 DOI: 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}

引用次数: 0

Combination of two FSI methods and their validation based on artificial wind gusts impacting a flexible T-structure 两种FSI方法的结合及其基于人工阵风对柔性t型结构影响的验证

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-03 DOI: 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}

引用次数: 0

Hierarchical radial basis functions method for solving the unsteady Navier–Stokes equations 求解非定常Navier-Stokes方程的分层径向基函数法

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-04-03 DOI: 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}

引用次数: 0

A 3D numerical strategy for the computations of shock-induced bubble collapse near a wall 一种计算近壁激波气泡崩塌的三维数值策略

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-03-26 DOI: 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}

引用次数: 0

Multi-fidelity modeling of net power savings of an actuated turbulent boundary layer 驱动湍流边界层净功率节约的多保真度建模

IF 2.5 3区工程技术

Computers & Fluids Pub Date : 2025-03-24 DOI: 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}

引用次数: 0