Effects of Prandtl and Jacob Numbers and Dimensionless Thermal Conductivity on Velocity Profiles in Media (Porous and Liquid)

Abstract

In this present study, we analyzed the effects of Prandtl and Jacob numbers and dimensionless thermal conductivity on the velocity profiles in media (porous and liquid). The transfers in the porous medium and the liquid film are described respectively by the improved Wooding model and the classical boundary layer equations. The mesh of the digital domain is considered uniform in the transverse and longitudinal directions. The advection and diffusion terms are discretized with a back-centered and centered scheme respectively. The coupled systems of algebraic equations thus obtained are solved numerically using an iterative line-by-line relaxation method of the Gauss-Seidel type. The results show that the parameters relating to the thermal problem (the dimensionless thermal conductivity, the Prandtl (Pr) and Jacob (Ja) numbers) have no influence on the dimensionless speed, although the thermal and hydrodynamic problems are coupled. Via the heat balance equation. The results obtained show that the parameters relating to the thermal problem have no influence on the dimensionless speed, although the thermal and hydrodynamic problems are coupled via the heat balance equation. So, at first approximation with the chosen constants, we can solve the hydrodynamic problem independently of the thermal problem.

Share and Cite:

Ndiaye, M. , Sene, M. and Ndiaye, G. (2024) Effects of Prandtl and Jacob Numbers and Dimensionless Thermal Conductivity on Velocity Profiles in Media (Porous and Liquid). Journal of Modern Physics, 15, 1603-1614. doi: 10.4236/jmp.2024.1510068.

1. Introduction

The study of condensation in porous media is a topic of growing interest due to its importance in many technological fields, such as heat exchangers, fuel cells, geothermal energy, energy storage, etc. …It also finds its application in certain technological problems presenting a two-phase zone, called a pasty zone, represented by a saturated porous medium and a binary fluid. Numerous studies have been carried out using: digital experimentation (numerical simulation) and/or practical experimentation (laboratory experiment), as well as theoretical developments. We plan to study the effects of Prandtl and Jacob numbers and dimensionless thermal conductivity on velocity profiles in liquid and porous media.

NDIAYE M. et al. [1] proposed a model for the numerical study of the condensation of pure saturated steam of the thin film type in forced convection on a wall covered with porous material. The transfers in the porous medium and the liquid film are described respectively by the Darcy-Brinkman model and the classical equations of the boundary layer. The dimensionless equations are solved by an implicit finite difference method combined with an iterative Gauss-Seidel method. They analyzed the influences of Prandtl and Froude numbers on transfers in the liquid phase. The parameters relating to the thermal problem (Prandtl number) have no influence on the dimensionless speed although the thermal and hydrodynamic problems are coupled via the heat balance equation.

The study showed that heat transfer increases when Pr ≤ 0.2, however, for values Pr > 0.2, there is no change in heat transfer. However, the values of FrK ≥ 0.001 have no effect on the hydrodynamic field, for FrK < 0.001, the speed increases significantly in the liquid film and the effects of inertia can no longer be ignored (a reduction in the Froude number leads to an increase in the permeability K of the porous medium, thus the substance becomes more and more permeable).

NDIAYE M. et al. [2] also proposed a model for the numerical study of the condensation of pure saturated steam of the thin film type in forced convection on a wall covered with porous material and they analyzed the influences of the Reynolds and Jacob numbers, the thickness dimensions of the porous layer and the dimensionless heat transfer conductivity on transfers in the liquid phase. They used the same resolutions as before. Their results showed that the thermal field is not influenced by dynamic problems despite their coupling at the interface. The parameters linked to the thermal problem (thermal conductivity ratio and Jacob number) have no effect on the longitudinal speed. It also appears that the Reynolds number has little influence on the thermal field. The longitudinal speed increases with the thickness of the porous layer. The decrease in the dimensionless thickness of the porous layer and the thermal conductivity ratio and an increase in the Jacob number leads to an increase in temperature.

