Search results

The lattice Boltzmann simulation of fluid flow in partial porous geometries with curved porous-fluid interfaces has not been investigated yet. It is mainly because of the lack of a method in the lattice Boltzmann framework to model the hydrodynamic compatibility conditions at curved porous-fluid interfaces, which is required for the two-domain approach. Therefore, the purpose of this study is to develop such a method.

This research extends the non-equilibrium extrapolation lattice Boltzmann method for satisfying no-slip conditions at curved solid boundaries, to model hydrodynamic compatibility conditions at curved porous-fluid interfaces.

The proposed method is tested against the results available from conventional numerical methods via the problem of fluid flow through and around a porous circular cylinder in crossflow. As such, streamlines, geometrical characteristics of recirculating wakes and drag coefficient are validated for different Reynolds (5 ≤ Re ≤ 40) and Darcy (10⁻⁵ ≤ Da ≤ 5 × 10⁻¹) numbers. It is also shown that without applying any compatibility conditions at the interface, the predicted flow structure is not satisfactory, even for a very fine mesh. This result highlights the importance of the two-domain approach for lattice Boltzmann simulation of the fluid flow in partial porous geometries with curved porous-fluid interfaces.

No research is found in the literature for applying the hydrodynamic compatibility conditions at curved porous-fluid interfaces in the lattice Boltzmann framework.

This study aims to investigate lattice Boltzmann (LB) simulation of the fluid flow and heat transfer characteristics of a heated porous elliptic cylinder in uniform flow based on the two-domain scheme. In the present research, the effect of axis ratio (1 ≤ AR ≤ 2), Reynolds number (5 ≤ Re ≤ 40) and Darcy number (10⁻⁴ ≤ Da ≤ 10⁻²) are studied.

To perform the LB simulation based on the two-domain scheme, the nonequilibrium extrapolation method is modified to model the heat transfer interfacial conditions required at the curved interface.

The results show that the axis ratio as well as Reynolds and Darcy numbers significantly affect the fluid flow and heat transfer characteristics of the porous elliptic cylinder. It is shown that for AR > 1, the phenomenon of detached recirculating zone occurs at much higher Darcy numbers compared with the case of the porous circular cylinder (AR = 1). The results show that the location of maximum temperature within the cylinder moves downstream when the Reynolds number, Darcy number and axis ratio increase. It is also concluded that the average Nusselt number of a porous elliptic cylinder is always lower than that of a porous circular cylinder.

The LB simulation of forced convection from a porous cylinder in uniform flow with a curved interface based on the two-domain scheme has not been studied yet.

The purpose of this paper is to implement the numerical analysis based on finite volume method to compare the effects of stress-jump (SJ) and stress-continuity (SC) conditions on flow structure around and through a porous circular cylinder.

In this study, a steady flow of a viscous, incompressible fluid around and through a porous circular cylinder of diameter “D,” using Darcy-Brinkman-Forchheimer’s equation in the porous region, is discussed. The SJ condition proposed by Ochoa-Tapia and Whitaker is applied at the porous-fluid interface and compared with the traditional interfacial condition based on the SC condition in fluid and porous media. Equations with the relevant boundary conditions are numerically solved using a finite volume approach. In this study, Reynolds and Darcy numbers are varied within the ranges of 1 < Re < 40 and 10-7 < Da < 10-2, respectively, and the porosities are e=0.45, 0.7 and 0.95.

Results show that the SJ condition leads to a much smaller boundary layer within porous medium near the interface as compared to the SC condition. Two interfacial conditions yield similar results with decrease in porosity.

There is no published research in the literature about the effects of important parameters, such as Porosity and Darcy numbers on different fluid-porous interface conditions for a porous cylinder and comparison the effects of SJ and SC conditions on flow structure around and through a porous circular cylinder.

The purpose of this paper is to develop a coupled lattice Boltzmann model for the simulation of the freezing process in unsaturated porous media.

In the developed model, the porous structure with complexity and disorder was generated by using a stochastic growth method, and then the Shan-Chen multiphase model and enthalpy-based phase change model were coupled by introducing a freezing interface force to describe the variation of phase interface. The pore size of porous media in freezing process was considered as an influential factor to phase transition temperature, and the variation of the interfacial force formed with phase change on the interface was described.

The larger porosity (0.2 and 0.8) will enlarge the unfrozen area from 42 mm to 70 mm, and the rest space of porous medium was occupied by the solid particles. The larger specific surface area (0.168 and 0.315) has a more fluctuated volume fraction distribution.

The concept of interfacial force was first introduced in the solid–liquid phase transition to describe the freezing process of frozen soil, enabling the formulation of a distribution equation based on enthalpy to depict the changes in the water film. The increased interfacial force serves to diminish ice formation and effectively absorb air during the freezing process. A greater surface area enhances the ability to counteract liquid migration.

This paper aims to consider numerical analysis of laminar double-diffusive natural convection inside a non-homogeneous closed medium composed of a saturated porous matrix and a clear binary fluid under spatial sinusoidal heating/cooling on one side wall and uniform salting.

The domain of interest is a partially square porous enclosure with sinusoidal wall heating and cooling. The fluid flow, heat and mass transfer dimensionless governing equations associated with the corresponding boundary conditions are discretized using the finite volume method. The resulting algebraic equations are solved by an in-house FORTRAN code and the SIMPLE algorithm to handle the non-linear character of conservation equations. The validity of the in-house FORTRAN code is checked by comparing the current results with previously published experimental and numerical works. The effect of the porous layer thickness, the spatial frequency of heating and cooling, the Darcy number, the Rayleigh number and the porous to fluid thermal conductivity ratio is analyzed.

