The numerical study of flow around a pair of spheres and a square array of spheres is investigated by using a direct‐forcing immersed boundary method. Using high resolution three‐dimensional computations, we analyzed the flow around several configurations: a sphere, a pair of spheres in a tandem arrangement with center‐to‐center streamwise ratio L/D ranging from 1 to 6, and a square array with 9 spheres in a uniform arrangement. In the latter case, we explore the ratio of array diameter (DG) to sphere diameter (D) at 4, 5, 6 and 7. The center‐to‐center streamwise and transverse pitch is the same, varied from L/D = 1.5, 2, 2.5 to 3, and they were arranged in a square periodic array to allow uniform distribution within the array. Based on the effective direct‐forcing immersed boundary projection method, the fractional time marching methodology is applied for solving four field variables involving three velocities and one pressure component. The pressure Poisson equation is advanced in space by using the fast Fourier transform (FFT) and a tridiagonal matrix algorithm (TDMA), effectively solving for the diagonally dominant tridiagonal matrix equations. A direct‐forcing immersed boundary method is involved to treat the interfacial terms by adding the appropriate sources as force function at the boundary, separating the phases. Geometries featuring the stationary solid obstacles in the flow are embedded in the Cartesian grid with special discretizations near the embedded boundary using a discrete Dirac delta function to ensure the accuracy of the solution in the cut cells. An important characteristic of flow over the multiple spheres is devised by comparing with the drag and lift coefficients, as well as vortex shedding.
{"title":"An efficient direct‐forcing immersed boundary method for flow around a pair of spheres","authors":"Der Chang Lo, Katherine Lee, Pao‐Lan Shen","doi":"10.1002/fld.5326","DOIUrl":"https://doi.org/10.1002/fld.5326","url":null,"abstract":"The numerical study of flow around a pair of spheres and a square array of spheres is investigated by using a direct‐forcing immersed boundary method. Using high resolution three‐dimensional computations, we analyzed the flow around several configurations: a sphere, a pair of spheres in a tandem arrangement with center‐to‐center streamwise ratio L/D ranging from 1 to 6, and a square array with 9 spheres in a uniform arrangement. In the latter case, we explore the ratio of array diameter (<jats:italic>D</jats:italic><jats:sub>G</jats:sub>) to sphere diameter (D) at 4, 5, 6 and 7. The center‐to‐center streamwise and transverse pitch is the same, varied from L/D = 1.5, 2, 2.5 to 3, and they were arranged in a square periodic array to allow uniform distribution within the array. Based on the effective direct‐forcing immersed boundary projection method, the fractional time marching methodology is applied for solving four field variables involving three velocities and one pressure component. The pressure Poisson equation is advanced in space by using the fast Fourier transform (FFT) and a tridiagonal matrix algorithm (TDMA), effectively solving for the diagonally dominant tridiagonal matrix equations. A direct‐forcing immersed boundary method is involved to treat the interfacial terms by adding the appropriate sources as force function at the boundary, separating the phases. Geometries featuring the stationary solid obstacles in the flow are embedded in the Cartesian grid with special discretizations near the embedded boundary using a discrete Dirac delta function to ensure the accuracy of the solution in the cut cells. An important characteristic of flow over the multiple spheres is devised by comparing with the drag and lift coefficients, as well as vortex shedding.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141866391","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}
Granular flow problems characterized by large deformations are widespread in various applications, including coastal and geotechnical engineering. The paper deals with the application of a rigid‐perfectly plastic two‐phase model extended by the Drucker–Prager yield criterion to simulate granular media with a finite volume flow solver (FV). The model refers to the combination of a Bingham fluid and an Eulerian strain measure to assess the failure region of granular dam slides. A monolithic volume‐of‐fluid (VoF) method is used to distinguish between the air and granular phases, both governed by the incompressible Navier–Stokes equations. The numerical framework enables modeling of large displacements and arbitrary shapes for large‐scale applications. The displayed validation and verification focuses on the rigid‐perfectly plastic material model for noncohesive and cohesive materials with varying angles of repose. Results indicate a good agreement of the predicted soil surface and strain results with experimental and numerical data.
{"title":"A two‐phase volume of fluid approach to model rigid‐perfectly plastic granular materials","authors":"W. Düsterhöft‐Wriggers, S. Schubert, T. Rung","doi":"10.1002/fld.5323","DOIUrl":"https://doi.org/10.1002/fld.5323","url":null,"abstract":"Granular flow problems characterized by large deformations are widespread in various applications, including coastal and geotechnical engineering. The paper deals with the application of a rigid‐perfectly plastic two‐phase model extended by the Drucker–Prager yield criterion to simulate granular media with a finite volume flow solver (FV). The model refers to the combination of a Bingham fluid and an Eulerian strain measure to assess the failure region of granular dam slides. A monolithic volume‐of‐fluid (VoF) method is used to distinguish between the air and granular phases, both governed by the incompressible Navier–Stokes equations. The numerical framework enables modeling of large displacements and arbitrary shapes for large‐scale applications. The displayed validation and verification focuses on the rigid‐perfectly plastic material model for noncohesive and cohesive materials with varying angles of repose. Results indicate a good agreement of the predicted soil surface and strain results with experimental and numerical data.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141866463","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}
The classical staggered 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.
经典交错谱元法(SEM)被重新审视并扩展到准静态磁流体动力学(MHD)流。在磁雷诺数消失的极限条件下,动量方程中洛伦兹力的评估需要确定受欧姆定律和安培定则电荷守恒条件支配的电流密度。一旦使用 SEM 进行离散化,这就转化为解决涉及所谓一致泊松算子的电动势的额外问题。该方法非常适合复杂几何结构中的全三维流动。除了分辨率要求的变化之外,考虑电磁量估计会使 MHD 的相关计算成本比流体力学增加约 40%。该方案的准确性和功能已在 MHD 文献中的一组常见流体上得到验证。电流密度的多项式阶指数收敛得到了证实。
{"title":"A spectral element discretization for quasi‐static magnetohydrodynamic flows","authors":"Mattias Brynjell‐Rahkola","doi":"10.1002/fld.5321","DOIUrl":"https://doi.org/10.1002/fld.5321","url":null,"abstract":"The classical staggered 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.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141744841","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}
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.
{"title":"Large Courant–Friedrichs–Lewy explicit scheme for one‐dimensional hyperbolic conservation laws","authors":"Vincent Guinot, Antoine Rousseau","doi":"10.1002/fld.5322","DOIUrl":"https://doi.org/10.1002/fld.5322","url":null,"abstract":"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.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-05","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}
SummaryIn 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.
{"title":"A hybrid marching cubes based IsoAlpha method for interface reconstruction","authors":"G.S. Abhishek, Shyamprasad Karagadde","doi":"10.1002/fld.5320","DOIUrl":"https://doi.org/10.1002/fld.5320","url":null,"abstract":"SummaryIn 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.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-04","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}
Marina G. Kontou, Xenofon S. Trompoukis, Varvara G. Asouti, Kyriakos C. Giannakoglou
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 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.
{"title":"The continuous adjoint method to the γ−R˜eθt$$ gamma -tilde{R}{e}_{theta t} $$ transition model coupled with the Spalart–Allmaras model for compressible flows","authors":"Marina G. Kontou, Xenofon S. Trompoukis, Varvara G. Asouti, Kyriakos C. Giannakoglou","doi":"10.1002/fld.5319","DOIUrl":"https://doi.org/10.1002/fld.5319","url":null,"abstract":"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 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.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-07-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"141503919","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}
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.
{"title":"A modified lattice Boltzmann approach based on radial basis function approximation for the non‐uniform rectangular mesh","authors":"X. Hu, J. M. Bergadà, D. Li, W. M. Sang, B. An","doi":"10.1002/fld.5318","DOIUrl":"https://doi.org/10.1002/fld.5318","url":null,"abstract":"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.","PeriodicalId":50348,"journal":{"name":"International Journal for Numerical Methods in Fluids","volume":null,"pages":null},"PeriodicalIF":1.8,"publicationDate":"2024-06-25","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}
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