A case-study in open-source CFD code verification. Part II: Boundary condition non-orthogonal correction
Abstract This investigation constitutes the follow-up to a previous work in which we applied the manufactured solution procedure to verify the OpenFOAM diffusion operator on different types of meshes. Theoretical convergence orders were observed for Poisson’s equation solver on orthogonal hexahedral grids using several boundary conditions. We noticed that on non-orthogonal grids even slight mesh distortions reduce the theoretical second order convergence rate to first order. Our investigation showed that this loss in convergence order takes place when Dirichlet and/or Neumann boundary conditions are used on non-orthogonal meshes. In this paper, we introduce ways to achieve theoretical second order convergence accuracy, mainly based on applying non-orthogonal correction to Dirichlet and/or Neumann boundary condition schemes. A peculiarity was discovered for the Least-Squares gradient scheme that slightly affects second-order convergence, and which we attribute to the non-orthogonal correction at boundaries in the Least-Squares vector reconstruction library. The verification procedure is performed using the method of manufactured solutions for the Poisson equation and an analytical solution, the SIMPLE solver, for the Navier–Stokes equations. The average and maximum error norms were used to calculate the convergence rate. Comparative results are presented.
Verification of RANS solvers with manufactured solutions
This paper discusses code verification of Reynolds-Averaged Navier Stokes (RANS) solvers with the method of manufactured solutions (MMS). Examples of manufactured solutions (MSs) for a two-dimensional, steady, wall-bounded, incompressible, turbulent flow are presented including the specification of the turbulence quantities incorporated in several popular eddy-viscosity turbulence models. A wall-function approach for the MMS is also described. The flexiblity and usefulness of the MS is illustrated with calculations performed in three different exercises: the calculation of the flow field using the manufactured eddy-viscosity; the calculation of the eddy-viscosity using the manufactured velocity field; the calculation of the complete flow field coupling flow and turbulence variables. The results show that the numerical performance of the flow solvers is model dependent and that the solution of the complete problem may exhibit different orders of accuracy than in the exercises with no coupling between the flow and turbulence variables.
A new method of near field analysis
The problem of the accurate determination of the total electromagnetic field produced by a current distribution on a circular disc type aperture is examined. For comparison of the near and Fresnel region solutions, an exact Maxwellian integral is used as the starting point. Aperture and space coordinates are separated by the use of Bessel function addition theorems, allowing the field to be expressed as a sum of terms containing radial and angular integrals. The radial integrals involve only the disc (aperture) distribution and coordinates while the angular integrals involve only the elevation angle. The distance from origin appears implicitly in the series coefficients. These integrals are suitable for high speed computations and have been evaluated for a number of cases. The results are compared with the small angle and general Fresnel formula values for various observation angles, distances from the aperture, and aperture distributions.