Pub Date : 2023-01-02DOI: 10.1080/10618562.2023.2242271
J. Muralha, L. Eça, C. Klaij
This paper presents Solution Verification exercises with the pressure-based compressible flow solver ReFRESCO for five test cases available in the NASA Turbulence Modeling Resource: the two-dimensional flows over a flat plate, a bump-in-channel, a DSMA661 airfoil and a multi-element airfoil and the three-dimensional flow of a bump-in-channel. Simulations are performed with the Favre-averaged continuity and Navier-Stokes equations using the Spalart & Allmaras turbulence model. ReFRESCO results are compared with reference data from density-based compressible flow solvers (CFL3D and FUN3D). two aspects of the implementation of the turbulence model are addressed: the calculation of the distance to the wall and the discretization scheme used in the convective terms of the turbulence model transport equation. Results of this study show perfect consistency with the reference data for the test cases that are not affected by the determination of the distance to the wall.
{"title":"Verification of a Pressure-Based Compressible Flow Solver","authors":"J. Muralha, L. Eça, C. Klaij","doi":"10.1080/10618562.2023.2242271","DOIUrl":"https://doi.org/10.1080/10618562.2023.2242271","url":null,"abstract":"This paper presents Solution Verification exercises with the pressure-based compressible flow solver ReFRESCO for five test cases available in the NASA Turbulence Modeling Resource: the two-dimensional flows over a flat plate, a bump-in-channel, a DSMA661 airfoil and a multi-element airfoil and the three-dimensional flow of a bump-in-channel. Simulations are performed with the Favre-averaged continuity and Navier-Stokes equations using the Spalart & Allmaras turbulence model. ReFRESCO results are compared with reference data from density-based compressible flow solvers (CFL3D and FUN3D). two aspects of the implementation of the turbulence model are addressed: the calculation of the distance to the wall and the discretization scheme used in the convective terms of the turbulence model transport equation. Results of this study show perfect consistency with the reference data for the test cases that are not affected by the determination of the distance to the wall.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2023-01-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81391573","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":"Influence of the Downstream Vehicle Length on Train Aerodynamics Subjected to Crosswind","authors":"Zhuang Tianci, Li Wenhui, Liu Tanghong","doi":"10.36959/717/661","DOIUrl":"https://doi.org/10.36959/717/661","url":null,"abstract":"","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80753220","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":"Thermal Radiation, Chemical Reaction and Viscous Dissipation Effects on MHD Mixed Convection Flow of Micro Polar Fluid with Stretching Surface in the Presence of Heat Generation/Absorption","authors":"Zigta Binyam","doi":"10.36959/717/662","DOIUrl":"https://doi.org/10.36959/717/662","url":null,"abstract":"","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-12-31","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"84918205","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}
Pub Date : 2022-11-26DOI: 10.1080/10618562.2023.2225416
Aref Zanjani, A. Tahsini, Kimia Sadafi, Fatemeh Ghavidel Mangodeh
Shape optimisation of supersonic nozzles is of crucial importance in designing propulsion systems and space thrusters. In order to optimise the profile of a supersonic nozzle, the properties of the flow inside the nozzle should be obtained. This paper proposes and verifies a new methodology for analysing flows and designing supersonic nozzles. Flow analysis has been conducted using the method of characteristics, Ansys Fluent and convolutional neural networks. It is shown that deep convolutional neural networks can reach high levels of accuracy in predicting supersonic flow behaviour inside the nozzle. Also, shape optimisation of the supersonic nozzle has been conducted using the genetic algorithm in Ansys Fluent and artificial neural networks. The proposed ANN can optimise the shape of a supersonic nozzle for the given throat diameter, outlet diameter and nozzle length with high accuracy.
{"title":"Shape Optimization and Flow Analysis of Supersonic Nozzles Using Deep Learning","authors":"Aref Zanjani, A. Tahsini, Kimia Sadafi, Fatemeh Ghavidel Mangodeh","doi":"10.1080/10618562.2023.2225416","DOIUrl":"https://doi.org/10.1080/10618562.2023.2225416","url":null,"abstract":"Shape optimisation of supersonic nozzles is of crucial importance in designing propulsion systems and space thrusters. In order to optimise the profile of a supersonic nozzle, the properties of the flow inside the nozzle should be obtained. This paper proposes and verifies a new methodology for analysing flows and designing supersonic nozzles. Flow analysis has been conducted using the method of characteristics, Ansys Fluent and convolutional neural networks. It is shown that deep convolutional neural networks can reach high levels of accuracy in predicting supersonic flow behaviour inside the nozzle. Also, shape optimisation of the supersonic nozzle has been conducted using the genetic algorithm in Ansys Fluent and artificial neural networks. The proposed ANN can optimise the shape of a supersonic nozzle for the given throat diameter, outlet diameter and nozzle length with high accuracy.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87690820","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}
Pub Date : 2022-11-26DOI: 10.1080/10618562.2023.2221645
Dawei Peng, Lanhao Zhao, Chuanyuan Zhou, Jia Mao
Pile safety has received increasing attention in marine engineering, especially in the field of local scour. In this paper, a finite element numerical model is established for local scour around a cylinder in steady currents. The flow is described by unsteady Reynolds–averaged Navier–Stokes equations with a traditional turbulent closure model. The proposed scour model takes bed load into account. The Exner equation is solved to determine the bed variation and the moving mesh approach is used to capture the evolution of the bed. When the resulting slope exceeds the angle of repose, a novel sand-slide model based on Rodrigues' rotation formula is used to prevent simulation distortion. All the equations are discretized by the two-step Taylor–Galerkin algorithm, and the resulting approach is fast to implement with second-order accuracy in space. The numerical results are found to be in good agreement with the experimental data.
桩的安全问题在海洋工程领域,尤其是局部冲刷领域受到越来越多的关注。本文建立了稳定流条件下圆柱局部冲刷的有限元数值模型。用非定常reynolds - average Navier-Stokes方程和传统的湍流闭包模型来描述流动。提出的冲刷模型考虑了河床荷载。通过求解Exner方程确定床层的变化,采用移动网格法捕捉床层的演变。当产生的坡度超过休止角时,采用基于Rodrigues旋转公式的新型滑坡模型来防止模拟失真。采用两步Taylor-Galerkin算法对所有方程进行离散化,得到的方法在空间上具有二阶精度,实现速度快。数值计算结果与实验数据吻合较好。
{"title":"Finite Element Numerical Simulation of Local Scour of a Three-Dimensional Cylinder under Steady Flow","authors":"Dawei Peng, Lanhao Zhao, Chuanyuan Zhou, Jia Mao","doi":"10.1080/10618562.2023.2221645","DOIUrl":"https://doi.org/10.1080/10618562.2023.2221645","url":null,"abstract":"Pile safety has received increasing attention in marine engineering, especially in the field of local scour. In this paper, a finite element numerical model is established for local scour around a cylinder in steady currents. The flow is described by unsteady Reynolds–averaged Navier–Stokes equations with a traditional turbulent closure model. The proposed scour model takes bed load into account. The Exner equation is solved to determine the bed variation and the moving mesh approach is used to capture the evolution of the bed. When the resulting slope exceeds the angle of repose, a novel sand-slide model based on Rodrigues' rotation formula is used to prevent simulation distortion. All the equations are discretized by the two-step Taylor–Galerkin algorithm, and the resulting approach is fast to implement with second-order accuracy in space. The numerical results are found to be in good agreement with the experimental data.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"75219110","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}
Pub Date : 2022-11-26DOI: 10.1080/10618562.2023.2237898
Xiaogang Li, Tian Xia, Yuxi Deng, Siqi Yang, Yonbin Ge
In this work, a new, improved third-order finite difference weighted essentially non-oscillatory scheme is presented for one- and two-dimensional hyperbolic conservation laws and associated problems. The parameter p which is regulate dissipation is introduced in the nonlinear weights in the framework of the conventional WENO-Z scheme, and the higher-order global smoothness indicator is obtained by the idea of Wang [Wang, Y. H., Y. L. Du, K. L. Zhao and L. Yuan. 2020. ‘A Low-dissipation Third-order Weighted Essentially Nonoscillatory Scheme with a New Reference Smoothness Indicator’. International Journal for Numerical Methods in Fluids. 92 (9): 1212–1234.], the sufficient condition of nonlinear weights is proved by using Taylor expansions. Finally, the value range of parameter p is obtained. The proposed scheme is verified to achieve the optimal order near critical points by linear convection equations with different initial values, and the high-resolution characteristic of the present scheme is proved on a variety of one- and two- dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with the other third-order WENO schemes.
在这项工作中,提出了一种新的改进的三阶有限差分加权本质非振荡格式,用于一维和二维双曲守恒律和相关问题。在传统WENO-Z格式框架的非线性权值中引入了调节耗散参数p,并采用Wang的思想得到了高阶全局平滑指标[j] . Wang, yyh ., Du yl ., Zhao kl ., Yuan L. 2020。一种具有新的参考平滑指标的低耗散三阶加权基本非振荡格式。流体力学与工程学报,2009(9):1212-1234。],用泰勒展开式证明了非线性权值存在的充分条件。最后得到参数p的取值范围。通过不同初值的线性对流方程,验证了该格式在临界点附近的最优阶,并在各种一维和二维标准数值算例上证明了该格式的高分辨率特性。数值结果表明,与其他三阶WENO算法相比,该算法具有更好的性能。
{"title":"A New Third-Order Finite Difference WENO Scheme to Improve Convergence Rate at Critical Points","authors":"Xiaogang Li, Tian Xia, Yuxi Deng, Siqi Yang, Yonbin Ge","doi":"10.1080/10618562.2023.2237898","DOIUrl":"https://doi.org/10.1080/10618562.2023.2237898","url":null,"abstract":"In this work, a new, improved third-order finite difference weighted essentially non-oscillatory scheme is presented for one- and two-dimensional hyperbolic conservation laws and associated problems. The parameter p which is regulate dissipation is introduced in the nonlinear weights in the framework of the conventional WENO-Z scheme, and the higher-order global smoothness indicator is obtained by the idea of Wang [Wang, Y. H., Y. L. Du, K. L. Zhao and L. Yuan. 2020. ‘A Low-dissipation Third-order Weighted Essentially Nonoscillatory Scheme with a New Reference Smoothness Indicator’. International Journal for Numerical Methods in Fluids. 92 (9): 1212–1234.], the sufficient condition of nonlinear weights is proved by using Taylor expansions. Finally, the value range of parameter p is obtained. The proposed scheme is verified to achieve the optimal order near critical points by linear convection equations with different initial values, and the high-resolution characteristic of the present scheme is proved on a variety of one- and two- dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with the other third-order WENO schemes.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80247655","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}
Pub Date : 2022-11-26DOI: 10.1080/10618562.2023.2229250
N. Pham-Sy, C. Tran
A non-overlapping domain decomposition-based parallel algorithm coupled with a compact local integrated radial basis function (CLIRBF) method is developed for solving Navier-Stokes equations. For this approach, a problem is divided into subdomains. In each sub-domain, a CLIRBF scheme is applied to solve the Navier-Stokes equations of flows. A relaxation factor is used at the interface between sub-domains to ensure the quick convergence of the present method. The Bitmap termination detection technique is introduced to complete the global termination. The present approach is verified using two fluid flow problems: the lid-driven cavity and the natural convection in concentric annuli flow. The numerical results have demonstrated the efficiency of the present parallel method compared with the corresponding sequential one and other published methods. Especially, super-linear speed-up was achieved for several CPUs. In terms of accuracy, the obtained results are in very good agreement with benchmark results.
{"title":"Parallel Computation Using Non-Overlapping Domain Decomposition Coupled with Compact Local Integrated RBF for Navier–Stokes Equations","authors":"N. Pham-Sy, C. Tran","doi":"10.1080/10618562.2023.2229250","DOIUrl":"https://doi.org/10.1080/10618562.2023.2229250","url":null,"abstract":"A non-overlapping domain decomposition-based parallel algorithm coupled with a compact local integrated radial basis function (CLIRBF) method is developed for solving Navier-Stokes equations. For this approach, a problem is divided into subdomains. In each sub-domain, a CLIRBF scheme is applied to solve the Navier-Stokes equations of flows. A relaxation factor is used at the interface between sub-domains to ensure the quick convergence of the present method. The Bitmap termination detection technique is introduced to complete the global termination. The present approach is verified using two fluid flow problems: the lid-driven cavity and the natural convection in concentric annuli flow. The numerical results have demonstrated the efficiency of the present parallel method compared with the corresponding sequential one and other published methods. Especially, super-linear speed-up was achieved for several CPUs. In terms of accuracy, the obtained results are in very good agreement with benchmark results.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-11-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91306707","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}
Pub Date : 2022-10-21DOI: 10.1080/10618562.2023.2189704
G. Zhao, Chengwen Zhong, Sha Liu, Yong Wang, Congshan Zhuo
A gas-kinetic scheme (GKS) with kinetic boundary condition based on unstructured mesh is present here. In the GKS method, the solid wall boundary conditions can be constructed by virtue of the gas distribution function, which is similar to the diffuse-scattering rule used in the other kinetic schemes. The kinetic boundary condition has a concise form and easy to implement. The use of unstructured mesh expands the adaptability of GKS to simulate the flows with complex geometry. The kinetic boundary condition can recover to the non-slip boundary condition in the continuum regime. In the slip regime, the slip velocity can be accurately predicted by kinetic boundary condition, which turns into the slip boundary condition. The use of kinetic boundary condition improves the calculation results of GKS in near-continuum flow. A series of test cases, from incompressible to compressible flow with a wide range of Knudsen number, are investigated to demonstrate the performance of kinetic boundary condition in near-continuum flow, which can provide a reference for the construction and optimisation for GKS-based multi-scale hybrid algorithms.
{"title":"Application of Gas-Kinetic Scheme for Continuum and Near-Continuum Flow on Unstructured Mesh","authors":"G. Zhao, Chengwen Zhong, Sha Liu, Yong Wang, Congshan Zhuo","doi":"10.1080/10618562.2023.2189704","DOIUrl":"https://doi.org/10.1080/10618562.2023.2189704","url":null,"abstract":"A gas-kinetic scheme (GKS) with kinetic boundary condition based on unstructured mesh is present here. In the GKS method, the solid wall boundary conditions can be constructed by virtue of the gas distribution function, which is similar to the diffuse-scattering rule used in the other kinetic schemes. The kinetic boundary condition has a concise form and easy to implement. The use of unstructured mesh expands the adaptability of GKS to simulate the flows with complex geometry. The kinetic boundary condition can recover to the non-slip boundary condition in the continuum regime. In the slip regime, the slip velocity can be accurately predicted by kinetic boundary condition, which turns into the slip boundary condition. The use of kinetic boundary condition improves the calculation results of GKS in near-continuum flow. A series of test cases, from incompressible to compressible flow with a wide range of Knudsen number, are investigated to demonstrate the performance of kinetic boundary condition in near-continuum flow, which can provide a reference for the construction and optimisation for GKS-based multi-scale hybrid algorithms.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"76193943","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}
Pub Date : 2022-10-21DOI: 10.1080/10618562.2023.2189247
D. J. Lusher, G. Coleman
ABSTRACT Direct numerical simulation is used to determine the turbulent Prandtl number above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal wall scaling proposed in 1995 by [Huang, P. G., G. N. Coleman, and P. Bradshaw. 1995. “Compressible Turbulent Channel Flows: DNS Results and Modelling.” Journal of Fluid Mechanics 305: 185–218. doi:10.1017/S0022112095004599.], with only a weak dependence on . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.
摘要采用直接数值模拟方法确定了一类低声速通道中冷(等温)壁面和热(绝热)壁面上的湍流普朗特数。考虑了一个平均温度/密度变化范围,对应于有效/边缘马赫数在1.1到2.2之间,基于壁面变量的雷诺数在73到3800之间。绝热条件是一个特别有趣的新特征。在等温壁面和绝热壁面以上,离壁面的差值均接近0.85。[Huang, P. G. Coleman, and P. Bradshaw. 1995]在1995年提出的被动标量模拟中,近壁面剖面的变化作为半局部壁面结垢的函数而崩溃。可压缩湍流通道流动:DNS结果和建模。流体力学学报(自然科学版);doi: 10.1017 / S0022112095004599。对…只有微弱的依赖。这导致了一个相当简单的传热模型的提议,附属于涡流粘度模型。
{"title":"Numerical Study of Compressible Wall-Bounded Turbulence – the Effect of Thermal Wall Conditions on the Turbulent Prandtl Number in the Low-Supersonic Regime","authors":"D. J. Lusher, G. Coleman","doi":"10.1080/10618562.2023.2189247","DOIUrl":"https://doi.org/10.1080/10618562.2023.2189247","url":null,"abstract":"ABSTRACT Direct numerical simulation is used to determine the turbulent Prandtl number above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal wall scaling proposed in 1995 by [Huang, P. G., G. N. Coleman, and P. Bradshaw. 1995. “Compressible Turbulent Channel Flows: DNS Results and Modelling.” Journal of Fluid Mechanics 305: 185–218. doi:10.1017/S0022112095004599.], with only a weak dependence on . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88663788","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}
Pub Date : 2022-10-21DOI: 10.1080/10618562.2023.2202391
Kai Yang, T. Aoki
A consistent and conservative formulation for mass and momentum transport is proposed in the context of simulating incompressible two-phase flows by using weakly compressible method. Combined with the evolving pressure projection method to prevent oscillation of the solution induced by the acoustic wave, this solver aims at a robust and accurate computation of violent two-phase flows with a high density ratio, while taking advantage of fully explicit time integration of the weakly compressible Navier–Stokes equations. Coupled with the volume of fluid method for capturing interfaces, the mass and momentum fluxes are evaluated in a consistent manner using the finite volume method. In addition, a special implementation of the pressure projection is devised to avoid velocity-pressure decoupling on a collocated grid. The solver's accuracy and stability are demonstrated through various two-phase flow simulations, including dam break and liquid jet atomization scenarios, emphasizing its momentum-conserving properties.
{"title":"A Momentum-Conserving Weakly Compressible Navier–Stokes Solver for Simulation of Violent Two-Phase Flows with High Density Ratio","authors":"Kai Yang, T. Aoki","doi":"10.1080/10618562.2023.2202391","DOIUrl":"https://doi.org/10.1080/10618562.2023.2202391","url":null,"abstract":"A consistent and conservative formulation for mass and momentum transport is proposed in the context of simulating incompressible two-phase flows by using weakly compressible method. Combined with the evolving pressure projection method to prevent oscillation of the solution induced by the acoustic wave, this solver aims at a robust and accurate computation of violent two-phase flows with a high density ratio, while taking advantage of fully explicit time integration of the weakly compressible Navier–Stokes equations. Coupled with the volume of fluid method for capturing interfaces, the mass and momentum fluxes are evaluated in a consistent manner using the finite volume method. In addition, a special implementation of the pressure projection is devised to avoid velocity-pressure decoupling on a collocated grid. The solver's accuracy and stability are demonstrated through various two-phase flow simulations, including dam break and liquid jet atomization scenarios, emphasizing its momentum-conserving properties.","PeriodicalId":56288,"journal":{"name":"International Journal of Computational Fluid Dynamics","volume":null,"pages":null},"PeriodicalIF":1.3,"publicationDate":"2022-10-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"88643613","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}