Grid-Free Vortex Methods for Natural Convection in Two-Dimensional Domains

Author(s):  
Issam Lakkis

Vortex methods for simulating natural convection of an ideal gas in unbounded two-dimensional domains are presented. In particular, the redistribution method for diffusion is extended to enable simulation of nonlinear diffusion of an ideal gas in isobaric conditions encountered in unbounded low-Mach number flows. We also address the problem of handling source terms in grid-free vortex methods and propose a fast, accurate, and physically motivated method for solving the associated inverse problems. Examples include generation of baroclinic vorticity in non-reacting buoyancy driven flows, and in addition, generation of internal energy and species in buoyant reacting flows. Accuracy and speed of the proposed algorithms for nonlinear diffusion and vorticity generation are investigated separately. Simulations of natural convection of a “thermal patch” for Grashof number ranging from to 1562.5 to 25000 are presented.

Author(s):  
Leila Issa ◽  
Issam Lakkis

We present reduced order models of unsteady low Mach number isothermal ideal gas flows in two-dimensional rectangular microchannels subject to first order slip boundary conditions. The Navier-Stokes equations are simplified using Low Mach Number expansions of the pressure and velocity fields. This approximation allows decoupling the density from spatial pressure variations, thus simplifying the momentum equation. The resulting diffusion equation and the subsequent pressure-flow-rate relationship enables modeling the flow using analog circuit components. The accuracy of the proposed models is investigated for steady and unsteady flows in a two-dimensional channel for different values of Reynolds and Knudsen numbers.


Author(s):  
K. M. Akyuzlu ◽  
Y. Pavri ◽  
A. Antoniou

A two-dimensional, mathematical model is adopted to investigate the development of buoyancy driven circulation patterns and temperature contours inside a rectangular enclosure filled with a compressible fluid (Pr=1.0). One of the vertical walls of the enclosure is kept at a higher temperature then the opposing vertical wall. The top and the bottom of the enclosure are assumed insulated. The physics based mathematical model for this problem consists of conservation of mass, momentum (two-dimensional Navier-Stokes equations) and energy equations for the enclosed fluid subjected to appropriate boundary conditions. The working fluid is assumed to be compressible through a simple ideal gas relation. The governing equations are discretized using second order accurate central differencing for spatial derivatives and first order forward finite differencing for time derivatives where the computation domain is represented by a uniform orthogonal mesh. The resulting nonlinear equations are then linearized using Newton’s linearization method. The set of algebraic equations that result from this process are then put into a matrix form and solved using a Coupled Modified Strongly Implicit Procedure (CMSIP) for the unknowns (primitive variables) of the problem. A numerical experiment is carried out for a benchmark case (driven cavity flow) to verify the accuracy of the proposed solution procedure. Numerical experiments are then carried out using the proposed compressible flow model to simulate the development of the buoyancy driven circulation patterns for Rayleigh numbers between 103 and 105. Finally, an attempt is made to determine the effect of compressibility of the working fluid by comparing the results of the proposed model to that of models that use incompressible flow assumptions together with Boussinesq approximation.


2002 ◽  
Vol 124 (4) ◽  
pp. 791-799 ◽  
Author(s):  
Achintya Mukhopadhyay ◽  
Xiao Qin ◽  
Suresh K. Aggarwal ◽  
Ishwar K. Puri

A new formulation for extending the concept of heatlines and masslines to reacting flows through use of conserved scalars has been proposed. The formulation takes into account the distinct diffusion coefficients of different species. Results have been obtained for a number of two-dimensional nonreacting and reacting free shear flows under normal and zero gravity. For nonreacting flows, total enthalpy and elemental mass fractions have been used as the transported conserved scalars. For reacting flows, mixture fractions, defined as normalized elemental mass fractions and enthalpy, have been employed. The results show this concept to be a useful tool for obtaining better insights into the global qualitative picture of scalar transport for both nonreacting and reacting flows.


Sign in / Sign up

Export Citation Format

Share Document