In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether t...In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether those r-factors would be appropriate and effective for the vertex-centered（VC） FVM. In the paper, we collected five kinds of r-factor formulae and found out that only three of those, respectively by Bruner（1996）, Darwish and Moukalled（2003） and Cassuli and Zanolli（2005） can be formally extended to a context of the VC FVM. Numerical tests indicate that the TVD schemes and r-factors, after being extended and introduced to a context of the VC FVM, maintained their similar characteristics as in a context of the CC FVM. However, when the gradient-based r-factors and the SUPERBEE scheme were applied simultaneously, non-physical oscillations near the sharp step would appear. In the transient case, the oscillations were weaker in a context of the VC FVM than those in a context of the CC FVM, while the effect was reversed in the steady case. To eliminate disadvantages in the gradient-based r-factor formula, a new modification method by limiting values on the virtual node, namely Фu in the paper, was validated by the tests to effectively dissipate spurious oscillations.展开更多
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.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.41306078 and 41301414)the National Engineering Research Center for Inland Waterway Regulation and Key Laboratory of Hydraulic and Waterway Engineering of the Ministry of Education Program(Grant No.SLK2016B03)the Key Laboratory of the Inland Waterway Regulation of the Ministry of Transportation Program(Grant No.NHHD-201514)
文摘In most TVD schemes, the r-factors were proposed according to the cell-centered（CC） finite volume method（FVM） framework for the numerical approximation to the convective term. However, it is questionable whether those r-factors would be appropriate and effective for the vertex-centered（VC） FVM. In the paper, we collected five kinds of r-factor formulae and found out that only three of those, respectively by Bruner（1996）, Darwish and Moukalled（2003） and Cassuli and Zanolli（2005） can be formally extended to a context of the VC FVM. Numerical tests indicate that the TVD schemes and r-factors, after being extended and introduced to a context of the VC FVM, maintained their similar characteristics as in a context of the CC FVM. However, when the gradient-based r-factors and the SUPERBEE scheme were applied simultaneously, non-physical oscillations near the sharp step would appear. In the transient case, the oscillations were weaker in a context of the VC FVM than those in a context of the CC FVM, while the effect was reversed in the steady case. To eliminate disadvantages in the gradient-based r-factor formula, a new modification method by limiting values on the virtual node, namely Фu in the paper, was validated by the tests to effectively dissipate spurious oscillations.
基金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.
