International Journal for Numerical Methods in Fluids最新文献

{"title":"Description and Convergence Order Analysis of the Finite Element-Volume Spatial Discretization Method","authors":"Maria Adela Puscas,&nbsp;Pierre-Emmanuel Angeli,&nbsp;Nathalie Nouaime,&nbsp;Erell Jamelot","doi":"10.1002/fld.5399","DOIUrl":"https://doi.org/10.1002/fld.5399","url":null,"abstract":"<p>This article reviews the Finite Element-Volume spatial discretization method on tetrahedral meshes implemented in the TrioCFD code. TrioCFD is a computational fluid dynamics software specifically designed for simulating turbulent flows and heat transfer, with a primary focus on nuclear engineering applications. Its versatility also extends to a wide range of engineering applications. This article presents the principles of Finite Element-Volume discretization and conducts an analysis of its properties and convergence orders. The time scheme can be explicit or semi-implicit, and the velocity and pressure are updated based on a projection-correction scheme. The discretization ensures local mass conservation, second-order convergence for velocity, and first-order convergence for pressure. Moreover, a super-convergence is achieved in two-dimension when the source term is a gradient. The accuracy and convergence rate of the numerical method is rigorously assessed in a variety of two- and three-dimensional test cases, including steady Stokes and Navier–Stokes problems with manufactured solutions. The classical lid-driven cavity flow is also addressed, and the results of comparisons between the solutions from the present numerical method and the results of reference literature are provided.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1189-1208"},"PeriodicalIF":1.8,"publicationDate":"2025-05-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5399","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767488","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"PolyMAC: A Staggered Finite Volume Method on General Meshes for Incompressible Navier-Stokes Equations","authors":"Pierre-Loïc Bacq,&nbsp;Antoine Gerschenfeld,&nbsp;Michael Ndjinga","doi":"10.1002/fld.5398","DOIUrl":"https://doi.org/10.1002/fld.5398","url":null,"abstract":"<div>\u0000 \u0000 <p>We consider the numerical resolution of the incompressible Navier-Stokes equations. We present a new compatible Finite Volume discretisation that generalises the famous Marker-and-Cell (MAC) method to polyhedral meshes that we call PolyMAC. In the first part of the paper, we recall the principles of compatible schemes and detail the key operators of the discretisation. The convergence and robustness of PolyMAC are assessed numerically first on a benchmark from the FVCA conferences. We then consider a problem of industrial complexity that allows us to confirm the robustness of PolyMAC on more complex problems. The second part of the article is dedicated to the efficient numerical resolution of the resulting linear system. Concretely, we use a PISO-like prediction-correction approach and develop efficient preconditioners for linear systems. In particular, we show that the saddle-point system arising from the correction step is very challenging for iterative methods on distorted meshes. In this work, we develop a robust preconditioner based on an algebraic transformation of the system. In particular, this new preconditioner shows impressive convergence on problems of industrial complexity.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1171-1188"},"PeriodicalIF":1.8,"publicationDate":"2025-04-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144768087","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Retraction: An Efficient Parallel Algorithm for Three-Dimensional Analysis of Subsidence Above Gas Reservoirs","authors":"","doi":"10.1002/fld.5396","DOIUrl":"https://doi.org/10.1002/fld.5396","url":null,"abstract":"<p>\u0000 <span>B.A. Schrefler</span>, <span>X. Wang</span>, <span>V.A. Salomoni</span>, and <span>G. Zuccolo</span>, “ <span>An Efficient Parallel Algorithm for Three-Dimensional Analysis of Subsidence Above Gas Reservoirs</span>,” <i>International Journal for Numerical Methods in Fluids</i> <span>31</span>, no. <span>1</span> (<span>1999</span>): <span>247</span>–<span>260</span>, \u0000https://doi.org/10.1002/(SICI)1097-0363(19990915)31:1&lt;247::AID-FLD966&gt;3.0.CO;2-D.</p><p>The above article, published online on 14 September 1999 in Wiley Online Library (wileyonlinelibrary.com), has been retracted by agreement between the authors; the journal Editor-in-Chief, Alina Bruma; and John Wiley &amp; Sons Ltd. The retraction has been agreed due to the authors' discovery that the proper permissions for use of Table 1 and Figures 4, 5, 6, 7, and 8 were not obtained prior to publication. As it was not possible to obtain retrospective permission, the article must therefore be retracted.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 8","pages":""},"PeriodicalIF":1.7,"publicationDate":"2025-04-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5396","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144525185","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Vortex Panel Method in Axisymmetric Cylindrical Coordinates: Inviscid Formulation","authors":"Suguru Shiratori,&nbsp;Kosuke Kimata,&nbsp;Masaya Katoh,&nbsp;Hideaki Nagano,&nbsp;Kenjiro Shimano","doi":"10.1002/fld.5400","DOIUrl":"https://doi.org/10.1002/fld.5400","url":null,"abstract":"<p>This study addresses the vortex panel method, which is an efficient solution for inviscid flow around a body. For Cartesian coordinates, the panel method has been developed and widely applied. However, a formulation has not been proposed for axisymmetric cylindrical coordinates. This study derives a Green's function corresponding to the governing equation of the potential flow in axisymmetric cylindrical coordinates by modifying the Green's function provided by Cohl and Tohline [doi:10.1086/308062]. The derived Green's function <span></span><math>\u0000 \u0000 <semantics>\u0000 <mrow>\u0000 <mi>G</mi>\u0000 </mrow>\u0000 <annotation>$$ G $$</annotation>\u0000 </semantics></math> is significantly compact and is composed of a single term of the half-integer degree Legendre function of the second type. The derivation of Green's function <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>G</mi>\u0000 </mrow>\u0000 <annotation>$$ G $$</annotation>\u0000 </semantics></math> is confirmed by analytically evaluating the requirement <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℒ</mi>\u0000 <mi>G</mi>\u0000 <mo>=</mo>\u0000 <mn>0</mn>\u0000 </mrow>\u0000 <annotation>$$ mathit{mathcal{L}G}=0 $$</annotation>\u0000 </semantics></math>, where <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>ℒ</mi>\u0000 </mrow>\u0000 <annotation>$$ mathcal{L} $$</annotation>\u0000 </semantics></math> is a linear operator of the governing equation. Under the derived <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>G</mi>\u0000 </mrow>\u0000 <annotation>$$ G $$</annotation>\u0000 </semantics></math>, the vortex panel method is formulated by discretizing the vorticity distribution along the body surface. The validity of the constructed panel method is confirmed through calculations of the flow past a sphere, an ellipsoid, a torus, and a teardrop-like object by comparison with analytical or numerical solutions.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1161-1170"},"PeriodicalIF":1.8,"publicationDate":"2025-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5400","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767933","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Comparative Analysis of Reinforcement Learning Agents for Optimizing Airfoil Shapes","authors":"Harsh H. Sawant,&nbsp;Rahul Gujar,&nbsp;Neeta Mandhare,&nbsp;M. J. Sable,&nbsp;Prashant K. Ambadekar,&nbsp;S. H. Gawande","doi":"10.1002/fld.5395","DOIUrl":"https://doi.org/10.1002/fld.5395","url":null,"abstract":"<div>\u0000 \u0000 <p>This work investigates the optimization of airfoil shapes using various reinforcement learning (RL) algorithms, including Deep Deterministic Policy Gradient (DDPG), Twin Delayed Deep Deterministic Policy Gradient (TD3), and Trust Region Policy Optimization (TRPO). The primary objective is to enhance the aerodynamic performance of airfoils by maximizing lift forces across different angles of attack (AoA). The study compares the optimized airfoils against the standard NACA 2412 airfoil. The DDPG-optimized airfoil demonstrated superior performance at lower and moderate AoAs, while the TRPO-optimized airfoil excelled at higher AoAs. In contrast, the TD3-optimized airfoil consistently underperformed. The results indicate that RL algorithms, particularly DDPG and TRPO, can effectively improve airfoil designs, offering substantial benefits in lift generation. This paper underscores the potential of RL techniques in aerodynamic shape optimization, presenting significant implications for aerospace and related industries.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 8","pages":"1142-1156"},"PeriodicalIF":1.7,"publicationDate":"2025-03-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144524550","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Modified CIP-Soroban Method and Its Application in Implosion Process of Inertial Confinement Fusion","authors":"Zhehao Lin,&nbsp;Kazumasa Takahashi,&nbsp;Toru Sasaki,&nbsp;Takashi Kikuchi,&nbsp;Atsushi Sunahara","doi":"10.1002/fld.5392","DOIUrl":"https://doi.org/10.1002/fld.5392","url":null,"abstract":"<div>\u0000 \u0000 <p>The CIP-Soroban method is an excellent adaptive meshless method capable of solving advection problems with 3rd-order accuracy by combining the Constrained Interpolation Profile/Cubic Interpolated Pseudo-particle (CIP) method. This study proposes a modified version of the CIP-Soroban method specifically designed to address severe compressible hydrodynamic scenarios. The proposed method includes a material distinguishing approach, incorporates a modified form of monitoring functions for grid generation, utilizes a staggered grid arrangement, incorporates the Maximum and minimum Bounds method, solves non-advection terms using a finite difference method, and employs an adjusted procedure for stably solving the governing equations. We applied the modified CIP-Soroban method to simulate the implosion process in inertial confinement fusion (ICF), which is commonly modeled by compressible fluid and has the problems of large gradients of physical values and strong nonlinearity for stable and accurate numerical analysis. Implosion simulations were performed using a series of grids with increasing resolutions, ranging from coarse to fine grid settings, as one of the application examples. The results indicated that compared to the conventional uniform grid CIP method, the modified CIP-Soroban method reduced computational costs (calculation time, memory occupancy, and grid number) for obtaining the same precision results.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 8","pages":"1120-1141"},"PeriodicalIF":1.7,"publicationDate":"2025-03-11","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144524883","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"A GPU Accelerated Mixed-Precision Finite Difference Informed Random Walker (FDiRW) Solver for Strongly Inhomogeneous Diffusion Problems","authors":"Zirui Mao,&nbsp;Shenyang Hu,&nbsp;Ang Li","doi":"10.1002/fld.5394","DOIUrl":"https://doi.org/10.1002/fld.5394","url":null,"abstract":"<div>\u0000 \u0000 <p>In nature, many complex multi-physics coupling problems exhibit significant diffusivity inhomogeneity, where one process occurs several orders of magnitude faster than others temporally. Simulating rapid diffusion alongside slower processes demands intensive computational resources due to the necessity for small time steps. To address these computational challenges, we have developed an efficient numerical solver named Finite Difference informed Random Walker (FDiRW). In this study, we propose a GPU-accelerated, mixed-precision configuration for the FDiRW solver to maximize efficiency through GPU multi-threaded parallel computation and lower precision computation. Numerical evaluation results reveal that the proposed GPU-accelerated mixed-precision FDiRW solver can achieve a 117× speedup over the CPU baseline, while an additional 1.75× speedup is achieved by employing lower precision GPU computation. Notably, for large model sizes, the GPU-accelerated mixed-precision FDiRW solver demonstrates strong scaling with the number of nodes used in simulation. When simulating radionuclide absorption processes by porous wasteform particles with a medium-sized model of 192 × 192 × 192, this approach reduces the total computational time to 10 min, enabling the simulation of larger systems with strongly inhomogeneous diffusivity.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 8","pages":"1104-1119"},"PeriodicalIF":1.7,"publicationDate":"2025-03-08","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144524736","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Non-Hydrostatic Model for Simulating Moving Bottom-Generated Waves: A Shallow Water Extension With Quadratic Vertical Pressure Profile","authors":"Kemal Firdaus,&nbsp;Jörn Behrens","doi":"10.1002/fld.5393","DOIUrl":"https://doi.org/10.1002/fld.5393","url":null,"abstract":"<p>We formulate a depth-averaged non-hydrostatic model to solve wave equations with generation by a moving bottom. This model is built on the shallow water equations, which are widely used in tsunami wave modelling. An extension leads to two additional unknowns to be solved: vertical momentum and non-hydrostatic pressure. We show that a linear vertical velocity assumption turns out to give us a quadratic pressure relation, which is equivalent to Boussinesq-type equations, the Green-Naghdi equations specifically, making it suitable for weakly dispersive cases. However, this extension involves a time derivative of an unknown parameter, rendering the solution by a projection method ambiguous. In this study, we derive an alternative form of the elliptic system of equations to avoid such ambiguity. The new set of equations satisfies the desired solubility property, while also consistently representing the non-flat moving topography wave generation. Validations are performed using several test cases based on the previous experiments and a high-fidelity simulation. First, we show the efficiency of our model in solving a vertical movement, which represents an undersea earthquake-generated tsunami. Following that, we demonstrate the accuracy of the model for landslide-generated waves. Finally, we compare the performance of our novel set of equations with the linear and simplified quadratic pressure profiles.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 8","pages":"1093-1103"},"PeriodicalIF":1.7,"publicationDate":"2025-03-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5393","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144524994","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"An Extended Height-Function Method for 3D VOF Simulations of Wetting Phenomena on Super-Hydrophilic and Hydrophobic Surfaces","authors":"Kenya Kitada,&nbsp;Ryoichi Kurose","doi":"10.1002/fld.5391","DOIUrl":"https://doi.org/10.1002/fld.5391","url":null,"abstract":"<div>\u0000 \u0000 <p>An extended height-function (HF) method that can be consistently utilized for 3D volume of fluid (VOF) simulations of wetting phenomena on super-hydrophilic and super-hydrophobic surfaces, is proposed. First, the standard HF method is briefly explained. Then, 2D and 3D HF methods that reflect the contact angles reported so far are described, with their limitations discussed. Finally, specific treatments of contact line identification and HF construction reflecting the contact angle boundary condition, required to overcome such limitations, are presented in detail. Numerical tests for a sessile droplet reveal that the contact line identification and HF construction are conducted appropriately with respect to the imposed contact angles ranging from <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>1</mn>\u0000 <msup>\u0000 <mrow>\u0000 <mn>5</mn>\u0000 </mrow>\u0000 <mrow>\u0000 <mo>∘</mo>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ 1{5}^{circ } $$</annotation>\u0000 </semantics></math> to <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>16</mn>\u0000 <msup>\u0000 <mrow>\u0000 <mn>5</mn>\u0000 </mrow>\u0000 <mrow>\u0000 <mo>∘</mo>\u0000 </mrow>\u0000 </msup>\u0000 </mrow>\u0000 <annotation>$$ 16{5}^{circ } $$</annotation>\u0000 </semantics></math> in the proposed numerical scheme. Additionally, the present method shows approximately first- or second-order convergence of the curvature at the contact line for a wide range of contact angles. Moreover, simulations of droplet spreading driven by surface tension reveal that the proposed method can reasonably reproduce the behavior of a droplet reaching an equilibrium state defined by an imposed contact angle.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 7","pages":"1073-1091"},"PeriodicalIF":1.7,"publicationDate":"2025-03-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144197161","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}