NDIAYE M. et al. [3], using the same numerical model, analyzed the influence of Reynolds numbers and the dimensionless thickness of the porous layer on transfers in the liquid phase and the porous medium as well as the thickness of the liquid film. The thickness of the condensate film is determined by the heat balance equation at the liquid-vapor interface which is solved using an iterative Gauss-Seidel type procedure. Their results showed that the Reynolds number does not influence the thermal field, the temperature increases with the decrease in the thickness of the dimensionless porous layer. The increase in the Reynolds number and the dimensionless thickness of the porous layer leads to an increase in the longitudinal speed and a reduction in the thickness of the liquid film. Increases in Reynolds numbers and dimensionless thickness of the porous layer have adverse effects on condensation (decrease in liquid film thickness). The study also showed that inertia effects can no longer be ignored when the Reynolds number calculated from the permeability coefficient is greater than 7.

PT Ndiaye et al. [4] [5] presented the numerical study of the condensation in thin layers in forced convection of a pure saturated vapor in a channel whose walls are covered with a porous material. The generalized Darcy-Brinkman-Forchheimer (DBF) model is used to describe the flow in the porous medium while the classical boundary layer equations have been exploited in the case of a pure liquid. The dimensionless equations are solved by an implicit finite difference method and the Gauss-Seidel iterative method. This study made it possible to examine and highlight the role of parameters such as the Reynolds number and the Prandtl number on the longitudinal speed and on the temperature in the porous medium and in the pure liquid, and the transfer rate heat (local Nusselt number). The increase in Reynolds and Prandtl numbers’ results in an increase in longitudinal velocity and temperature. The tangents of the velocity curves at the porous interface on the medium side are smaller than those obtained on the liquid side with low Reynolds and Prandtl numbers. We also note that increasing the Reynolds and Prandtl numbers improves the thermal performance of the interface. The Prandtl number improves heat exchanges at the interface of the porous medium and the liquid film.

Dharma Rao et al. [6] developed a model for the study of water vapor condensation in the presence of humid air by forced laminar convection in a vertical tube. The model is developed by considering the classical boundary layer equations but neglecting the convective terms in the liquid phase equations. The temperature at the gas-liquid interface and the condensate thickness are calculated by considering the thermal and mass balances at the interface. The equations are solved numerically by the finite difference method. Calculations were carried out for a wide range of parameters to determine the influence of fluid inlet conditions on transfers. The authors show that the Nusselt and liquid Reynolds numbers increase with temperature, relative humidity and the Reynolds number of the gas mixture at the tube inlet, but decrease as the mixture inlet pressure increases. However, the effects of these parameters on the Nusselt number of the gas phase are less significant.

Maheshwari et al. [7] developed a theoretical model to determine the local heat transfer coefficient during vapor condensation in the presence of a non-condensable gas, flowing in a vertical tube for a wide range of the Reynolds number of the mixture. The equations for transfer in the mixture were developed using the analogy relationships between heat and mass transfer. The flow of the liquid film is governed by the Nusselt model modified to take into account the undulatory character 26 of the film at the interface. The results show that the thermal resistance of the layer of non-condensable gases at the mixture-liquid interface is greater than that of the condensate, at low intake Reynolds number of the mixture. When this increases, the thermal resistance of the liquid film becomes preponderant.

Asbik M. et al. [8]-[10] were interested in a semi analytical and numerical study of condensation by forced convection laminar film on a vertical plate covered by a porous layer. These authors showed the significant contribution of the porous substrate and the thermal effect of the dispersion on a flat wall causing an increase in the transfer rate.

Chaynane R. et al. [11] Also considered by the condensation of a forced convection laminar sublayer of a pure on a saturated porous plate is inclined relative to the vertical steam. The Darcy-Brinkman model is used to describe the flow in the porous medium, while the classical equations of the boundary layer have been used in the liquid. The problem was solved by analytical and numerical means. Results are mainly presented as the dimensionless thickness of the liquid film, profiles of speed, temperature and heat transfer coefficients represented by the Nusselt number. The results obtained were compared with those of experimental Renken et al. The effects of various influencing parameters such as the angle, the effective viscosity, Reynolds number, the dimensionless thickness of the porous substrate and the H * dimensionless thermal conductivity, the heat transfer and flow are illustrated.

