A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing e...A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.展开更多
Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisi...Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisit the question of non-oscillatory behavior of high-resolution central schemes in terms of the slope limiter proposed by van Albada et al.(Astron Astrophys 108:76–84,1982).The van Albada(vA)limiter is smoother near extrema,and consequently,in many cases,it outperforms the results obtained using the standard minmod limiter.In particular,we prove that the vA limiter ensures the one-dimensional Total-Variation Diminishing(TVD)stability and demonstrate that it yields noticeable improvement in computation of one-and two-dimensional systems.展开更多
Flow transition from laminar to turbulent is prerequisite to decide whereabouts to apply surface flow control techniques. This appears missing in a number of works in which the control effects were merely investigated...Flow transition from laminar to turbulent is prerequisite to decide whereabouts to apply surface flow control techniques. This appears missing in a number of works in which the control effects were merely investigated without getting insight into alteration of transition position. The aim of this study is to capture the correct position of transition over NACA0012 aerofoil at different angles of attack. Firstly, an implicit, time marching, high resolution total variation diminishing （TVD） scheme was developed to solve the governing Navier-Stokes equations for compressible fluid flows around aerofoil sections to obtain velocity profiles around the aerofoil surfaces. Secondly, the linear instability solver based on the Orr-Sommerfeld equations and the eg methods were developed to calculate the onset of transition over the aerofoil surfaces. For the low subsonic Mach number of 0.16, the accuracy of the compressible solutions was assessed by some available experimental results of low speed incompressible flows. In all cases, transition positions were accurately predicted which shows applicability and superiority of the present work to be extended for higher Mach number compressible flows. Here, transition prediction methodology is described and the results of this analysis without active flow control or separation are presented.展开更多
A finite-volume Total Variation Diminishing (TVD) scheme is presented formodeling dam-break flows in open channels. This method is used for solving the 2D shallow waterequations on arbitrary quadrilateral meshes, base...A finite-volume Total Variation Diminishing (TVD) scheme is presented formodeling dam-break flows in open channels. This method is used for solving the 2D shallow waterequations on arbitrary quadrilateral meshes, based upon a second-order hybrid TVD scheme with anoptimum-selected limiter in the space discretization and a two-step Runge-Kutta approach in the timediscretization. Verification for a circular dam-break problem is carried out by comparing thepresent results with others and very good agreement is shown. The present algorithm is then used topredict dam-break flow characteristics in open channels such as in furcated channels. Morecomplicated unsteady flow characteristics in these furcated channels than in the regular channelsstudied previously can observed in this work.展开更多
Flows around a circular cylinder displaying an unsteady vortex shedding process at the Reynolds numbers of 1000,3900 and 1×104 are studied using a finite-volume Total Variation Diminishing（TVD） scheme for solvi...Flows around a circular cylinder displaying an unsteady vortex shedding process at the Reynolds numbers of 1000,3900 and 1×104 are studied using a finite-volume Total Variation Diminishing（TVD） scheme for solving the Unsteady Reynolds-Averaged Navier-Stokes（URANS） equations.An Elemental Velocity Vector Transformation（EVVT） approach is proposed for the local normal and tangential velocity transformation at the interfaces of main and satellite elements.The presented method is validated by comparing with the available experimental data and numerical results.It is shown that the two-dimensional TVD finite volume method with the Renormalization Group（RNG） turbulence model can be used to determine hydrodynamic forces and captures vortex shedding characteristics very well.展开更多
A finite-difference Total Variation Diminishing （TVD） numerical simulation model for coupling the Reynolds Averaged Navier-Stokes （RANS） equations, pressure-relative continuity equation and various k-εturbulence ...A finite-difference Total Variation Diminishing （TVD） numerical simulation model for coupling the Reynolds Averaged Navier-Stokes （RANS） equations, pressure-relative continuity equation and various k-εturbulence models was developed to solve the incompressible flow based on the pseudo-compressibility method. The hyperbolicity of all these equations was studied and the discretization of the fully coupling equations with all the primal variables and source terms were made in this article. Numerical simulation for modeling the flow around a ground-mounted square rib was implemented and validated by comparing with the published wind tunnel experimental data. It is shown that such a numerical simulation method with a proper turbulence model has a very good accuracy to simulate the flow around a surface-mounted rib. It is concluded that the Renormalization Group （RNG） and Chen-Kim k-εturbulence models have much better ability to predict the characteristics of the vortex structure and flow separation than the standard k-εmodel.展开更多
为研究气体离心机取料支臂附近的流动状况,数值模拟了取料支臂附近的流场。采用矢量分裂的有限体积法求解不同出流条件下的N av ier-S tokes方程。采用二阶总变差减小原理,捕捉激波,同时将支臂外部和内部的流动结合起来进行计算。讨论...为研究气体离心机取料支臂附近的流动状况,数值模拟了取料支臂附近的流场。采用矢量分裂的有限体积法求解不同出流条件下的N av ier-S tokes方程。采用二阶总变差减小原理,捕捉激波,同时将支臂外部和内部的流动结合起来进行计算。讨论了激波和气体粘性对流场的影响。得到了不同取料情况下的流场分布图像。结果表明:所采用的计算格式很好地捕捉到了激波和边界层,验证了计算方法的有效性,为数值模拟实际离心机支臂附近的流动状况和支臂损耗的计算打下了基础。展开更多
基金supported by the National Natural Science Foundation of China(Grant No.51579034)the Open Fund of the Key Laboratory of Ocean Circulation and Waves,Chinese Academy of Sciences(Grant No.KLOCW1502)
文摘A total variation diminishing-weighted average flux （TVD-WAF）-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer （HLL） Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws （MUSCL）. The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.
基金Research was supported in part by the ONR Grant N00014-2112773.
文摘Slope limiters play an essential role in maintaining the non-oscillatory behavior of high-resolution methods for nonlinear conservation laws.The family of minmod limiters serves as the prototype example.Here,we revisit the question of non-oscillatory behavior of high-resolution central schemes in terms of the slope limiter proposed by van Albada et al.(Astron Astrophys 108:76–84,1982).The van Albada(vA)limiter is smoother near extrema,and consequently,in many cases,it outperforms the results obtained using the standard minmod limiter.In particular,we prove that the vA limiter ensures the one-dimensional Total-Variation Diminishing(TVD)stability and demonstrate that it yields noticeable improvement in computation of one-and two-dimensional systems.
文摘Flow transition from laminar to turbulent is prerequisite to decide whereabouts to apply surface flow control techniques. This appears missing in a number of works in which the control effects were merely investigated without getting insight into alteration of transition position. The aim of this study is to capture the correct position of transition over NACA0012 aerofoil at different angles of attack. Firstly, an implicit, time marching, high resolution total variation diminishing （TVD） scheme was developed to solve the governing Navier-Stokes equations for compressible fluid flows around aerofoil sections to obtain velocity profiles around the aerofoil surfaces. Secondly, the linear instability solver based on the Orr-Sommerfeld equations and the eg methods were developed to calculate the onset of transition over the aerofoil surfaces. For the low subsonic Mach number of 0.16, the accuracy of the compressible solutions was assessed by some available experimental results of low speed incompressible flows. In all cases, transition positions were accurately predicted which shows applicability and superiority of the present work to be extended for higher Mach number compressible flows. Here, transition prediction methodology is described and the results of this analysis without active flow control or separation are presented.
文摘A finite-volume Total Variation Diminishing (TVD) scheme is presented formodeling dam-break flows in open channels. This method is used for solving the 2D shallow waterequations on arbitrary quadrilateral meshes, based upon a second-order hybrid TVD scheme with anoptimum-selected limiter in the space discretization and a two-step Runge-Kutta approach in the timediscretization. Verification for a circular dam-break problem is carried out by comparing thepresent results with others and very good agreement is shown. The present algorithm is then used topredict dam-break flow characteristics in open channels such as in furcated channels. Morecomplicated unsteady flow characteristics in these furcated channels than in the regular channelsstudied previously can observed in this work.
基金supported by the National High Technology Research and Development Program of China (863 Program,Grant No. 2008AA09Z310)the Important National Scienceand Technology Specific Sub-Project (Grant No.2008ZX05026-001)
文摘Flows around a circular cylinder displaying an unsteady vortex shedding process at the Reynolds numbers of 1000,3900 and 1×104 are studied using a finite-volume Total Variation Diminishing（TVD） scheme for solving the Unsteady Reynolds-Averaged Navier-Stokes（URANS） equations.An Elemental Velocity Vector Transformation（EVVT） approach is proposed for the local normal and tangential velocity transformation at the interfaces of main and satellite elements.The presented method is validated by comparing with the available experimental data and numerical results.It is shown that the two-dimensional TVD finite volume method with the Renormalization Group（RNG） turbulence model can be used to determine hydrodynamic forces and captures vortex shedding characteristics very well.
文摘A finite-difference Total Variation Diminishing （TVD） numerical simulation model for coupling the Reynolds Averaged Navier-Stokes （RANS） equations, pressure-relative continuity equation and various k-εturbulence models was developed to solve the incompressible flow based on the pseudo-compressibility method. The hyperbolicity of all these equations was studied and the discretization of the fully coupling equations with all the primal variables and source terms were made in this article. Numerical simulation for modeling the flow around a ground-mounted square rib was implemented and validated by comparing with the published wind tunnel experimental data. It is shown that such a numerical simulation method with a proper turbulence model has a very good accuracy to simulate the flow around a surface-mounted rib. It is concluded that the Renormalization Group （RNG） and Chen-Kim k-εturbulence models have much better ability to predict the characteristics of the vortex structure and flow separation than the standard k-εmodel.
文摘为研究气体离心机取料支臂附近的流动状况,数值模拟了取料支臂附近的流场。采用矢量分裂的有限体积法求解不同出流条件下的N av ier-S tokes方程。采用二阶总变差减小原理,捕捉激波,同时将支臂外部和内部的流动结合起来进行计算。讨论了激波和气体粘性对流场的影响。得到了不同取料情况下的流场分布图像。结果表明:所采用的计算格式很好地捕捉到了激波和边界层,验证了计算方法的有效性,为数值模拟实际离心机支臂附近的流动状况和支臂损耗的计算打下了基础。