{"title":"Cavitation Implementation Algorithms Based on Pressure Projection Method for Incompressible Flows With Three-Phase Interactions","authors":"Zhaoyuan Wang,&nbsp;Sungtek Park,&nbsp;Frederick Stern","doi":"10.1002/fld.5389","DOIUrl":"https://doi.org/10.1002/fld.5389","url":null,"abstract":"<div>\u0000 \u0000 <p>In the present study, a cavitation implementation algorithm is developed using a pressure-based method for incompressible flows with three-phase interactions. Central to this implementation algorithm is the treatment of the velocity jump due to the phase change, which is included in both the cavitation transport and pressure equations. The velocity jump, as a function of the phase change rate, is added as a source term to the pressure Poisson equation. A non-conservative form of the vapor transport equation is derived, and the velocity divergence is replaced by a term related to the mass phase change rate. An algorithm for the three-phase (air, water, and vapor) interactions is also developed. The VOF method is modified and used to identify the ‘dry’ (air) phase and the ‘wet’ (water/vapor mixture) phase, since the cavitation can only occur inside the water phase. The liquid volume fraction is used to distinguish water and vapor phases. The numerical results of the 2D NACA66MOD and 3D Delft Twist 11 hydrofoils show good agreement with the experimental measurement. The forced unsteady cavitation flows are investigated using a pitching foil with the results compared with the experimental observations. Air–water interface effect on the cavitation is investigated using the NACA66MOD hydrofoil. The code is applied to simulate a surface piercing super cavitating hydrofoil with both ventilation and cavitation involved.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 7","pages":"1057-1072"},"PeriodicalIF":1.7,"publicationDate":"2025-03-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144197405","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":4,"RegionCategory":"工程技术","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}