Computers & Fluids最新文献_第4页

Development of ISPH-FVM coupling method by embedded MHD model to simulate bubble rising in conductive liquids under magnetic field 基于嵌入式MHD模型的ISPH-FVM耦合方法在磁场作用下模拟导电液体气泡上升

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-05 DOI: 10.1016/j.compfluid.2025.106941

Yixiang Xu , Gang Yang , Yulin Xing , Dean Hu

{"title":"Development of ISPH-FVM coupling method by embedded MHD model to simulate bubble rising in conductive liquids under magnetic field","authors":"Yixiang Xu , Gang Yang , Yulin Xing , Dean Hu","doi":"10.1016/j.compfluid.2025.106941","DOIUrl":"10.1016/j.compfluid.2025.106941","url":null,"abstract":"<div><div>The motion of bubbles in conductive fluid under the external magnetic field is affected by multiple physical fields, so that the flow characteristics of bubbles are usually difficult to predict, which has gradually become a hot topic in magnetohydrodynamics research. In this paper, the dynamic behavior of bubble rising in conductive liquids under magnetic field is simulated by three-dimensional (3D) ISPH-FVM coupling method. The bubble motion is captured by the ISPH particles, while the FVM mesh is used to calculate the flow field. For the combination of magnetohydrodynamics model and ISPH-FVM coupling method, FVM grid is used to solve Maxwell equations combining electric field and magnetic field, and the calculated Lorentz force is introduced into the momentum equation. To verify the accuracy and stability of the present ISPH-FVM coupling model, the Shercliff case, bubble rising under horizontal magnetic field and two horizontal bubbles coalescence are tested. Subsequently, the single bubble rising under different magnetic field directions and intensities is simulated, and the dynamic mechanism behind the morphological changes and rising velocity of the bubble is deeply analyzed. Finally, the influence of different magnetic field directions and intensities on the terminal velocity and bubble shape of rising bubble with different sizes are discussed.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106941"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145787872","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

Curve boundary algorithms based on interpolation for multiphase lattice Boltzmann method 基于插值的多相晶格玻尔兹曼法曲线边界算法

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-26 DOI: 10.1016/j.compfluid.2025.106954

Dingyu Song , Fu Ling , Yonggang Zhang , Binghai Wen

引用次数: 0

Numerical investigation of flow induced by plasma actuators around a circular cylinder in quiescent air under duty-cycle actuation 占空比作用下静空气中等离子体作动器绕圆柱体流动的数值研究

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-24 DOI: 10.1016/j.compfluid.2025.106953

Zihao Zhu, Pau Fradera-Soler, Yalu Zhu, Feng Liu

{"title":"Numerical investigation of flow induced by plasma actuators around a circular cylinder in quiescent air under duty-cycle actuation","authors":"Zihao Zhu, Pau Fradera-Soler, Yalu Zhu, Feng Liu","doi":"10.1016/j.compfluid.2025.106953","DOIUrl":"10.1016/j.compfluid.2025.106953","url":null,"abstract":"<div><div>The flow induced by a pair of Dielectric-Barrier-Discharge (DBD) plasma actuators symmetrically-mounted on a circular cylinder in quiescent air is simulated by solving the Unsteady Reynolds-Averaged Navier-Stokes (URANS) equations. A new body-force model for the plasma actuator is developed, which yields significantly improved agreement with experiments. In addition to velocity and vorticity fields, the computations provide time-instantaneous and time-averaged pressure field and pressure and skin friction distributions over the cylinder, which have not been available from existing experiments for this problem. The computational results are analyzed to reveal the relations between the plasma body force, pressure field, velocity, vorticity, forces on the cylinder, and momentum imparted to the flow as they impact the flow control effectiveness for different duty-cycle ratio and frequency of the plasma actuation input. A pair of cumulative vortices is discovered outside the experimental measurement window for the first time by the computations. Regardless of the level of duty-cycle frequency and the detailed near field vortex patterns, the vortices appear to eventually accumulate some distance downstream of the cylinder and move very slowly downstream. The present study is of direct relevance to problems involving vortex shedding and noise production from many circular-shaped aeronautical and civil structures. It helps to provide the necessary fundamental understanding of the flow physics and guidance for future design and optimization of DBD plasma actuators for flow control.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106953"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880628","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

An all-Mach cell-centered multi-dimensional finite volume numerical scheme for the Euler equations 欧拉方程的全马赫中心多维有限体积数值格式

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-25 DOI: 10.1016/j.compfluid.2025.106951

Alessia Del Grosso , Wasilij Barsukow , Raphaël Loubère , Pierre-Henri Maire

{"title":"An all-Mach cell-centered multi-dimensional finite volume numerical scheme for the Euler equations","authors":"Alessia Del Grosso , Wasilij Barsukow , Raphaël Loubère , Pierre-Henri Maire","doi":"10.1016/j.compfluid.2025.106951","DOIUrl":"10.1016/j.compfluid.2025.106951","url":null,"abstract":"<div><div>In the context of the numerical approximation of Euler equations, great efforts have been devoted to developing schemes that can accurately reproduce solutions in low Mach number flows. Solutions of classic Finite Volume (FV) schemes are usually plagued by an excessive diffusion as the numerical scheme is not consistent with the limit equations for the Mach number that tends to zero. Instead, a numerical scheme that satisfies such a property is called Asymptotic-Preserving (AP). In this paper, we propose an AP FV scheme for the multi-dimensional Euler equations. In classic FV methods, the numerical approximation of the face flux is obtained by means of a two-state 1D approximate Riemann Solver (RS) in the normal direction to the face. Here, we rely on a node-based flux approximation that exploits a particular RS involving a nodal quantity which depends on all the cells around a given node. Such an idea has been exploited by Barsukow et al. (2023) for the linear acoustic equations. Their method is vorticity-preserving, but its extension to the Euler equations proved to be far from trivial. For such a reason, a change of perspective is needed in the definition of the RS.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106951"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880043","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

Fine-tuning physics-informed neural networks for cavity flows using coordinate transformation 利用坐标变换对腔体流动的物理信息神经网络进行微调

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-27 DOI: 10.1016/j.compfluid.2025.106957

Ryuta Takao , Satoshi Ii

{"title":"Fine-tuning physics-informed neural networks for cavity flows using coordinate transformation","authors":"Ryuta Takao , Satoshi Ii","doi":"10.1016/j.compfluid.2025.106957","DOIUrl":"10.1016/j.compfluid.2025.106957","url":null,"abstract":"<div><div>Physics-informed neural networks (PINNs) have attracted attention as an alternative approach to solve partial differential equations using a deep neural network (DNN). Their simplicity and capability allow them to solve inverse problems for many applications. Despite the versatility of PINNs, it remains challenging to reduce their training cost. Using a DNN pre-trained with an arbitrary dataset with transfer learning or fine-tuning is a potential solution. However, a pre-trained model using a different geometry and flow condition than the target may not produce suitable results. This paper proposes a fine-tuning approach for PINNs with coordinate transformation, modelling lid-driven cavity flows with various shapes. We formulate the inverse problem, where the reference data inside the domain and wall boundary conditions are given. A pre-trained PINN model with an arbitrary Reynolds number and shape is used to initialize a target DNN. To reconcile the reference shape with different targets, governing equations as a loss of the PINNs are given with coordinate transformation using a deformation gradient tensor. Numerical examples for various cavity flows with square, rectangular, shear deformed and inflated geometries demonstrate that the proposed fine-tuning approach improves the training convergence compared with a randomly-initialized model. A pre-trained model with a similar geometry to the target further increases training efficiency. These findings are useful for real-world applications such as modelling intra-aneurysmal blood flows in clinical use.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106957"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145880627","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

Artificial viscosity based on the velocity gradient components in a 3D Lagrangian Lax-Wendroff scheme 三维拉格朗日Lax-Wendroff格式中基于速度梯度分量的人工粘度

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-13 DOI: 10.1016/j.compfluid.2025.106944

Matěj Klíma, Milan Kuchařík, Richard Liska

引用次数: 0

Aerodynamic shape optimization with discretization error control for high-order compressible flow simulations 高阶可压缩流动模拟的离散误差控制气动形状优化

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-16 DOI: 10.1016/j.compfluid.2025.106949

Kush Pandya, Aravind Balan

{"title":"Aerodynamic shape optimization with discretization error control for high-order compressible flow simulations","authors":"Kush Pandya, Aravind Balan","doi":"10.1016/j.compfluid.2025.106949","DOIUrl":"10.1016/j.compfluid.2025.106949","url":null,"abstract":"<div><div>The shape of an aerospace vehicle significantly influences its performance, and Aerodynamic Shape Optimization (ASO) refines the shape to achieve objectives such as drag reduction and lift improvement. ASO typically combines Computational Fluid Dynamics (CFD) simulations with gradient-based optimization. Adjoint methods can compute gradients independently of the number of shape parameters, making it computationally efficient when optimizing for a large number of shape parameters. For high reliability on the results of ASO, we need accurate evaluation of the objective functions and their gradients. High-order numerical methods, such as the Discontinuous Galerkin (DG) method, provide superior accuracy for compressible flow simulations at a computational cost comparable to traditional methods, making them ideal for ASO. The present study highlights the impact of discretization errors on optimal geometry, emphasizing the need for controlling them in an efficient way during optimization. A novel strategy is proposed that efficiently integrates adjoint-based mesh optimization within the ASO framework that uses a high-order DG method for the flow simulations. The approach ensures that discretization errors in the objective function remain bounded throughout the ASO process by dynamically adapting the meshes based on adjoint-based error estimates. The mesh adaptation is formulated in such a way that one gets optimum meshes with minimal number of elements for a prescribed error tolerance. The strategy results in a multi-fidelity approach where there is an effective utilization of computation cost by progressively refining the mesh as ASO progresses. The effectiveness of the proposed methodology is demonstrated through a comparative analysis with conventional ASO performed on a fixed fine mesh, using a benchmark test case from the AIAA Aerodynamic Design Optimization Discussion Group (ADODG). The results highlight significant improvements in the optimization process’s reliability and computational efficiency.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106949"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145836957","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 Markov Matrix iterative splitting algorithm for natural convection 自然对流的马尔可夫矩阵迭代分裂算法

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-05 DOI: 10.1016/j.compfluid.2025.106933

Manira E. Narvaez-Saucedo, Carlos A. López-Villalobos, Manuel A. Ramírez-Cabrera, Eduardo Ramos, Patricio Javier Valades-Pelayo

{"title":"A Markov Matrix iterative splitting algorithm for natural convection","authors":"Manira E. Narvaez-Saucedo, Carlos A. López-Villalobos, Manuel A. Ramírez-Cabrera, Eduardo Ramos, Patricio Javier Valades-Pelayo","doi":"10.1016/j.compfluid.2025.106933","DOIUrl":"10.1016/j.compfluid.2025.106933","url":null,"abstract":"<div><div>This work introduces a novel Markov Matrix Iterative Splitting Algorithm (MISA), a probabilistic operator splitting Monte Carlo approach designed for simulating natural convection in incompressible flows within the continuum regime. This approach effectively captures the nonlinear convective transport and stochastic diffusion while enforcing global incompressibility, without requiring mesh generation or matrix inversions common to deterministic solvers. Variance reduction techniques and a Multi-Level Monte Carlo framework improve convergence and stability. MISA naturally handles complex geometries with internal obstacles, incorporated effortlessly by excluding walkers inside them and imposing boundary conditions on their surfaces. Benchmark natural convection simulations in square cavities demonstrate MISA’s accuracy and computational efficiency compared to traditional finite-volume solvers. The method reliably resolves convective flow structures and temperature fields, highlighting its promising potential for complex fluid-thermal transport problems in irregular domains</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106933"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145920697","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

Non-orthogonal plane-marching parabolized stability equations for the secondary instability of crossflow vortices 横流涡旋二次不稳定性的非正交平面推进抛物化稳定性方程

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-17 DOI: 10.1016/j.compfluid.2025.106947

Koen J. Groot , Jordi Casacuberta , Stefan Hickel

{"title":"Non-orthogonal plane-marching parabolized stability equations for the secondary instability of crossflow vortices","authors":"Koen J. Groot , Jordi Casacuberta , Stefan Hickel","doi":"10.1016/j.compfluid.2025.106947","DOIUrl":"10.1016/j.compfluid.2025.106947","url":null,"abstract":"<div><div>A detailed derivation, analysis, and verification is given for the non-orthogonal, plane-marching Parabolized Stability Equations (PSE) approach. In applying the approach to a flow distorted by a medium-amplitude crossflow vortex, we determine its linear secondary instability mechanisms. We show that converged solutions can be achieved for a broad frequency range with an existing stabilization method for the line-marching PSE approach. We verify that 1) solutions converge versus grid size in all dimensions, 2) primary disturbance solutions agree with line-marching PSE results, and 3) secondary disturbance solutions match amplitude and growth-rate evolution of reference Direct Numerical Simulation (DNS) results. We show how and why the <em>type</em>-II instability displays a delayed neutral point when modeled with the plane-marching approach versus the considered local stability approaches, whether the streamwise evolution of the distorted base flow is accounted for or not. This may explain why the <em>type</em>-II disturbance is scarcely captured by DNS in the literature.</div></div>","PeriodicalId":287,"journal":{"name":"Computers & Fluids","volume":"306 ","pages":"Article 106947"},"PeriodicalIF":3.0,"publicationDate":"2026-02-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"145836980","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

High-order gas-kinetic scheme with two-derivative-based time discretization for compressible flows 可压缩流动的高阶气体动力学双导数时间离散格式

IF 3 3区工程技术

Computers & Fluids Pub Date : 2026-02-15 Epub Date: 2025-12-18 DOI: 10.1016/j.compfluid.2025.106950

Qiushuo Qin , Jie Wu , Lan Jiang

引用次数: 0