This work examines the feasibility of a novel high-order numerical method, which has been termed Flux Correction. It has been given this name because it "corrects" the flux terms of an established numerical method, cancelling various error terms in the fluxes and making the method higher-order. In this work, this change is made to a traditionally second-order finite volume Galerkin method. To accomplish this, higher-order gradients of solution variables, as well as gradients of the fluxes are introduced to the method. Gradients are computed using lagrange interpolations in a fashion reminiscent of Finite Element techniques. For the Euler Equations, Flux Correction is compared against Flux Reconstruction, a derivative of the high-order Discontinuous Galerkin and Spectral Difference methods, both of which are currently popular areas of research in high-order numerical methods. Flux Correction is found to compare favorably in terms of accuracy, and exceeds expectations for convergence rates. For the full Navier-Stokes Equations, the effect of curved elements on Flux Correction are examined. Flux Correction is found to react negatively to curved elements due to the gradient procedure's poor handling of high-aspect ratio elements.

A novel high-order finite volume scheme using flux correction methods in conjunction with structured finite difference schemes is extended to low Mach and incompressible flows on strand grids. Flux correction achieves high-order by explicitly canceling low-order truncation error terms in the finite volume cell. The flux correction method is applied in unstructured layers of the strand grid. The layers are then coupled together using a source term containing the derivatives in the strand direction. Proper source term discretization is verified. Strand-direction derivatives are obtained by using summation-by-parts operators for the first and second derivatives. A preconditioner is used to extend the method to low Mach and incompressible flows. We further extend the method to turbulent flows with the Spalart Allmaras model. We verify high-order accuracy via the method of manufactured solutions, method of exact solutions, and physical problems. Results obtained compare well to analytical solutions, numerical studies, and experimental data. It is foreseen that future application in the Naval field will be possible.