International Journal for Numerical Methods in Fluids最新文献

{"title":"A spectral element discretization for quasi-static magnetohydrodynamic flows","authors":"Mattias Brynjell-Rahkola","doi":"10.1002/fld.5321","DOIUrl":"10.1002/fld.5321","url":null,"abstract":"<p>The classical staggered <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ℙ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 </mrow>\u0000 </msub>\u0000 <mi>-</mi>\u0000 <msub>\u0000 <mrow>\u0000 <mi>ℙ</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>N</mi>\u0000 <mo>−</mo>\u0000 <mn>2</mn>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ {mathbb{P}}_Nhbox{-} {mathbb{P}}_{N-2} $$</annotation>\u0000 </semantics></math> spectral element method (SEM) is revisited and extended to quasi-static magnetohydrodynamic (MHD) flows. In this realm, which is valid in the limit of vanishing magnetic Reynolds number, the evaluation of the Lorentz force in the momentum equation requires the electric current density, governed by Ohm's law and a charge conservation condition derived from Ampère's law, to be determined. Once discretized with the SEM, this translates into solving one additional problem for the electric potential involving the so-called consistent Poisson operator. The method is well suited for fully three-dimensional flows in complex geometries. Changes in resolution requirements aside, consideration of the electromagnetic quantities is estimated to increase the computational cost associated with MHD by about 40% relative to hydrodynamics. The accuracy and the capabilities of the scheme is demonstrated on a set of common flows from the MHD literature. Exponential convergence with polynomial order is confirmed for the electric current density.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1795-1812"},"PeriodicalIF":1.7,"publicationDate":"2024-07-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5321","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141744841","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":"Large Courant–Friedrichs–Lewy explicit scheme for one-dimensional hyperbolic conservation laws","authors":"Vincent Guinot,&nbsp;Antoine Rousseau","doi":"10.1002/fld.5322","DOIUrl":"10.1002/fld.5322","url":null,"abstract":"<p>A large Courant–Friedrichs–Lewy (CFL) algorithm is presented for the explicit, finite volume solution of hyperbolic systems of conservation laws, with a focus on the shallow water equations. The Riemann problems used in the flux computation are determined using averaging kernels that extend over several computational cells. The usual CFL stability constraint is replaced with a constraint involving the kernel support size. This makes the method unconditionally stable with respect to the size of the computational cells, allowing the computational mesh to be refined locally to an arbitrary degree without altering solution stability. The practical implementation of the method is detailed for the shallow water equations with topographical source term. Computational examples report applications of the method to the linear advection, Burgers and shallow water equations. In the case of sharp bottom discontinuities, the need for improved, well-balanced discretisations of the geometric source term is acknowledged.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1760-1794"},"PeriodicalIF":1.7,"publicationDate":"2024-07-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141578024","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 hybrid marching cubes based IsoAlpha method for interface reconstruction","authors":"G.S. Abhishek,&nbsp;Shyamprasad Karagadde","doi":"10.1002/fld.5320","DOIUrl":"10.1002/fld.5320","url":null,"abstract":"<div>\u0000 \u0000 <p>In modelling two-phase flows, accurate representation of interfaces is crucial. A class of methods for interface reconstruction are based on isosurface extraction, which involves a non-iterative, interpolation based approach. These approaches have been shown to be faster by an order of magnitude than the conventional PLIC schemes. In this work, we present a new isosurface extraction based interface reconstruction scheme based on the Marching Cubes algorithm (MC), which is commonly used in computer graphics for visualizing isosurfaces. The MC algorithm apriori lists and categorizes all possible interface configurations in a single grid cell into a Look Up Table (LUT), which makes this approach fast and robust. We also show that for certain interface configurations, the inverse problem of obtaining the isovalue from the cell volume fraction is not surjective, and a special treatment is required while handling these cases. We then demonstrate the capabilities of the method through benchmark cases for 2D and 3D static/dynamic interface reconstruction.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1741-1759"},"PeriodicalIF":1.7,"publicationDate":"2024-07-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141550771","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":"The continuous adjoint method to the \u0000 \u0000 \u0000 γ\u0000 \u0000 −\u0000 \u0000 \u0000 \u0000 R\u0000 \u0000 ˜\u0000 \u0000 \u0000 \u0000 e\u0000 \u0000 \u0000 θ\u0000 t\u0000 \u0000 \u0000 \u0000 $$ gamma -tilde{R}{e}_{theta t} $$\u0000 transition model coupled with the Spalart–Allmaras model for compressible flows","authors":"Marina G. Kontou,&nbsp;Xenofon S. Trompoukis,&nbsp;Varvara G. Asouti,&nbsp;Kyriakos C. Giannakoglou","doi":"10.1002/fld.5319","DOIUrl":"10.1002/fld.5319","url":null,"abstract":"<p>The continuous adjoint method for transitional flows of compressible fluids is developed and assessed, for the first time in the literature. The gradient of aerodynamic objective functions (aerodynamic forces) with respect to design variables, in problems governed by the compressible Navier–Stokes equations coupled with the Spalart–Allmaras turbulence model and the <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>γ</mi>\u0000 <mspace></mspace>\u0000 <mo>−</mo>\u0000 <mspace></mspace>\u0000 <mover>\u0000 <mrow>\u0000 <mi>R</mi>\u0000 </mrow>\u0000 <mo>˜</mo>\u0000 </mover>\u0000 <msub>\u0000 <mrow>\u0000 <mi>e</mi>\u0000 </mrow>\u0000 <mrow>\u0000 <mi>θ</mi>\u0000 <mi>t</mi>\u0000 </mrow>\u0000 </msub>\u0000 </mrow>\u0000 <annotation>$$ gamma -tilde{R}{e}_{theta t} $$</annotation>\u0000 </semantics></math> transition model (in three, non-smooth and smooth, variants of it), is computed based on the continuous adjoint method. The development of the adjoint to the smooth transition model variant proved to be beneficial. The accuracy of the computed sensitivity derivatives is verified against finite differences. Programming is performed in an in-house, vertex-centered finite-volume code, efficiently running on GPUs. The proposed continuous adjoint method is used in 2D and 3D aerodynamic shape optimization problems, namely the constrained optimization of the NLF(1)–0416 isolated airfoil and that of the ONERA M6 wing. The impact of “frozen transition” (assumption according to which the adjoint to the transition model equations are not solved) or “frozen turbulence” (by additionally ignoring the adjoint to the turbulence model) are evaluated; it is shown that both lead to inaccurate sensitivities.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1715-1740"},"PeriodicalIF":1.7,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1002/fld.5319","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503919","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 modified lattice Boltzmann approach based on radial basis function approximation for the non-uniform rectangular mesh","authors":"X. Hu,&nbsp;J. M. Bergadà,&nbsp;D. Li,&nbsp;W. M. Sang,&nbsp;B. An","doi":"10.1002/fld.5318","DOIUrl":"10.1002/fld.5318","url":null,"abstract":"<p>We have presented a novel lattice Boltzmann approach for the non-uniform rectangular mesh based on the radial basis function approximation (RBF-LBM). The non-uniform rectangular mesh is a good option for local grid refinement, especially for the wall boundaries and flow areas with intensive change of flow quantities. Which allows, the total number of grid cells to be reduced and so the computational cost, therefore improving the computational efficiency. But the grid structure of the non-uniform rectangular mesh is no longer applicable to the classic lattice Boltzmann method (CLBM), which is based on the famous BGK collision-streaming evolution. This is why the present study is inspired by the idea of the interpolation-supplemented LBM (ISLBM) methodology. The ISLBM algorithm is improved in the present manuscript and developed into a novel LBM approach through the radial basis function approximation instead of the Lagrangian interpolation scheme. The new approach is validated for both steady states and unsteady periodic solutions. The comparison between the radial basis function approximation and the Lagrangian interpolation is discussed. It is found that the novel approach has a good performance on computational accuracy and efficiency. Proving that the non-uniform rectangular mesh allows grid refinement while obtaining precise flow predictions.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 11","pages":"1695-1714"},"PeriodicalIF":1.7,"publicationDate":"2024-06-24","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503920","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":"Decoupled and unconditionally stable iteration method for stationary Navier–Stokes equations","authors":"Jianhua Chen,&nbsp;Yingying Jiang,&nbsp;Guo-Dong Zhang","doi":"10.1002/fld.5317","DOIUrl":"https://doi.org/10.1002/fld.5317","url":null,"abstract":"<p>It is well known the Oseen iteration for the stationary Navier–Stokes equations is unconditionally stable. However, it is a coupled type scheme where the velocity <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>u</mi>\u0000 </mrow>\u0000 <annotation>$$ boldsymbol{u} $$</annotation>\u0000 </semantics></math> and pressure <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>p</mi>\u0000 </mrow>\u0000 <annotation>$$ p $$</annotation>\u0000 </semantics></math> are coupled together at each iteration. By treating pressure <span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>p</mi>\u0000 </mrow>\u0000 <annotation>$$ p $$</annotation>\u0000 </semantics></math> explicitly would lead to a decoupled iteration, but this treatment is unstable. In this article, we construct a decoupled and unconditionally stable iteration method to solve the stationary Navier–Stokes equations by adopting the pressure projection method to the temporal disturbed Navier–Stokes system whose solution approximates the steady state solution over time (<span></span><math>\u0000 <semantics>\u0000 <mrow>\u0000 <mi>t</mi>\u0000 <mo>→</mo>\u0000 <mo>+</mo>\u0000 <mi>∞</mi>\u0000 </mrow>\u0000 <annotation>$$ tto +infty $$</annotation>\u0000 </semantics></math>). We also rigorously prove its unconditional stability. Numerical simulations demonstrate that our iterative method is more efficient and stable than the extensively used T-S and Oseen iterations, and could solve the fluid flow with high Reynolds number.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 10","pages":"1680-1693"},"PeriodicalIF":1.7,"publicationDate":"2024-06-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142137867","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 general pressure equation based method for incompressible two-phase flows","authors":"Hormuzd Bodhanwalla,&nbsp;Dheeraj Raghunathan,&nbsp;Y. Sudhakar","doi":"10.1002/fld.5316","DOIUrl":"https://doi.org/10.1002/fld.5316","url":null,"abstract":"<div>\u0000 \u0000 <p>We present a fully-explicit, iteration-free, weakly-compressible method to simulate immiscible incompressible two-phase flows. To update pressure, we circumvent the computationally expensive Poisson equation and use the general pressure equation which is solved explicitly. In addition, a less diffusive algebraic volume-of-fluid approach is used as the interface capturing technique and in order to facilitate improved parallel computing scalability, the technique is discretised temporally using the operator-split methodology. Our method is fully-explicit and stable with simple local spatial discretization, and hence, it is easy to implement. Several two- and three-dimensional canonical two-phase flows are simulated. The qualitative and quantitative results prove that our method is capable of accurately handling problems involving a range of density and viscosity ratios and surface tension effects.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 10","pages":"1653-1679"},"PeriodicalIF":1.7,"publicationDate":"2024-06-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"142130384","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 modified fifth-order WENO-Z scheme based on the weights of the reformulated adaptive order WENO scheme","authors":"Yize Wang,&nbsp;Kunlei Zhao,&nbsp;Li Yuan","doi":"10.1002/fld.5314","DOIUrl":"10.1002/fld.5314","url":null,"abstract":"<p>A modified fifth-order WENO-Z scheme is proposed by analogy with the non-normalized weights of the reformulated fifth-order adaptive order (AO) WENO scheme. We show that if the original fifth-order WENO-AO scheme is rewritten as the form of the conventional WENO combination, the resulting non-normalized weights can be divided into three parts: a constant one term, a local stencil smoothness measure term and a global stencil smoothness measure term. In order to make use of the latter two terms for constructing a modified WENO-Z scheme with enhanced performance, we change the form of the third term and introduce an adaptive scaling factor to adjust the contributions from the second and third terms. Numerical examples show that the modified fifth-order WENO-Z scheme has the advantage of high resolution in smooth regions and sharp capturing of discontinuities, and it can obtain evidently better results for shocked flows with small-scale structures compared with the recently developed WENO-Z+ and WENO-Z+M schemes.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 10","pages":"1631-1652"},"PeriodicalIF":1.7,"publicationDate":"2024-06-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141387648","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":"Research on the construction of multi objective coupling model and optimization method of ship form","authors":"Jie Liu,&nbsp;Baoji Zhang,&nbsp;Yuyang Lai,&nbsp;Liqiao Fang","doi":"10.1002/fld.5315","DOIUrl":"10.1002/fld.5315","url":null,"abstract":"<p>Multi-objective optimization of ship form can effectively reduce ship energy consumption, and is one of the important research topics of green ships. However, the computational cost of numerical simulation based on computational fluid dynamics (CFD) theory is relatively high, which affects the efficiency of optimization. Traditional subjective weighting methods mostly rely on expert's experience, which affects the scientificity of optimization. This paper effectively integrates the CFD method, the improved multi-objective optimization algorithm and the objective weighting method to build a ship form multi-objective optimization framework. Conduct multi-objective optimization research on resistance and seakeeping performance of a very large crude oil carrier (KVLCC) ship. The improved bare-bones multi-objective particle swarm optimization (IBBMOPSO) algorithm is used to obtain the pareto front, and the kernel principal component analysis (KPCA) method is used to objectively assign the weight of each target. Finally, the optimal ship form scheme with high satisfaction was obtained. The multi-objective optimization framework constructed in this paper can provide a certain theoretical basis and technical support for the development of ship greening and digital transformation.</p>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 10","pages":"1617-1630"},"PeriodicalIF":1.7,"publicationDate":"2024-06-03","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141271050","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":"Simplified weighting formulations of weighted compact nonlinear schemes for compressible flows","authors":"Jinwei Bai,&nbsp;Meiliang Mao,&nbsp;Yankai Ma,&nbsp;Zhen-Guo Yan,&nbsp;Yaobing Min","doi":"10.1002/fld.5311","DOIUrl":"10.1002/fld.5311","url":null,"abstract":"<div>\u0000 \u0000 <p>Weighted compact nonlinear schemes (WCNSs) are a popular family of high-resolution shock-capturing schemes for simulating compressible flows, of which the nonlinear interpolation procedure is dominant for the performance. In this work, a simplified weighting strategy is introduced for the nonlinear interpolation procedure. Firstly, an equivalent weighting formulation of WCNS is presented by explicitly including the whole-point stencil into the set of candidate stencils. Secondly, motivated by the reorganization of WCNS, the WCNS-CU6 scheme is achieved in a more straightforward way. Thirdly, by introducing a TENO selection procedure in the framework of WCNS-CU6-Simplified, a TCNS6-Simplified scheme is proposed, the resolution of which is comparable with the excellent TENO6 scheme, while the computational cost is much lower. The simplified schemes exhibit more outstanding, at least comparable, fidelity than the original schemes, however, with superior characteristics in terms of efficiency and simplicity. A variety of benchmark test problems are studied to demonstrate the behaviour of the simplified weighting strategy.</p>\u0000 </div>","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":"96 9","pages":"1592-1616"},"PeriodicalIF":1.7,"publicationDate":"2024-05-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141189985","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}