K. Renken J. et al. [12]-[15] investigated the condensation laminar thin film on a vertical surface with a porous coating. The model simulates the condensation two dimensions within a very thin and very highly permeable porous conductive layer. Model Darcy-Brinkman-Forchheimer is used to describe the flow field in the porous layer while the classical equations of the boundary layer are used in the region of pure condensate. Numerical results, which describe in detail the dependence of the rate of heat transfer of the temperature field and the parameters (e.g., the Reynolds number, and the Prandtl Darcy, the thermal dispersion coefficient and the thickness the porous coating and the rate of thermal conductivity) are calculated using a finite difference scheme.

They also developed experiments on condensation in forced thin, porous coatings convection. This study presents the results of experiments of heat transfer by forced convection and condensation on plates with a thin porous layer. The composite system consists of a very thin porous conductive and permeable material, adhered to a cold condensation surface isotherm moves parallel to the flow of saturated steam. The results of this study provide significant predictions of the heat transfer of condensation that can be used in a number of thermal engineering applications that require improved heat transfer.

2. Mathematical Formulation

2.1. Physical Model

The research carried out in this article generally focuses on the modeling of flows and transfers in porous and liquid media. To carry out a simple formulation of the mathematical model describing the phenomenon of thin film condensation, we consider a saturated porous medium confined on a vertical plate, of thickness H, permeability K and porosity ε (Figure 1). This vertical flat plate of length L is placed in a flow of pure and saturated steam, with longitudinal speed U0. The steam condenses on the wall of the plate maintained at the temperature (Tw) lower than that of saturation (Ts) of the steam. The condensate film flows under the effect of gravity and viscous friction forces. There are three zones: Zone (1) is the porous medium saturated by the liquid. Zone (2) corresponds to the liquid film while zone (3) relates to the saturated vapor. Let (x, y) and (u, v) respectively be the Cartesian coordinates and the components of the speed in the porous medium and of the liquid in the reference frame associated with the model.

Figure 1. Physical model geometry and coordinate system.

2.2. Assumptions

We accept the following simplifying assumptions:

The porous matrix is isotropic and homogeneous;

The fluid saturating the porous medium is Newtonian and incompressible;

The flow generated is laminar and two-dimensional;

The work, induced by the viscous and pressure forces, is negligible;

The thermo-physical properties of the fluids and those of the porous matrix are assumed to be constant;

The improved Wooding model is used to describe the flow in the porous layer;

The effective dynamic and kinematic viscosities of the porous material are equal to those of the condensate film;

Condensation occurs as a thin film;

The porous matrix is in local equilibrium with the condensate;

The liquid-vapor interface is in thermodynamic equilibrium and the shear stress is assumed to be negligible;

Steam and film are separated by a distinct boundary;

The transverse and longitudinal variation of pressure is not taken into account.

2.3. Nondimensionalization and Transformation Equations

The equations governing transfers in the fields (1) and (2) defined above and the boundary conditions associated with them, were dimensionless using the following variables and parameters (Abik M et al. [2003]), we ask:

H = H K (1-a)

x = x K (1-b)

y = y K (1-c)

θ ξ = T ξ T w T s T w (2-a)

u = u ξ u r (2-b)

u r = K ν eff g (2-c)

δ = δ K (2-d)

To pinpoint the liquid-vapor interface and avoid the problems caused by non-uniformity of the mesh in the vicinity of the interface we make the following change of variable:

X= x ; η=coef y H +( 1coef ){ 1+ y H δ H } (3)

with coef = 1 if we are in the porous layer ( 0 y H ) and coef = 0 when we are in the pure liquid ( H y δ ).

coef = coefficient

So the plan ( x , y ) is transformed into a rectangular domain ( X,η ) .

