Curvilinear Coordinate System for Mathematical Analysis of Inclined Buoyant Jets Using the Integral Method

• Main Authors: Aristeidis A. Bloutsos, Panayotis C. Yannopoulos
• Format: Article
• Language: English
• Published: Hindawi Limited 2018-01-01
• Series: Mathematical Problems in Engineering
• Online Access: http://dx.doi.org/10.1155/2018/3058425
• ISSN: 1024-123X; 1563-5147

The development of a local system of orthogonal curvilinear coordinates, which is appropriate to monitor the flow of an inclined buoyant jet with reference to the basic Cartesian coordinate system is presented. Such a system is necessary for the correct application of the integral method, since the well-known Gaussian profiles should be integrated on the cross-sectional area of inclined buoyant jet, where they are valid. This is the major advantage of the present work compared to all other integral methods using Cartesian coordinate systems. Consequently, the flow and mixing governing partial differential equations (PDE), i.e., continuity, momentum, buoyancy, and/or tracer conservation, are written in the local orthogonal curvilinear coordinate system and, then, the Reynolds substitution regarding mean and fluctuating components of all dependent variables is applied. After averaging with respect to time, the mean flow PDEs are taken, omitting second-order terms, as the dynamic pressure and molecular viscosity, compared to the mean flow and mixing contributions of turbulent terms. The latter are introduced through empirical coefficients. The Boussinesq’s approximation regarding small density differences is taken into consideration. The system of PDEs is closed by assuming known spreading coefficients along with Gaussian similarity profiles. The methodology is applied in the inclined two-dimensional buoyant jet; thus, PDEs are integrated on the jet cross-sectional area resulting in ordinary differential equations (ODE), which are appropriate to be solved by applying the 4th order Runge-Kutta algorithm coded in either FORTRAN or EXCEL. The numerical solution of ODEs, concerning trajectory of the inclined two-dimensional buoyant jet, as well as longitudinal variations of the mean axial velocity, mean concentration, minimum dilution, and entrainment velocity or entrainment coefficient, occurs quickly, saving computer memory and effort. The satisfactory agreement of results with experimental data available in the literature empowers the usefulness of the proposed methodology in inclined buoyant jets.