International Journal for Numerical Methods in Fluids最新文献

{"title":"Three-Directional Orthogonality Preserving Method for Hyperbolic Grid Generation","authors":"Siyuan Pi,&nbsp;Hongyuan Lin,&nbsp;Shuyao Hu,&nbsp;Chongwen Jiang,&nbsp;Chun-Hian Lee","doi":"10.1002/fld.5407","DOIUrl":"https://doi.org/10.1002/fld.5407","url":null,"abstract":"<div>\u0000 \u0000 <p>The hyperbolic grid generation method is widely used for generating computational grids. Because of conflicts arising from various grid constraints, the traditional hyperbolic grid generation method faces challenges in guaranteeing the fulfillment of all orthogonal constraints among three directions during the grid generation. A new three-directional orthogonality preserving method (TDOP) is introduced in the present work to enhance the orthogonality of the computational grid during the grid generation process. Unlike the traditional grid generation method, TDOP takes all three orthogonal constraints into consideration, establishes a function to quantify the overall grid orthogonality, and subsequently derives new governing equations for grid generation by solving a constrained optimization problem. Compared with the traditional method, TDOP exhibits enhanced control over the orthogonality among three directions, thereby enabling the generation of a computational grid with better orthogonality. Three application cases are employed to demonstrate the effectiveness and superiority of TDOP in hyperbolic grid generation. Results indicate that, compared with the traditional method, TDOP can effectively prevent the emergence of highly skewed grids and enables enhanced optimization of orthogonality in the advancing front layer. Consequently, TDOP can generate a computational grid with better orthogonality and higher quality than the traditional method.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 10","pages":"1321-1342"},"PeriodicalIF":1.8,"publicationDate":"2025-06-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144927610","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":"Higher Order Method Based on the Combination of Local Discontinuous Galerkin and Spectral Deferred Correction Method for the Rotating Navier–Stokes Equations","authors":"Xuewei Zhang,&nbsp;Demin Liu","doi":"10.1002/fld.5408","DOIUrl":"https://doi.org/10.1002/fld.5408","url":null,"abstract":"<div>\u0000 \u0000 <p>In this article, the spatial local discontinuous Galerkin (LDG) method and the temporal spectral deferred correction (SDC) method are combined to construct the higher-order approximating method for the unsteady rotating Navier–Stokes equations on the triangular mesh. First, the artificially compressible method is used to circumvent the incompressibility constraint, and the rotating Navier–Stokes equations are transformed into the artificially compressible rotating Navier–Stokes equations. Then, based on equal LDG interpolation and repeated temporal SDC, the higher-order fully discrete method is presented. Theoretically, the stability analysis of the second-order fully discrete method is provided, and it is shown that the time step <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>τ</mi>\u0000 </mrow>\u0000 <annotation>$$ tau $$</annotation>\u0000 </semantics></math> is stable within the upper bound constraints. Numerical examples are presented to demonstrate the effectiveness of the proposed method.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1303-1318"},"PeriodicalIF":1.8,"publicationDate":"2025-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767639","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":"An Extended HLLD Riemann Solver for the Numerical Simulation of Magneto-Hydrodynamics","authors":"Xiaocheng Guo,&nbsp;Xinyue Xi,&nbsp;Chi Wang","doi":"10.1002/fld.5405","DOIUrl":"https://doi.org/10.1002/fld.5405","url":null,"abstract":"<div>\u0000 \u0000 <p>By revisiting the derivation of multi-state HLL approximate Riemann solver for the ideal magneto-hydrodynamics, an extended HLLD Riemann solver is constructed based on the assumption that the normal velocity is constant over the Riemann fan, which is bounded by two fast waves, and separated by two compound waves and a middle contact wave. Compared with the HLLD solver, the slow waves are allowed to persist inside the Riemann fan, so that the two Alfvén waves are replaced by the two compound waves that are the merging product of the Alfvén and slow waves. Conseq uently, the corresponding wave speeds are chosen to be an interpolation between the Alfvén and slow waves for simplicity. The numerical tests showed that the extended HLLD solver (called HLLD-P) has better performance for the capture of slow waves than the HLLD solver, and exhibits overall better accuracy in some situations where the slow waves exist. However, the new solver does not capture the Alfvén wave as well as the HLLD solver once the estimated speeds of compound waves deviate from the Alfvén wave speeds. Overall, the HLLD-P solver is fully compatible with the HLLD solver as long as the compound waves degenerate to the Alfvén waves inside the Riemann fan. It is indicated that the HLLD-P solver can be used for the various applications of MHD simulation, especially for those cases where the slow waves are expected to be generated.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1289-1302"},"PeriodicalIF":1.8,"publicationDate":"2025-06-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767609","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":"Comparison of Numerical Methods for Geometric Warpage Compensation","authors":"Steffen Tillmann,&nbsp;Stefan Basermann,&nbsp;Stefanie Elgeti","doi":"10.1002/fld.5404","DOIUrl":"https://doi.org/10.1002/fld.5404","url":null,"abstract":"<p>In injection molding processes, shrinkage and warpage cause deviations in the size and shape of produced parts compared to the cavity shape. While shrinkage is due to the change of material density during solidification, warpage is caused by uneven cooling and internal stresses within the part. One approach to mitigate these effects is by adjusting the cavity shape to anticipate the deformation. While finding the optimal cavity shape is often experience-based in practice, numerical design optimization can greatly assist in this process. In this study, we evaluate various numerical algorithms from existing literature to identify the optimal cavity shape. Each method is briefly outlined regarding how it adapts the geometry, and we discuss their respective strengths and weaknesses for different scenarios. We conduct comparisons using 3D geometries of varying complexity. Our findings demonstrate that, for geometric warpage compensation, the node-based reverse geometry method yields the least warpage and is computationally cost-effective. Furthermore, it is straightforward to implement and consistently performs well across different geometries.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1280-1288"},"PeriodicalIF":1.8,"publicationDate":"2025-05-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5404","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144768129","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":"Extrapolation Boundary Conditions for 2-D Smoothed Particle Hydrodynamics","authors":"Hossein Mahdizadeh,&nbsp;Colin D. Rennie,&nbsp;Benedict D. Rogers,&nbsp;Abolghasem Pilechi","doi":"10.1002/fld.5397","DOIUrl":"https://doi.org/10.1002/fld.5397","url":null,"abstract":"<p>This paper presents a new robust treatment for smoothed particle hydrodynamics (SPH) open (inflow/outflow) and solid boundary conditions, avoiding the unphysical fluctuations and numerical noise present in existing techniques. By novel use of concepts from finite volume methods, the fluid properties from sequential dynamic particles with different normal distances to the boundaries are extrapolated to ghost particles. No so-called mirror points are required, making the method computationally efficient and easy to implement. The new methodology is validated through a series of progressively challenging test cases. The effectiveness of the wall and inflow-outflow boundaries is evaluated for 2-D Poiseuille laminar flow. The performance of the wall boundary for complex geometries is demonstrated using a hydrostatic tank with a triangular wedge, followed by a conventional 2-D dam-break problem to capture impact pressures. A range of challenging vertical inflows rarely explored using SPH, with varying efflux velocities, demonstrate highly accurate performance of the boundary treatment, with results compared to STAR-CCM+. Finally, the robust performance is demonstrated for flow past circular and square cylinders over a range of Reynolds numbers, showing excellent results compared to reference results.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1248-1279"},"PeriodicalIF":1.8,"publicationDate":"2025-05-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5397","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767912","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":"Linear Discontinuity Sharpening for Highly Resolved and Robust Magnetohydrodynamics Simulations","authors":"Tomohiro Mamashita,&nbsp;Gaku Fukushima,&nbsp;Keiichi Kitamura","doi":"10.1002/fld.5402","DOIUrl":"https://doi.org/10.1002/fld.5402","url":null,"abstract":"<p>This study applies a reconstruction scheme, “hybrid MUSCL–THINC” for finite volume methods developed by Chiu et al., to magnetohydrodynamics (MHD) simulations. The scheme is a hybrid of monotone upstream-centered schemes for conservation law (MUSCL) and a tangent of hyperbola interface capturing (THINC) scheme. THINC sharply captures discontinuous distributions of physical quantities by using a hyperbolic tangent function. Our investigation reveals that hybrid MUSCL–THINC is more oscillatory in MHD simulations than in gas dynamics simulations, owing to the greater number of physical variables and associated complex waves in MHD. Analytical results demonstrate that artificial compression by THINC is excessive for MHD shock waves, whereas it is effective for linear discontinuities, such as contact discontinuities. Therefore, we propose a modification in which the artificial compression by THINC is weakened in the vicinity of nonlinear discontinuities and applied only to linear regions. The new scheme is tested using one- and two-dimensional MHD problems, and the results demonstrate that the scheme sharply captures linear discontinuities while avoiding numerical oscillations due to excessive artificial compression.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1226-1247"},"PeriodicalIF":1.8,"publicationDate":"2025-05-14","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5402","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767915","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":"A Framework Unifying Three-Cell-Based Scale-Invariant Exponential and Trigonometric WENO Weighting Functions With Optimal Shape Parameters","authors":"Xi Deng,&nbsp;Zhen-hua Jiang,&nbsp;Bin Xie,&nbsp;Chao Yan","doi":"10.1002/fld.5401","DOIUrl":"https://doi.org/10.1002/fld.5401","url":null,"abstract":"<div>\u0000 \u0000 <p>Exponential and trigonometric functions have been extensively employed as the kernel of reconstruction operators within numerous WENO (Weighted Essentially Non-oscillatory) schemes to accelerate the convergence rate. However, most of them are scale-dependent, compromising the robustness required for multi-scale flow simulations. Thus, this work aims to develop novel three-cell-based scale-invariant WENO schemes that use exponential and trigonometric functions as the kernel of non-linear weights. First, to achieve the scale-invariant property, this work reformulates the newly proposed scale-invariant ROUND (Reconstruction Operators in Unified Normalized-variable Diagram) schemes into the form of WENO weighting functions, thereby facilitating the design of scale-invariant WENO schemes. Then, this work proposes new WENO non-linear weights using exponential and trigonometric functions—such as Gaussian, hyperbolic, and cosine functions—to enhance the accuracy of the three-cell-based WENO scheme. The proposed WENO weights contain a shape parameter that controls the errors between the non-linear weight and the ideal weight. As the value of the shape parameter increases, the non-linear weight converges towards the ideal weight but also becomes more likely to produce numerical oscillations. To approximate the optimal value of the shape parameter, the WENO reconstruction operator is projected into normalized variable space, and the shape parameters are fine-tuned to ensure the normalized reconstruction operator falls into the CBC (Convection Bounded Criterion) region of UND (Unified Normalized-variable Diagram). The accuracy analysis reveals that the proposed weighting functions outperform classical WENO schemes, particularly when the smooth function contains a first-order critical point. The accuracy and shock-capturing properties of the proposed schemes are further validated through benchmark tests. Thus, this work demonstrates using the ROUND framework to design scale-invariant three-cell-based WENO schemes with exponential and trigonometric functions and optimal shape parameters.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"97 9","pages":"1209-1225"},"PeriodicalIF":1.8,"publicationDate":"2025-05-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"144767490","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":"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}