A posterior error estimate for finite volume methods of the second order elliptic problem

2010 ◽  
Vol 27 (5) ◽  
pp. 1165-1178 ◽  
Author(s):  
Xiu Ye
2021 ◽  
Vol 2021 ◽  
pp. 1-6
Author(s):  
Xiong Liu ◽  
Wenming He

In this paper, we will investigate a multiscale homogenization theory for a second-order elliptic problem with rapidly oscillating periodic coefficients of the form ∂ / ∂ x i a i j x / ε , x ∂ u ε x / ∂ x j = f x . Noticing the fact that the classic homogenization theory presented by Oleinik has a high demand for the smoothness of the homogenization solution u 0 , we present a new estimate for the homogenization method under the weaker smoothness that homogenization solution u 0 satisfies than the classical homogenization theory needs.


2019 ◽  
Vol 40 (2) ◽  
pp. 1577-1600
Author(s):  
Gang Chen ◽  
Jintao Cui

Abstract Hybridizable discontinuous Galerkin (HDG) methods retain the main advantages of standard discontinuous Galerkin (DG) methods, including their flexibility in meshing, ease of design and implementation, ease of use within an $hp$-adaptive strategy and preservation of local conservation of physical quantities. Moreover, HDG methods can significantly reduce the number of degrees of freedom, resulting in a substantial reduction of computational cost. In this paper, we study an HDG method for the second-order elliptic problem with discontinuous coefficients. The numerical scheme is proposed on general polygonal and polyhedral meshes with specially designed stabilization parameters. Robust a priori and a posteriori error estimates are derived without a full elliptic regularity assumption. The proposed a posteriori error estimators are proved to be efficient and reliable without a quasi-monotonicity assumption on the diffusion coefficient.


Author(s):  
Nicole M. W. Poe ◽  
D. Keith Walters

Finite volume methods on structured and unstructured meshes often utilize second-order, upwind-biased linear reconstruction schemes to approximate the convective terms, in an attempt to improve accuracy over first-order methods. Limiters are employed to reduce the inherent variable over- and under-shoot of these schemes; however, they also can significantly increase the numerical dissipation of a solution. This paper presents a novel non-local, non-monotonic (NLNM) limiter developed by enforcing cell minima and maxima on dependent variable values projected to cell faces. The minimum and maximum values for a cell are determined primarily through the recursive reference to the minimum and maximum values of its upwind neighbors. The new limiter is implemented using the User Defined Function capability available in the commercial CFD solver Ansys FLUENT. Various simple test cases are presented which exhibit the NLNM limiter’s ability to eliminate non-physical oscillations while maintaining relatively low dissipation of the solution. Results from the new limiter are compared with those from other limited and unlimited second-order upwind (SOU) and first-order upwind (FOU) schemes. For the cases examined in the study, the NLNM limiter was found to improve accuracy without significantly increasing solution convergence rate.


Author(s):  
Nicole M. W. Poe ◽  
D. Keith Walters

Finite volume methods employing second-order gradient reconstruction schemes are often utilized to computationally solve the governing equations of transport. These reconstruction schemes, while not as dissipative as first-order schemes, frequently produce either dispersive or oscillatory solutions, especially in regions of discontinuities, and/or unsatisfactory levels of dissipation in smooth regions of the variable field. A novel gradient reconstruction scheme is presented in this work which shows significant improvement over traditional second-order schemes. This Optimization-based Gradient REconstruction (OGRE) scheme works to minimize an objective function based on the mismatch between local reconstructions at midpoints between cell stencil neighbors, i.e. the degree to which the projected values of a dependent variable and its gradients in a given cell differ from each of these values in neighbor cells. An adjustable weighting parameter is included in the definition of the objective function that allows the scheme to be tuned towards greater accuracy or greater stability. This scheme is implemented using the User Defined Function capability available in the commercially available CFD solver, Ansys FLUENT. Various test cases are presented that demonstrate the ability of the new method to calculate superior predictions of both a scalar transported variable and its gradients. These cases include calculation of a discontinuous variable field, several sinusoidal variable fields and a non-uniform velocity field. Results for each case are determined on both structured and unstructured meshes, and the scheme is compared with existing standard first- and second-order upwind discretization methods.


Sign in / Sign up

Export Citation Format

Share Document