Porous layer: 0η1

u p X + 1 H v p η =0 (4)

u p u p X + v p H u p η = ε 2 R e K u p + ε 2 R e K ( 1 H 2 2 u p η 2 )+ ε 2 F r K (5)

u p θ p X + v p H θ P η = 1 p r R e K H 2 2 θ P η 2 (6)

Pure liquid: 1η2

( δ H ) u l X ( η1 ) d δ dX u l η + v l η =0 (7)

u l [ u l X η1 δ H d δ dX u l η ]+ v l δ H u l η = 1 R e K [ 1+ μ ( δ H ) 2 δ 2 u l δ η 2 ] (8)

u l [ θ l X η1 δ H d δ dX θ l η ]+ v l δ H θ l η = 1 R e K P r ( δ H ) 2 2 θ l η 2 (9)

Boundary conditions:

in the wall η=0

u p = v p =0 (10-a)

θ p =0 (10-b)

interface to the porous layer/pure liquid, η=1

θ l = θ p (11-a)

u l = u p (11-b)

μ H u p η = 1 δ H u l η (12-a)

1 H θ p η = λ δ H θ l η (12-b)

A liquid interface/steam,

θ l =1 (13-a)

u l η =0 (13-b)

The heat and mass balance satisfy the following expressions:

Ja ( Pe ) eff 1 H * θ p η | η=0 = d d x * [ H * 0 1 { 1+Ja( 1 θ p ) } u p * dη ] + d d x * [ 1 2 ( δ * H * ){ 1+Ja( 1 θ l ) } u l * dη ] (14)

with

( Pe ) eff = λ Pr R e K (15)

The mass flow rate:

H 0 1 u p dη +( δ H ) 1 2 u l dη = ρ v ρ l δ (16)

3. Digital Methodology

The transfer equations are discretized by an implicit finite difference method. The mesh of the digital domain is considered uniform in the transverse and longitudinal directions. The advection and diffusion terms are discretized with a back-centered and centered scheme respectively. The coupled systems of algebraic equations thus obtained are solved numerically using an iterative line-by-line relaxation method of the Gauss-Seidel type.

The equations are discretized transfer by an implicit finite difference method. The mesh of the domain is considered uniform in the transverse and longitudinal directions. The terms of advection and diffusion are discretized respectively with a rear and centered upwind scheme.

Thus, the discretization in the field of study of the equations of continuity, energy and movement leads to the following algebraic equations:

b.In t ( i,j )=amiIn t ( i1,j )+amjIn t ( i,j1 ) +apjIn t ( i,j+1 )+coef0 (17)

(i, j) and (Δx, Δη) respectively represent nodes and the step along x and η.

im and jm are the maximum nodes along x and η

2iimet2jjm1 (18)

with

b= coefx Δx + coefe Δη +2 coef2 Δ η 2 +coefint ; (19-a)

ami= coefx Δx ; (19-b)

amj= coefe Δη + coef2 Δ η 2 (19-c)

apj= coef2 Δ η 2 (20-a)

coef2= ε 2 F r K (20-b)

coef1= ε 2 / R e K (20-c)

In t ( i,j )=In t ( i,j1 )+co1co2 (21)

Discretization of the equation of heat balance:

δ ( i )= δ ( i1 )+ R δ (22)

Coefficients that fall within the terms of the general equation are shown in the appendix.

Coupled algebraic systems thus obtained are solved numerically using an iterative relaxation method line by line Gauss-Seidel.

4. Results and Discussion

Our study makes it possible to examine and highlight the role of parameters such as: Prandtl numbers, Jacob numbers and dimensionless thermal conductivity on the velocity profiles in media (porous and liquid). The results from the numerical simulations relate to:

ν = 10 7 , μ =1 ε=0.4 ( ν dimensionless kinematic viscosity, μ =1 dimensionless dynamic viscosity, ε = porosity). The study of the sensitivity of the mesh led us to choose im = 51 et jm = 101 (im horizontal sensitivity and jm vertical sensitivity). The convergence criterion is set at 106.

