International Journal for Numerical Methods in Fluids最新文献

{"title":"Robust Linearly-Implicit Backward Difference Formulas for Navier–Stokes Equations: Computations of Steady and Unsteady Flows","authors":"Kak Choon Loy, Yves Bourgault","doi":"10.1002/fld.70051","DOIUrl":"https://doi.org/10.1002/fld.70051","url":null,"abstract":"<div>\u0000 \u0000 <p>We propose a linearly-implicit method (called LBDFT) to solve the incompressible Navier–Stokes equations. Linearly-implicit methods have an algorithmic complexity that lies between fully-implicit and semi-implicit time-stepping schemes. In LBDFT, the nonlinear advection in the Navier–Stokes equations is split into three linear terms using a Taylor series expansion. One term is taken explicitly, and the other two are updated with the linear diffusion term at each time step. An additional variant of linearly-implicit methods, referred to as LBDFE, was also included in this study. It is based on an extrapolation formula similar to those used in semi-implicit methods. For the sake of comparison, we recall a set of fully-implicit and semi-implicit time-stepping methods, for the most part based on backward differentiation formulae (BDF). These methods are compared in terms of accuracy, stability, computing time, and ability to compute various flows. We first use a two-dimensional manufactured problem to assess the order convergence (in time) of the methods. The second test case is a two-dimensional unsteady flow around a cylinder at <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>R</mi>\u0000 <mi>e</mi>\u0000 <mo>=</mo>\u0000 <mn>100</mn>\u0000 </mrow>\u0000 <annotation>$$ mathit{operatorname{Re}}=100 $$</annotation>\u0000 </semantics></math>. We observed that linearly-implicit methods are more CPU efficient compared to fully-implicit BDF, both at second- and third-order accuracy. Our third test case explores the ability of the methods to compute steady flows at high Reynolds numbers, in our case a steady two-dimensional lid-driven cavity for <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>R</mi>\u0000 <mi>e</mi>\u0000 <mo>=</mo>\u0000 <mn>1,000</mn>\u0000 </mrow>\u0000 <annotation>$$ mathit{operatorname{Re}}=mathrm{1,000} $$</annotation>\u0000 </semantics></math> to <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mn>65</mn>\u0000 <mo>,</mo>\u0000 <mn>000</mn>\u0000 </mrow>\u0000 <annotation>$$ 65,000 $$</annotation>\u0000 </semantics></math>. Our LBDFT method was the most efficient for capturing these flows, with solutions almost identical to published results, while LBDFE methods do not work at all. LBDFT was able to compute steady cavity flows for <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>R</mi>\u0000 <mi>e</mi>\u0000 <mo>=</mo>\u0000 <mn>100</mn>\u0000 <mo>,</mo>\u0000 <mn>000</mn>\u0000 </mrow>\u0000 <annotation>$$ mathit{operatorname{Re}}=100,000 $$</annotation>\u0000 </semantics></math> to <span></span><math>\u0000 ","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 4","pages":"510-528"},"PeriodicalIF":1.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147562452","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":"Optimal Feedback Force Locations in Measurement-Integrated Simulation for Two-Dimensional Flow Around a Circular Cylinder","authors":"Kota Hirose,&nbsp;Suguru Miyauchi","doi":"10.1002/fld.70045","DOIUrl":"https://doi.org/10.1002/fld.70045","url":null,"abstract":"<div>\u0000 \u0000 <p>Measurement-integrated simulation, one of the data assimilation methods, is a simulation technique that reproduces an actual flow field by adding feedback forces, including measurement data, to the external force term of the Navier–Stokes equations. The analysis in the previous studies using wind tunnels was limited to flows around a square cylinder. There is no principle for determining the optimal location of the feedback force, and the optimal location of feedback forces for other shapes is unknown. This study aimed to identify the optimal location for applying the feedback force in the flow around a circular cylinder and explore its relationship with the flow field. Using pressure on the cylinder surface obtained from a numerical experiment with disturbed inflow containing random fluctuations as measurement data, several measurement-integrated simulations were conducted, each with a different feedback location. By comparing the results of the measurement-integrated simulation with those of the numerical experiment, the optimal position was identified at 120° from the stagnation point that most accurately reproduced the flow behind the cylinder in the numerical experiment. Furthermore, this location was identified as the flow separation point, suggesting that applying a feedback force at the separation point is optimal for reproducing the flow behind the object.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 4","pages":"349-359"},"PeriodicalIF":1.8,"publicationDate":"2026-03-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"147564471","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":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 6","pages":"766-782"},"PeriodicalIF":1.8,"publicationDate":"2026-02-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148066170","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":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 6","pages":"751-765"},"PeriodicalIF":1.8,"publicationDate":"2026-01-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148098742","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":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 6","pages":"709-730"},"PeriodicalIF":1.8,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148098568","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":"","authors":"","doi":"","DOIUrl":"","url":null,"abstract":"","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 6","pages":"731-750"},"PeriodicalIF":1.8,"publicationDate":"2026-01-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"148098569","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":"Analysis of Turbulent Modeling for Free-Surface Flows Using a Hybrid RANS-LES Model and Particle-Based Moving Particle Semi-Implicit Method","authors":"Fabio Kenji Motezuki,&nbsp;Lucas Soares Pereira,&nbsp;Liang-Yee Cheng,&nbsp;Fernando Akira Kurokawa","doi":"10.1002/fld.70031","DOIUrl":"10.1002/fld.70031","url":null,"abstract":"&lt;p&gt;Engineering problems often comprise free-surface flows in turbulent regime. Lagrangian mesh-free particle-based methods are well suited for the simulation of flows involving complex free-surface deformation. However, the analysis of turbulent modeling for particle-based methods is relatively scarce in the literature. In this work, an analysis of a hybrid RANS-LES turbulence model adapted for the Moving Particle Semi-implicit (MPS) method is performed. In the turbulence model, a zero-equation RANS is applied near the wall boundaries and a standard Smagorinsky LES model is applied elsewhere. Given that the eddy viscosity of the turbulent modeling depends on the distance between the fluid and the nearest wall particle, the calculation of the fluid-wall particle distance may demand a high computational cost due to undefined topology among moving particles. In this way, a method based on the cell-linked list is proposed to improve the nearest wall search for the turbulence model. The implementation is verified through simulation of a lid-driven flow with Reynolds number between &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mn&gt;10,000&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ mathrm{10,000} $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mn&gt;50,000&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ mathrm{50,000} $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;. The result shows that despite the overhead when the turbulence model is adopted, the time needed to reach steady state is shortened so that the overall computational costs are almost the same. In addition, the improvement due to the adoption of turbulence model is more evident for the highest Reynolds numbers. As an application, the flow around a submerged square cylinder near the surface with Reynolds number of &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mn&gt;25,000&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ mathrm{25,000} $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; is simulated. The influences of the cylinder submergence depths on the drag and lift coefficients are investigated for a range of depth-to-length ratios between &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mn&gt;0.3&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ 0.3 $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt; and &lt;span&gt;&lt;/span&gt;&lt;math&gt;\u0000 &lt;semantics&gt;\u0000 &lt;mrow&gt;\u0000 &lt;mn&gt;3.0&lt;/mn&gt;\u0000 &lt;/mrow&gt;\u0000 &lt;annotation&gt;$$ 3.0 $$&lt;/annotation&gt;\u0000 &lt;/semantics&gt;&lt;/math&gt;. When the turbulence model is applied, a smoother convergence tendency is obtained as the resolution increases. Moreover, the flow around the square cylinder is better represented","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 3","pages":"321-347"},"PeriodicalIF":1.8,"publicationDate":"2025-12-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.70031","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146193322","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 Mathematical Model for Two-Phase Flow in Confined Environments: Numerical Solution and Validation","authors":"Giuseppe Sciumè,&nbsp;Haohong Pi,&nbsp;Abdelaziz Omari,&nbsp;Thomas Lavigne","doi":"10.1002/fld.70030","DOIUrl":"10.1002/fld.70030","url":null,"abstract":"<p>This study presents a numerical framework for modeling two-phase flow in confined environments, focusing on the interplay between capillary and viscous forces. The model integrates the Cahn-Hilliard and Navier-Stokes (CH-NS) equations, utilizing a diffuse-interface approach to capture interfacial dynamics without the limitations of sharp-interface models. Implemented in the finite element platform <i>FEniCS</i>, the framework incorporates Dirichlet boundary conditions to model a fully non-wetting phase. The validation of the proposed model is achieved through two applications: The retraction of an oil droplet from a capillary tube and the drainage of water-wet microfluidic chips. Numerical results align with experimental data, demonstrating the framework's ability to replicate interfacial behaviors, including capillary-driven dynamics and fingering phenomena. This work provides a versatile computational tool for studying immiscible fluid flow, offering potential for advancements in fundamental research on microfluidics, enhanced oil recovery, and remediation of contaminated soil.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 3","pages":"306-320"},"PeriodicalIF":1.8,"publicationDate":"2025-12-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.70030","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146193449","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":"Buoyant Magneto-Convection in an Internally Heated Anisotropic Porous Cavity With Sinusoidal Boundary Flux","authors":"Manisha Jangir,&nbsp;Poosan Muthu,&nbsp;Jaya Krishna Devanuri","doi":"10.1002/fld.70028","DOIUrl":"10.1002/fld.70028","url":null,"abstract":"<div>\u0000 \u0000 <p>This study examines buoyancy-driven magneto-convection within an anisotropic porous cavity incorporating internal heat generation/absorption. The top and bottom boundaries are subjected to sinusoidal heat fluxes, whereas the vertical walls are thermally insulated. The flow and heat transfer behavior are numerically analyzed using the Darcy-Brinkman extended model, implemented via FVM and the semi-implicit method for pressure-linked equations (SIMPLE) algorithm. The influence of key parameters, including the periodicity parameter, permeability ratio, thermal conductivity ratio, Hartmann number, internal heat generation/absorption, and orientation angle on the flow structure and heat transfer efficiency within the system is analyzed. The findings show that the cavity exhibits a multicellular convective pattern, where the anisotropic permeability tilt induces sinusoidal flow features near the thermally active walls. In contrast, a strong magnetic field (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>H</mi>\u0000 <mi>a</mi>\u0000 <mo>=</mo>\u0000 <mn>100</mn>\u0000 </mrow>\u0000 <annotation>$$ Ha=100 $$</annotation>\u0000 </semantics></math>) suppresses the flow circulation, and higher levels of internal heat generation/absorption <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mo>(</mo>\u0000 <mi>Q</mi>\u0000 <mo>=</mo>\u0000 <mn>5</mn>\u0000 <mo>,</mo>\u0000 <mo>−</mo>\u0000 <mn>5</mn>\u0000 <mo>)</mo>\u0000 </mrow>\u0000 <annotation>$$ left(Q=5,-5right) $$</annotation>\u0000 </semantics></math> lead to reduced heat transfer efficiency. The study is conducted for a steady, two-dimensional, Darcy-Brinkman model. This study could be beneficial for the solar collector designs, thermal management systems, setups for room ventilation, and electronic cooling applications.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 3","pages":"290-305"},"PeriodicalIF":1.8,"publicationDate":"2025-12-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146193351","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":"Positivity-Preserving and Shock-Capturing via Adaptive Filtering Direct Flux Reconstruction Scheme for Solving Multi-Component Euler Equations","authors":"Raagvendra Singh,&nbsp;Abhishek M. Kalluri,&nbsp;V. K. Suman,&nbsp;Rakesh Kumar","doi":"10.1002/fld.70029","DOIUrl":"10.1002/fld.70029","url":null,"abstract":"<div>\u0000 \u0000 <p>An extension of the direct flux reconstruction scheme for simulating compressible multi-component flows using the Euler equations is presented via adaptive filtering and positivity preservation. It is not an unknown fact that high-order schemes are not robust in the vicinity of shocks and other flow discontinuities. Thus, in the presented work, an adaptive filtering procedure for shock capturing, which dissipates Gibbs oscillations around shocks, has been used. While shock capturing dissipates spurious oscillations associated with shocks, it doesn't guarantee that the physical quantities will not attain unphysical states. These unphysical states are not just due to high-order interpolations, but also when these physical values are already near zero. Thus, in the present work, we also use a simple positivity-preserving scheme that can be applied when solving the multi-component Euler equations. The developed flow solver is then validated using problems for both 1-dimensional and 2-dimensional flows.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"98 3","pages":"274-289"},"PeriodicalIF":1.8,"publicationDate":"2025-11-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"146193669","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}