The results demonstrate that for high values of the spatial frequency of heating and cooling (f = 7), temperature contours show periodic variations with positive and negative values providing higher temperature gradient near the thermally active wall. In this case, the temperature variation is mainly in the porous layer, while the temperature of the clear fluid region is practically the same as that imposed on the left vertical wall. This aspect can have a beneficial impact on thermal insulation. Besides, the porous to fluid thermal conductivity ratio, Rk, has practically no effect on Shhot wall, contrary to Nuinterface where a strong increase is observed as Rk is increased from 0.1 to 100, and much heat transfer from the hot wall to the clear fluid via the porous media is obtained.

The findings are useful for devices working on double-diffusive natural convection inside non-homogenous cavities.

The authors believe that the presented results are original and have not been published elsewhere.

The purpose of this study is to present a theoretical investigation of the effects of cavitation and couple stress on the squeeze film behavior between an anisotropic poroelastic rigid disc and a sinusoidally oscillating rigid disc.

Based on the microcontinuum theory of Vijay Kumar Stokes and the Elrod–Adam algorithm, the non-Newtonian Reynolds equation coupled with modified Darcy's law for lubricant flow through the porous disc is derived. This numerical study includes the continuity of tangential velocity at the porous–fluid interface and the effects of percolation of the polar additives into the anisotropic porous disc.

The effects of couple stress, oscillating amplitude, percolation additives, permeability and anisotropic permeability on the squeeze film characteristics are discussed. It is found that both the percolation effect of the lubricant additives and the anisotropic structure of the porous surface reduce the flow in the porous disc, resulting in a decrease in pressure. It is also observed that cavitation effects are more pronounced for Newtonian fluids than couple stress fluids.

The results of this study can be used to design a variety of engineering applications such as bearings, wet clutches and non-contact mechanical seals.

Thermally driven flow in a tall inclined cavity bounded by porous layers is studied analytically and numerically. A constant heat flux is applied for heating and cooling of two opposing walls of the cavity, while the other two are insulated. The Beavers—Joseph slip condition on velocity is applied at the interface between the fluid and porous layers. An analytical solution is obtained by assuming parallel flow in the core region of the cavity and a numerical solution by solving the complete governing equations. The flow and heat transfer variables are obtained in terms of the Rayleigh number, Ra, slip condition parameter N and angle of inclination of the cavity Φ. The critical Rayleigh numbers for the onset of convection in a layer heated from below are predicted for various hydrodynamic boundary conditions. The results for a fluid layer bounded by solid walls (N → ∞) and by free surfaces (N → 0) emerge from the present analysis as limiting cases.

The purpose of this paper is to numerically investigate the heat transfer enhancement in a tube filled partially with porous media under non-uniform porosity distribution and thermal dispersion effects. The optimum porous thickness ratio [R_(r,Nu)] for the heat transfer enhancement under these conditions with and without considering required pumping power is evaluated.

The local thermal non-equilibrium and Darcy–Brinkman–Forchheimer models are used to simulated thermal and flow fields in porous region. The tube wall and flow regime are assumed to be isothermal and laminar, respectively. The impacts of Darcy number (Da = 10^-6 - 10^-1) and inertia parameter (F = 0 − 2) on the Nusselt number and friction factor are studied for non-uniform porosity distribution.

First, the effect of Nusselt number indicates that there are two different behaviors with respect to uniform and non-uniform porosity for partially and fully filled porous pipe. Second, variable porosity in porous region has significant influence on the optimum thickness ratio with considering required pumping power. For negligible inertia term, it depends on the Darcy number, whereas it is 0.9 at F > 1. Third, the plug flow assumption cannot be valid even at lower Darcy number under non-uniform porosity, while this assumption is applicable at Da < 10^-3 for constant porosity distribution in porous region.

According to the best knowledge of authors, the optimum porous thickness ratio for the heat transfer enhancement considering the pressure loss effects under variable porosity has not reported up to now. Also the plug flow assumption in such physics is not discussed.

Our purpose was to theoretically analyse the effects of slip velocity and a magnetic fluid lubricant characterised by a material parameter on a parallel plate porous slider bearing. Various special cases were also discussed.

The analysis might be useful in designing better slider bearings useful in machines like gas turbines.

Computed values displayed in tabular form show that frictional force on the slider and the coefficient of friction decreased, the load capacity and the position of the centre of pressure unaltered when the slip parameter was increased. On the other hand, increase in the material parameter caused decrease in the frictional force and the coefficient of friction and, shift of the position of the centre of pressure towards the bearing inlet with no change in the load capacity of the bearing.

The parallel plate slider bearing cannot support a load with a conventional lubricant. However, it can support a load with a magnetic fluid lubricant with less frictional force on the slider.

The purpose of this paper is to examine the flow, heat transfer and entropy generation characteristics for an inclined channel of two immiscible micropolar fluids.

The flow region consists of two zones, the flow of the heavier fluid taking place in the lower zone. The flow is assumed to be governed by Eringen’s micropolar fluid flow equation. The resulting governing equations are then solved using the homotopy analysis method.

The following findings are concluded: first, the entropy generation rate is more near the plates in both the zones as compared to that of the interface. This indicates that the friction due to surface on the fluids increases entropy generation rate. Second, the entropy generation rate is more near the plate in Zone I than that of Zone II. This may be due to the fact that the fluid in Zone I is more viscous. This indicates the more the viscosity of the fluid is, the more the entropy generation. Third, Bejan number is the maximum at the interface of the fluids. This indicates that the amount of exergy (available energy) is maximum and irreversibility is minimized at the interface between the fluids. Fourth, as micropolarity increases, entropy generation rate near the plates decreases and irreversibility decreases. This indicates an important industrial application for micropolar fluids to use them as a good lubricant.

The problem is original as no work has been reported on entropy generation in an inclined channel with two immiscible micropolar fluids.