Free Convection in Porous Media Introduction This model

  • Slides: 9
Download presentation
Free Convection in Porous Media

Free Convection in Porous Media

Introduction • This model exemplifies the use of COMSOL Multiphysics for modeling of free

Introduction • This model exemplifies the use of COMSOL Multiphysics for modeling of free convection in porous media. • It shows the following COMSOL Multiphysics features: – Porous media flow – Multiphysics between fluid flow and heat transfer – Results that are in excellent agreement with published models in the research journals in the field • The model has applications mainly in the fields of: – Geophysics – Chemical engineering

Geometry, Heating and Cooling Surfaces Tc Th-(Th-Tc)*s Th Tc Th-(Th-Tc)*s • Enclosed domain with

Geometry, Heating and Cooling Surfaces Tc Th-(Th-Tc)*s Th Tc Th-(Th-Tc)*s • Enclosed domain with porous material • The walls of the domain are impervious to flow • The walls are either heating or cooling surfaces with linear temperature profiles uniting the cool and hot surfaces • The arc length s goes from zero to 1 along a boundary segment.

Domain Equations Momentum and mass balances Boussinesq buoyant lifting term links flow and heat

Domain Equations Momentum and mass balances Boussinesq buoyant lifting term links flow and heat Heat balance

Brinkman equations for porous media flow p u h k r g b. T

Brinkman equations for porous media flow p u h k r g b. T T Tc Solution technique: Parametric solver to increase b. T from zero to problem- specific value = pressure = vector of directional velocities = dynamic viscosity = permeability = fluid density = gravity = thermal expansion coefficient = temperature from heat transfer application = initial temperature

Convection and conduction T Ke CL cp u = temperature = effective thermal conductivity

Convection and conduction T Ke CL cp u = temperature = effective thermal conductivity of fluid and solid medium = fluid volumetric heat capacity… CL= cp r = fluid specific heat capacity = vector of directional fluid velocities from flow application

Boundary Conditions • Brinkman equations no slip so velocity drops to zero at wall

Boundary Conditions • Brinkman equations no slip so velocity drops to zero at wall for unique solution fix pressure at a point • Convection and conduction

Results • Dimensionless temperature • Velocity field

Results • Dimensionless temperature • Velocity field

Concluding Remarks • • The model is simple to define and solve in COMSOL

Concluding Remarks • • The model is simple to define and solve in COMSOL Multiphysics The results give excellent agreement with published scientific papers, see M. Anwar Hossain and Mike Wilson, Natural convection flow in a fluid-saturated porous medium enclosed by non-isothermal walls with heat generation, International Journal of Thermal Sciences, Int. J. Therm. Sci. 41 (2002) 447 – 454.