To validate our model, we compared the results of Asbik et al. [8]-[11] with those from our calculation code in which the inertia terms have been omitted. We find that the correspondence is acceptable.

Figures 2-4 show that the parameters relating to the thermal problem (the dimensionless thermal conductivity, the Prandtl and Jacob numbers) have no influence on the dimensionless speed although the thermal and hydrodynamic problems are coupled via the heat balance equation. So at first approximation with the chosen constants, we can solve the hydrodynamic problem independently of the thermal problem.

Figure 2. Variation of the longitudinal speed to the ordinate η for different values of λ* ReK = 45, H* = 2.10−3, FrK = 104, Ja = 10−3, Pr = 2.

Figure 3. Variation of the longitudinal speed to the ordinate η for different values of Pr ReK = 45, H* = 2.10−3, FrK = 10−4, Ja = 10−3, Pr = 2.9.

Figure 3. Variation of the longitudinal speed to the ordinate η for different values of Ja ReK = 45, H* = 2.10−3, FrK = 10−4, Ja = 10−3, Pr = 2.9.

5. Conclusion

We studied the forced convection condensation of pure and saturated vapor on a vertical wall covered with porous material. The equations were solved using the back-centered implicit finite difference method for the velocity profiles in the two media (porous and liquid). We analyzed the influence of Prandtl and Jacob numbers and dimensionless thermal conductivity on the velocity profiles in the two media (porous and liquid). The results obtained show that the parameters relating to the thermal problem have no influence on the dimensionless speed although the thermal and hydrodynamic problems are coupled via the heat balance equation. So at first approximation with the chosen constants, we can solve the hydrodynamic problem independently of the thermal problem.

Nomenclature

Greek symbols:

δ

thickness of the condensate, m

ε

porosity

θ

temperature dimensionless

λ

thermal conductivity, W·m−1·K−1

μ

dynamic viscosity, kg·m−1·s−1

ν

kinematic viscosity, m2·s−1

ρ

density, kg/m3

Indices exhibitor:

eff

efficiency value

i

porous substrate interface/pure liquid

l

liquid

p

porous

s

saturation

v

steam

w

wall

*

dimensionless quantity

Latin letters:

c p

specific heat, J. kg−1·K−1

F r K

Froude number based on K

g

acceleration of gravity, m·s−2

H

thickness of the porous layer, m

h fg

heat of evaporation, J·kg−1

Ja

Jacob number

K

hydraulic conductivity or permeability, m2

L

length of the plate, m

Pe

Peclet number

Pr

Prandtl number

R e K

Reynolds number based on K

T

temperature, K

U 0

velocity of free fluid (steam), m/s

u

velocity along x, m/s

u r

velocity reference, m/s

v

velocity along y, m/s

x, y

