This research presents novel algorithms for computing flow within an unstructured, collocated, finite volume solver in the presence of non-orthogonality and compressibility in order to extend the range of problems which can be modelled with the University's in-house CFD code: PHYSICA.
A new non-orthogonality diffusion correction relaxation parameter has been successfully introduced and tested with benchmarks from the literature. Cases involving geometries meshed with commercial packages have been successfully run with the diffusion correction methods, variable bounding and proper under-relaxation practices. The applicability of a pressure interpolation method has also been tested with these cases.
A procedure for solving compressible flow within a finite volume, pressure correction type scheme, has been devised and successfully implemented in different test cases. This method is however prone to numerical diffusion in the presence of shocks, but does work even in the presence of skewed meshes. The method was then tested with the case of an oxygen jet entering a heated furnace, for which experimental data is available for comparison. The method was successful in predicting the axial variables of the jet, and used to develop a turbulence modification model for such jets.
The method was finally used to model the deformation of a free surface impinged by a compressible jet, using a novel zonal method called zonal Gas And Liquid Analyser (GALA). Convergence was achieved with the method developed in this research, together with the application of the counter diffusion method to model the moving interface
