Suhas V. Patankar

Bio: Suhas V. Patankar is an academic researcher from University of Minnesota. The author has contributed to research in topics: Heat transfer & Turbulence. The author has an hindex of 42, co-authored 135 publications receiving 7165 citations. Previous affiliations of Suhas V. Patankar include Innovative Research Inc..

Ray effect and false scattering in the discrete ordinates method

Abstract: A discussion on the ray effect and false scattering occurring in discrete ordinates solution of the radiative transfer equation is presented in this article. Ray effect arises from the approximation of a continuously varying angular nature of radiation by a specified set of discrete angular directions. It is independent of the spatial discretization practice. False scattering, on the other hand, is a consequence of the spatial discretization practice and is independent of the angular discretization practice. In multidimensional computations, when a beam is not aligned with the grid line, false scattering smears the radiative intensity field. It reduces the appearance of unwanted bumps, but does not eliminate ray effect. An inappropriate view of false scattering is also presented. Four sample problems are used to explain these two effects.

Analysis of melting in the presence of natural convection in the melt region

Abstract: An analysis of multidimensional melting is performed which takes account of natural convection induced by temperature differences in the liquid melt. Consideration is given to the melt region created by a heated vertical tube embedded in a solid which is at its fusion temperature. Solutions were obtained by an implicit finite-difference scheme tailored to take account of the movement of the liquid-solid interface as melting progresses. The results differed decisively from those corresponding to a conventional pure-conduction model of the melting problem. The calculated heat transfer rate at the tube wall decreased at early times and attained a minimum, then increased and achieved a maximum, and subsequently decreased. This is in contrast to the pure conduction solution whereby the heat transfer rate decreases monotonically with time. The thickness of the melt region was found to vary along the length of the tube, with the greatest thickness near the top. This contrasts with the uniform thickness predicted by the conduction solution. These findings indicate that natural convection effects, although unaccounted for in most phase change analyses, are of importance and have to be considered.

Numerical Prediction of Flow and Heat Transfer in a Parallel Plate Channel With Staggered Fins

Abstract: Fluid flow and heat transfer in two-dimensional finned passages were analyzed for constant property laminar flow. The passage is formed by two parallel plates to which fins are attached in a staggered fashion. Both the plates are maintained at a constant temperature. Streamwise periodic variation of the cross-sectional area causes the flow and temperature fields to repeat periodically after a certain developing length. Computations were performed for different values of the Reynolds number, the Prandtl number, geometric parameters, and the fin-conductance parameter. The fins were found to cause the flow to deflect significantly and impinge upon the opposite wall so as to increase the heat transfer significantly. However, the associated increase in pressure drop was an order of magnitude higher than the increase in heat transfer. Streamline patterns and local heat transfer results are presented in addition to the overall results.

Airflow and Cooling in a Data Center

Abstract: This paper deals with the distribution of airflow and the resulting cooling in a data center. First, the cooling challenge is described and the concept of a raised-floor data center is introduced. In this arrangement, cooling air is supplied through perforated tiles. The flow rates of the cooling air must meet the cooling requirements of the computer servers placed next to the tiles. These airflow rates are governed primarily by the pressure distribution under the raised floor. Thus, the key to modifying the flow rates is to influence the fiow field in the under-floor plenum. Computational fluid dynamics (CFD) is used to provide insight into various factors affecting the airflow distribution and the corresponding cooling. A number of ways of controlling the airflow distribution are explored. Then attention is turned to the above-floor space, where the focus is on preventing the hot air from entering the inlets of computer serves. Different strategies for doing this are considered. The paper includes a number of comparisons of measurements with the results of CFD simulations.

A control volume-based finite-element method for solving the navier-stokes equations using equal-order velocity-pressure interpolation

Abstract: A control volume-based finite-element method for solving the Navier-Stokes equations using equal-order velocity-pressure interpolation is presented. Unlike the unequal-order-type methods that compute pressure at much fewer grid points than velocity, the proposed equal-order method calculates velocity and pressure at all the grid points in the domain. The validity of the method is established by applying it to solve some test problems

Enhancements of the simple method for predicting incompressible fluid flows

Abstract: Variations of the SIMPLE method of Patankar and Spalding have been widely used over the past decade to obtain numerical solutions to problems involving incompressible flows. The present paper shows several modifications to the method which both simplify its implementation and reduce solution costs. The performances of SIMPLE, SIMPLER, and SIMPLEC (the present method) are compared for two recirculating flow problems. The paper is addressed to readers who already have experience with SIMPLE or its variants.

Numerical solution of saddle point problems

Abstract: Large linear systems of saddle point type arise in a wide variety of applications throughout computational science and engineering. Due to their indefiniteness and often poor spectral properties, such linear systems represent a significant challenge for solver developers. In recent years there has been a surge of interest in saddle point problems, and numerous solution techniques have been proposed for this type of system. The aim of this paper is to present and discuss a large selection of solution methods for linear systems in saddle point form, with an emphasis on iterative methods for large and sparse problems.

Error analysis and estimation for the finite volume method with applications to fluid flows

Abstract: The accuracy of numerical simulation algorithms is one of main concerns in modern Computational Fluid Dynamics. Development of new and more accurate mathematical models requires an insight into the problem of numerical errors. In order to construct an estimate of the solution error in Finite Volume calculations, it is first necessary to examine its sources. Discretisation errors can be divided into two groups: errors caused by the discretisation of the solution domain and equation discretisation errors. The first group includes insufficient mesh resolution, mesh skewness and non-orthogonality. In the case of the second order Finite Volume method, equation discretisation errors are represented through numerical diffusion. Numerical diffusion coefficients from the discretisation of the convection term and the temporal derivative are derived. In an attempt to reduce numerical diffusion from the convection term, a new stabilised and bounded second-order differencing scheme is proposed. Three new methods of error estimation are presented. The Direct Taylor Series Error estimate is based on the Taylor series truncation error analysis. It is set up to enable single-mesh single-run error estimation. The Moment Error estimate derives the solution error from the cell imbalance in higher moments of the solution. A suitable normalisation is used to estimate the error magnitude. The Residual Error estimate is based on the local inconsistency between face interpolation and volume integration. Extensions of the method to transient flows and the Local Residual Problem error estimate are also given. Finally, an automatic error-controlled adaptive mesh refinement algorithm is set up in order to automatically produce a solution of pre-determined accuracy. It uses mesh refinement and unrefinement to control the local error magnitude. The method is tested on several characteristic flow situations, ranging from incompressible to supersonic flows, for both steady-state and transient problems.