cartesian coordinates along x and y, m

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Ndiaye, M., Mbow, C., Sarr, J., Zeghmati, B. and Faye, M. (2013) Numerical Investigation of Laminar Forced Thin Film Condensation of a Saturated Vapor along a Vertical Wall Covered with a Porous Material: Effect of Prandtl and Froude Numbers. International Journal on Heat and Mass Transfer Theory and Applications (IREHEAT), 1, 339-344.
[2] Ndiaye, M., Mbow, C., Sarr, J. and Zeghmati, B. (2013) Numerical Study of the Thin Film-Type Condensation of Saturated Forced into a Vertical Wall Covered with a Porous Material Vapor Convection. International Journal on Heat and Mass Transfer Theory and Applications (IREHEAT), 1, 330-338.
[3] Ndiaye, M., Mbow, C. and Sarr, J. (2014) Numerical Study of Thin Film Type Condensation in Forced Convection of a Saturated Vapor on a Vertical Wall Covered with a Porous Material. 3rd French-Speaking International Conference on Energy and Mechanics, Renewable Energies and Mechanics Applied to Industry, Moroni, 5-7 May 2014, 211-216.
[4] Ndiaye, P.T., Ndiaye, M., Mbow, C. and Ndiaye, G. (2020) Numerical Study of Thin Film Condensation in Forced Convection in a Canal Whose Walls Are Covered with a Porous Material: Influence of Jacob Number-Determination of Lengths of Entry. International Journal on Engineering Applications (IREA), 8, 125-132.
https://doi.org/10.15866/irea.v8i4.18687
[5] Ndiaye, P.T., Ndiaye, M., Mbow, C. and Ndiaye, G. (2020) Influence of Reynolds and Prandtl Numbers on Thin Film Condensation in Forced Convection in a Canal Covered with a Porous Material. International Journal on Engineering Applications (IREA), 8, 178-187.
https://doi.org/10.15866/irea.v8i5.18678
[6] Dharma Rao, V., Murali Krishna, V., Sharma, K.V. and Rao, P.V.J.M. (2008) Convective Condensation of Vapor in the Presence of a Non-Condensable Gas of High Concentration in Laminar Flow in a Vertical Pipe. International Journal of Heat and Mass Transfer, 51, 6090-6101.
https://doi.org/10.1016/j.ijheatmasstransfer.2008.03.027
[7] Maheshwari, N.K., Saha, D., Sinha, R.K. and Aritomi, M. (2004) Investigation on Condensation in Presence of a Noncondensable Gas for a Wide Range of Reynolds Number. Nuclear Engineering and Design, 227, 219-238.
https://doi.org/10.1016/j.nucengdes.2003.10.003
[8] Asbik, M., Zeghmati, B., Chaynane, R. and Bresson, J. (2002) Analytical Study of Laminar Film Condensation in Forced Convection of a Pure and Saturated Vapor on a Vertical Porous Wall: Effect of Thermal Dispersion. Proceedings of the SFT 2002 Congress, Vittel, 3-6 June 2002, 381-386.
[9] Asbik, M., Chaynane, R., Boushaba, H., Zeghmati, B. and Khmou, A. (2003) Analytical Investigation of Forced Convection Film Condensation on a Vertical Porous-Layer Coated Surface. Heat and Mass Transfer, 40, 143-155.
https://doi.org/10.1007/s00231-002-0406-8
[10] Asbik, M., Zeghmati, B., Louahlia-Gualous, H. and Yan, W.M. (2006) The Effect of Thermal Dispersion on Free Convection Film Condensation on a Vertical Plate with a Thin Porous Layer. Transport in Porous Media, 67, 335-352.
https://doi.org/10.1007/s11242-006-9028-9
[11] Chaynane, R., Asbik, M., Boushaba, H., Zeghmati, B. and Khmou, A. (2001) Study of Laminar Film Condensation in Forced Convection of a Pure and Saturated Vapor on the Porous Wall of an Inclined Plate. Mechanics & Industry, 5, 381-391.
https://doi.org/10.1051/meca:2004038
[12] Renken, K.J., Soltykiewicz, D.J. and Poulikakos, D. (1989) A Study of Laminar Film Condensation on a Vertical Surface with a Porous Coating. International Communications in Heat and Mass Transfer, 16, 181-192.
https://doi.org/10.1016/0735-1933(89)90019-5
[13] Renken, K.J. and Aboye, M. (1993) Analysis of Film Condensation within Inclined Thin Porous-Layer Coated Surfaces. International Journal of Heat and Fluid Flow, 14, 48-53.
https://doi.org/10.1016/0142-727x(93)90039-p
[14] Renken, K.J. and Meechan, K. (1995) Impact of Thermal Dispersion during Forced Convection Condensation in a Thin Porous/Fluid Composite System. Chemical Engineering Communications, 131, 189-205.
https://doi.org/10.1080/00986449508936291
[15] Renken, K.J. and Raich, M.R. (1996) Forced Convection Steam Condensation Experiments within Thin Porous Coatings. International Journal of Heat and Mass Transfer, 39, 2937-2945.
https://doi.org/10.1016/0017-9310(95)00367-3

Copyright © 2025 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.