Theses and Dissertations

Issuing Body

Mississippi State University


Walters, D. Keith

Committee Member

Janus, J. Mark

Committee Member

Kim, Seongjai

Committee Member

Luke, Edward A.

Date of Degree


Document Type

Dissertation - Open Access


Computational Engineering (Program)

Degree Name

Doctor of Philosophy


James Worth Bagley College of Engineering


Computational Engineering Program


A new implicit and compact optimization-based method is presented for high order derivative calculation for finite-volume numerical method on unstructured meshes. Highorder approaches to gradient calculation are often based on variants of the Least-Squares (L-S) method, an explicit method that requires a stencil large enough to accommodate the necessary variable information to calculate the derivatives. The new scheme proposed here is applicable for an arbitrary order of accuracy (demonstrated here up to 3rd order), and uses just the first level of face neighbors to compute all derivatives, thus reducing stencil size and avoiding stiffness in the calculation matrix. Preliminary results for a static variable field example and solution of a simple scalar transport (advection) equation show that the proposed method is able to deliver numerical accuracy equivalent to (or better than) the nominal order of accuracy for both 2nd and 3rd order schemes in the presence of a smoothly distributed variable field (i.e., in the absence of discontinuities). This new Optimization-based Gradient REconstruction (herein denoted OGRE) scheme produces, for the simple scalar transport test case, lower error and demands less computational time (for a given level of required precision) for a 3rd order scheme when compared to an equivalent L-S approach on a two-dimensional framework. For three-dimensional simulations, where the L-S scheme fails to obtain convergence without the help of limiters, the new scheme obtains stable convergence and also produces lower error solution when compared to a third order MUSCL scheme. Furthermore, spectral analysis of results from the advection equation shows that the new scheme is better able to accurately resolve high wave number modes, which demonstrates its potential to better solve problems presenting a wide spectrum of wavelengths, for example unsteady turbulent flow simulations.