scholarly journals Extending the global-direction stencil with “face-area-weighted centroid” to unstructured finite volume discretization from integral form

2020 ◽  
Vol 2 (1) ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

AbstractAccuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work [Kong, et al. Chin Phys B 29(10):100203, 2020], we explored a novel global-direction stencil and combined it with the face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Greatly inspired by the differential form, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization from integral form on both second and third-order finite volume solver. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has superiorities on both computational accuracy and convergence rate.

2020 ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

Abstract Accuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work [Chinese Physics B. 2020, In press], we explored a novel global-direction stencil and combined it with the face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Greatly inspired by the differential form, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization from integral form on both second and third-order finite volume solver. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has superiorities on both computational accuracy and convergence rate.


2020 ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

Abstract Accuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work [Chinese Physics B. 2020, In press] , we explored a novel global-direction stencil and combined it with the face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Greatly inspired by the differential form, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization from integral form on both second and third-order finite volume solver. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has superiorities on both computational accuracy and convergence rate.


2020 ◽  
Author(s):  
Lingfa Kong ◽  
Yidao Dong ◽  
Wei Liu ◽  
Huaibao Zhang

Abstract Accuracy of unstructured finite volume discretization is greatly influenced by the gradient reconstruction. For the commonly used k-exact reconstruction method, the cell centroid is always chosen as the reference point to formulate the reconstructed function. But in some practical problems, such as the boundary layer, cells in this area are always set with high aspect ratio to improve the local field resolution, and if geometric centroid is still utilized for the spatial discretization, the severe grid skewness cannot be avoided, which is adverse to the numerical performance of unstructured finite volume solver. In previous work, we explored a novel global-direction stencil and combine it with face-area-weighted centroid on unstructured finite volume methods from differential form to realize the skewness reduction and a better reflection of flow anisotropy. Note, however, that the differential form is hard to achieve higher-order accuracy, and in order to set stage for the method promotion on higher-order numerical simulation, in this research, we demonstrate that it is also feasible to extend this novel method to the unstructured finite volume discretization in integral form. Numerical examples governed by linear convective, Euler and Laplacian equations are utilized to examine the correctness as well as effectiveness of this extension. Compared with traditional vertex-neighbor and face-neighbor stencils based on the geometric centroid, the grid skewness is almost eliminated and computational accuracy as well as convergence rate is greatly improved by the global-direction stencil with face-area-weighted centroid. As a result, on unstructured finite volume discretization from integral form, the method also has a better numerical performance.


Author(s):  
Yaw Kyei

Efficient higher-order accurate finite volume schemes are developed for the threedimensional Poisson’s equation based on optimizations of an equation error expansion on local control volumes. A weighted quadrature of local compact fluxes and the flux integral form of the equation are utilized to formulate the local equation error expansions. Efficient quadrature weights for the schemes are then determined through a minimization of the error expansion for higher-order accurate discretizations of the equation. Consequently, the leading numerical viscosity coefficients are more accurately and completely determined to optimize the weight parameters for uniform higher-order convergence suitable for effective numerical modeling of physical phenomena. Effectiveness of the schemes are evaluated through the solution of the associated eigenvalue problem. Numerical results and analysis of the schemes demonstrate the effectiveness of the methodology.


Water ◽  
2021 ◽  
Vol 13 (11) ◽  
pp. 1597
Author(s):  
Susovan Pal ◽  
Prashanth Reddy Hanmaiahgari ◽  
Bryan W. Karney

Here, recent developments in the key numerical approaches to water hammer modelling are summarized and critiqued. This paper summarizes one-dimensional modelling using the finite difference method (FDM), the method of characteristics (MOC), and especially the more recent finite volume method (FVM). The discussion is briefly extended to two-dimensional modelling, as well as to computational fluid dynamics (CFD) approaches. Finite volume methods are of particular note, since they approximate the governing partial differential equations (PDEs) in a volume integral form, thus intrinsically conserving mass and momentum fluxes. Accuracy in transient modelling is particularly important in certain (typically more nuanced) applications, including fault (leakage and blockage) detection. The FVM, first advanced using Godunov’s scheme, is preferred in cases where wave celerity evolves over time (e.g., due to the release of air) or due to spatial changes (e.g., due to changes in wall thickness). Both numerical and experimental studies demonstrate that the first-order Godunov’s scheme compares favourably with the MOC in terms of accuracy and computational speed; with further advances in the FVM schemes, it progressively achieves faster and more accurate codes. The current range of numerical methods is discussed and illustrated, including highlighting both their limitations and their advantages.


2017 ◽  
Vol 22 (5) ◽  
pp. 1385-1412 ◽  
Author(s):  
Yilang Liu ◽  
Weiwei Zhang ◽  
Chunna Li

AbstractThis paper proposes a novel distance derivative weighted ENO (DDWENO) limiter based on fixed reconstruction stencil and applies it to the second- and highorder finite volume method on unstructured grids. We choose the standard deviation coefficients of the flow variables as the smooth indicators by using the k-exact reconstruction method, and obtain the limited derivatives of the flow variables by weighting all derivatives of each cell according to smoothness. Furthermore, an additional weighting coefficient related to distance is introduced to emphasize the contribution of the central cell in smooth regions. The developed limiter, combining the advantages of the slope limiters and WENO-type limiters, can achieve the similar effect of WENO schemes in the fixed stencil with high computational efficiency. The numerical cases demonstrate that the DDWENO limiter can preserve the numerical accuracy in smooth regions, and capture the shock waves clearly and steeply as well.


Sign in / Sign up

Export Citation Format

Share Document