首页 > 最新文献

International Journal for Numerical Methods in Fluids最新文献

英文 中文
An efficient direct‐forcing immersed boundary method for flow around a pair of spheres 绕一对球体流动的高效直接强迫沉浸边界法
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-31 DOI: 10.1002/fld.5326
Der Chang Lo, Katherine Lee, Pao‐Lan Shen
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.
我们采用直接强迫沉浸边界法对一对球体和一个方形球体阵列周围的流动进行了数值研究。通过高分辨率三维计算,我们分析了几种结构周围的流动:一个球体、一对串联排列的球体(中心到中心的流向比 L/D 从 1 到 6 不等)以及一个由 9 个球体均匀排列的正方形阵列。在后一种情况下,我们探讨了阵列直径 (DG) 与球体直径 (D) 在 4、5、6 和 7 之间的比率。中心到中心的流向和横向间距是相同的,从 L/D = 1.5、2、2.5 到 3 不等,它们排列成正方形周期阵列,以便在阵列内均匀分布。基于有效的直接强迫沉浸边界投影法,应用分数时间行进法求解涉及三个速度和一个压力分量的四个场变量。通过使用快速傅立叶变换(FFT)和三对角矩阵算法(TDMA)在空间推进压力泊松方程,有效求解对角主导三对角矩阵方程。直接强迫沉浸边界法通过在边界添加适当的源作为力函数来处理界面项,从而分离相位。以流动中的静止固体障碍物为特征的几何图形被嵌入笛卡尔网格中,嵌入边界附近使用离散 Dirac delta 函数进行特殊离散,以确保切割单元中求解的准确性。通过与阻力和升力系数以及涡流脱落进行比较,设计出了多个球体上流动的一个重要特征。
{"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}
引用次数: 0
A two‐phase volume of fluid approach to model rigid‐perfectly plastic granular materials 用两相流体体积法建立刚-完全塑性颗粒材料模型
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-31 DOI: 10.1002/fld.5323
W. Düsterhöft‐Wriggers, S. Schubert, T. Rung
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.
以大变形为特征的粒状流问题广泛应用于各种领域,包括海岸工程和岩土工程。本文讨论了应用由德鲁克-普拉格屈服准则扩展的刚-完全塑性两相模型,用有限体积流动求解器(FV)模拟颗粒介质的问题。该模型采用宾汉流体和欧拉应变测量相结合的方法来评估颗粒坝滑体的破坏区域。流体体积(VoF)整体法用于区分空气相和颗粒相,二者均受不可压缩纳维-斯托克斯方程支配。该数值框架可对大规模应用中的大位移和任意形状进行建模。所显示的验证和检验主要针对具有不同休止角的非粘性和粘性材料的刚性-完全塑性材料模型。结果表明,预测的土壤表面和应变结果与实验和数值数据非常吻合。
{"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}
引用次数: 0
A spectral element discretization for quasi‐static magnetohydrodynamic flows 准静态磁流体动力学流的谱元离散法
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-19 DOI: 10.1002/fld.5321
Mattias Brynjell‐Rahkola
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}
引用次数: 0
Large Courant–Friedrichs–Lewy explicit scheme for one‐dimensional hyperbolic conservation laws 一维双曲守恒定律的大型库朗-弗里德里希斯-路维显式方案
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-05 DOI: 10.1002/fld.5322
Vincent Guinot, Antoine Rousseau
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.
本文介绍了一种大型库朗-弗里德里希斯-路维(CFL)算法,用于双曲守恒定律系统的显式有限体积求解,重点是浅水方程。通量计算中使用的黎曼问题是通过扩展到多个计算单元的平均核确定的。通常的 CFL 稳定性约束被一个涉及核支持大小的约束所取代。这使得该方法在计算单元大小方面无条件稳定,允许在不改变求解稳定性的情况下对计算网格进行任意程度的局部细化。针对带有地形源项的浅水方程,详细介绍了该方法的实际应用。计算实例报告了该方法在线性平流、伯格斯和浅水方程中的应用。在急剧的底部不连续的情况下,需要对几何源项进行改进的、均衡的离散处理。
{"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}
引用次数: 0
A hybrid marching cubes based IsoAlpha method for interface reconstruction 基于 IsoAlpha 的混合行进立方体界面重建法
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-04 DOI: 10.1002/fld.5320
G.S. Abhishek, Shyamprasad Karagadde
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.
摘要 在模拟两相流时,准确表示界面至关重要。有一类界面重建方法是基于等值面提取,其中涉及一种非迭代、基于插值的方法。事实证明,这些方法比传统的 PLIC 方案要快一个数量级。在这项工作中,我们提出了一种新的基于等值面提取的界面重建方案,该方案基于行进立方体算法(MC),该算法通常用于计算机图形学中的等值面可视化。MC 算法先将单个网格单元中所有可能的界面配置列出并归类到一个查找表(LUT)中,从而使这种方法既快速又稳健。我们还证明,对于某些界面配置,从单元体积分数获得等值的逆问题并不是可投射的,在处理这些情况时需要特殊处理。然后,我们通过二维和三维静态/动态界面重建的基准案例展示了该方法的能力。
{"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}
引用次数: 0
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 γ-R˜eθt$$ gamma -tilde{R}{e}_{theta t}$$ 过渡模型与可压缩流的 Spalart-Allmaras 模型耦合的连续邻接法
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-07-02 DOI: 10.1002/fld.5319
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.
在文献中首次开发并评估了可压缩流体过渡流的连续邻接法。在可压缩 Navier-Stokes 方程与 Spalart-Allmaras 湍流模型和过渡模型(非平滑和平滑三种变体)耦合的问题中,基于连续邻接法计算了空气动力目标函数(空气动力)相对于设计变量的梯度。事实证明,开发平稳过渡模型变体的邻接法是有益的。计算出的灵敏度导数的准确性与有限差分法进行了验证。编程在内部顶点为中心的有限体积代码中进行,可在 GPU 上高效运行。提出的连续邻接法被用于二维和三维气动外形优化问题,即 NLF(1)-0416 孤立翼面和 ONERA M6 机翼的约束优化。评估了 "冻结过渡"(假设不求解过渡模型的邻接方程)或 "冻结湍流"(另外忽略湍流模型的邻接方程)的影响;结果表明,这两种情况都会导致敏感度不准确。
{"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}
引用次数: 0
A modified lattice Boltzmann approach based on radial basis function approximation for the non‐uniform rectangular mesh 基于非均匀矩形网格径向基函数近似的改进型格子波尔兹曼方法
IF 1.8 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-06-25 DOI: 10.1002/fld.5318
X. Hu, J. M. Bergadà, D. Li, W. M. Sang, B. An
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.
我们提出了一种基于径向基函数近似(RBF-LBM)的非均匀矩形网格的新型格子波尔兹曼方法。非均匀矩形网格是局部网格细化的良好选择,特别是对于壁边界和流动量变化密集的流动区域。这样可以减少网格单元总数,从而降低计算成本,提高计算效率。但非均匀矩形网格的网格结构已不再适用于基于著名的 BGK 碰撞流演化的经典晶格玻尔兹曼方法(CLBM)。因此,本研究受到插值补充 LBM(ISLBM)方法的启发。本手稿对 ISLBM 算法进行了改进,通过径向基函数近似代替拉格朗日插值方案,将其发展为一种新型 LBM 方法。新方法在稳态和非稳态周期解中都得到了验证。讨论了径向基函数近似与拉格朗日插值之间的比较。结果发现,新方法在计算精度和效率方面都有很好的表现。证明了非均匀矩形网格允许网格细化,同时获得精确的流动预测。
{"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}
引用次数: 0
Decoupled and unconditionally stable iteration method for stationary Navier–Stokes equations 静态 Navier-Stokes 方程的解耦合和无条件稳定迭代法
IF 1.7 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-06-17 DOI: 10.1002/fld.5317
Jianhua Chen, Yingying Jiang, Guo-Dong Zhang

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 u$$ boldsymbol{u} $$ and pressure p$$ p $$ are coupled together at each iteration. By treating pressure p$$ p $$ 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 (t+$$ tto +infty $$). 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.

众所周知,静态纳维-斯托克斯方程的奥森迭代是无条件稳定的。然而,它是一种耦合型方案,每次迭代时速度 u $$ boldsymbol{u} $$ 和压力 p $$ p $$ 都耦合在一起。明确处理压力 p $$ p $ 会导致解耦迭代,但这种处理方法并不稳定。在本文中,我们通过采用压力投影法来求解时间扰动 Navier-Stokes 系统,构建了一种解耦且无条件稳定的迭代方法来求解静态 Navier-Stokes 方程,该方法的解近似于随时间变化的稳态解 ( t → + ∞ $$ tto +infty $$ )。我们还严格证明了它的无条件稳定性。数值模拟证明,与广泛使用的 T-S 和 Oseen 迭代法相比,我们的迭代法更有效、更稳定,而且可以解决高雷诺数的流体流动问题。
{"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":null,"pages":null},"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}
引用次数: 0
A general pressure equation based method for incompressible two-phase flows 基于一般压力方程的不可压缩两相流方法
IF 1.7 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-06-10 DOI: 10.1002/fld.5316
Hormuzd Bodhanwalla, Dheeraj Raghunathan, Y. Sudhakar

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.

我们提出了一种全显式、无迭代、弱可压缩的方法来模拟不相溶的不可压缩两相流。为了更新压力,我们避开了计算成本高昂的泊松方程,采用了显式求解的一般压力方程。此外,我们还使用了扩散性较低的流体体积代数方法作为界面捕捉技术,并且为了提高并行计算的可扩展性,我们使用算子分割方法对该技术进行了时间离散化。我们的方法是完全显式的,并且通过简单的局部空间离散化就能保持稳定,因此很容易实现。我们模拟了几种二维和三维典型两相流。定性和定量结果证明,我们的方法能够准确处理涉及一系列密度和粘度比以及表面张力效应的问题。
{"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":null,"pages":null},"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}
引用次数: 0
A modified fifth-order WENO-Z scheme based on the weights of the reformulated adaptive order WENO scheme 基于重构自适应阶 WENO 方案权重的改进型五阶 WENO-Z 方案
IF 1.7 4区 工程技术 Q3 COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS Pub Date : 2024-06-04 DOI: 10.1002/fld.5314
Yize Wang, Kunlei Zhao, Li Yuan

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.

通过类比重新制定的五阶自适应阶次(AO)WENO 方案的非归一化权值,我们提出了一种改进的五阶 WENO-Z 方案。我们的研究表明,如果将原始的五阶 WENO-AO 方案改写成传统 WENO 组合的形式,则所得到的非归一化权重可分为三个部分:常数项、局部模版平滑度测量项和全局模版平滑度测量项。为了利用后两个项构建性能更高的改进型 WENO-Z 方案,我们改变了第三项的形式,并引入了一个自适应缩放因子来调整第二项和第三项的贡献。数值示例表明,修正的五阶 WENO-Z 方案具有对平滑区域分辨率高、对不连续区域捕捉敏锐的优点,与最近开发的 WENO-Z+ 和 WENO-Z+M 方案相比,它在具有小尺度结构的冲击流中可以获得明显更好的结果。
{"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":null,"pages":null},"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}
引用次数: 0
期刊
International Journal for Numerical Methods in Fluids
全部 Acc. Chem. Res. ACS Applied Bio Materials ACS Appl. Electron. Mater. ACS Appl. Energy Mater. ACS Appl. Mater. Interfaces ACS Appl. Nano Mater. ACS Appl. Polym. Mater. ACS BIOMATER-SCI ENG ACS Catal. ACS Cent. Sci. ACS Chem. Biol. ACS Chemical Health & Safety ACS Chem. Neurosci. ACS Comb. Sci. ACS Earth Space Chem. ACS Energy Lett. ACS Infect. Dis. ACS Macro Lett. ACS Mater. Lett. ACS Med. Chem. Lett. ACS Nano ACS Omega ACS Photonics ACS Sens. ACS Sustainable Chem. Eng. ACS Synth. Biol. Anal. Chem. BIOCHEMISTRY-US Bioconjugate Chem. BIOMACROMOLECULES Chem. Res. Toxicol. Chem. Rev. Chem. Mater. CRYST GROWTH DES ENERG FUEL Environ. Sci. Technol. Environ. Sci. Technol. Lett. Eur. J. Inorg. Chem. IND ENG CHEM RES Inorg. Chem. J. Agric. Food. Chem. J. Chem. Eng. Data J. Chem. Educ. J. Chem. Inf. Model. J. Chem. Theory Comput. J. Med. Chem. J. Nat. Prod. J PROTEOME RES J. Am. Chem. Soc. LANGMUIR MACROMOLECULES Mol. Pharmaceutics Nano Lett. Org. Lett. ORG PROCESS RES DEV ORGANOMETALLICS J. Org. Chem. J. Phys. Chem. J. Phys. Chem. A J. Phys. Chem. B J. Phys. Chem. C J. Phys. Chem. Lett. Analyst Anal. Methods Biomater. Sci. Catal. Sci. Technol. Chem. Commun. Chem. Soc. Rev. CHEM EDUC RES PRACT CRYSTENGCOMM Dalton Trans. Energy Environ. Sci. ENVIRON SCI-NANO ENVIRON SCI-PROC IMP ENVIRON SCI-WAT RES Faraday Discuss. Food Funct. Green Chem. Inorg. Chem. Front. Integr. Biol. J. Anal. At. Spectrom. J. Mater. Chem. A J. Mater. Chem. B J. Mater. Chem. C Lab Chip Mater. Chem. Front. Mater. Horiz. MEDCHEMCOMM Metallomics Mol. Biosyst. Mol. Syst. Des. Eng. Nanoscale Nanoscale Horiz. Nat. Prod. Rep. New J. Chem. Org. Biomol. Chem. Org. Chem. Front. PHOTOCH PHOTOBIO SCI PCCP Polym. Chem.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
0
微信
客服QQ
Book学术公众号 扫码关注我们
反馈
×
意见反馈
请填写您的意见或建议
请填写您的手机或邮箱
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
现在去查看 取消
×
提示
确定
Book学术官方微信
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术
文献互助 智能选刊 最新文献 互助须知 联系我们:info@booksci.cn
Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。
Copyright © 2023 Book学术 All rights reserved.
ghs 京公网安备 11010802042870号 京ICP备2023020795号-1