Numerical modelling of viscous fluid flow in diverse porous media under transverse magnetic fields

Introduction

The study highlights the importance of multi-fluid transport in engineering and geophysics, particularly in petroleum recovery, MHD power generators, oil pipelines, desalination systems, and plasma control devices. Transverse magnetic fields, with or without heat transfer, significantly influence MHD flows. Previous research by Khedr et al. (2009) [3] examined MHD micropolar fluid flow along a vertical permeable plate, while Chamkha (1997) [3] and Magyari & Chamkha (2008) [5] explored MHD convection flows in various setups. The study explores the hydrodynamics of two immiscible fluids in a dual porous medium with oscillatory wall suction, inspired by technologies like magneto-hydrodynamic oil spill control. It develops and non-dimensionalizes laminar unsteady flow equations, solving them numerically using the Galerkin finite element method. The research analyzes the effects of key parameters, including Darcy and Forchheimer numbers, fluid viscosity ratio, Hartmann number, and Reynolds number, providing insights for optimized fluid flow control.

Background

The study of fluid dynamics in porous media is essential for applications in chemical engineering, environmental science, and geophysics. Darcy's law, extended by Brinkman (1947) to include viscous shear effects, provides a more detailed description of fluid flow. However, existing models often struggle to accurately represent flow in heterogeneous and anisotropic porous media, especially when thermal conductivity, mass transfer, and magnetic fields are considered.

Problem Statement

Current mathematical models for viscous fluid flow in porous media often assume homogeneity and isotropy, limiting their real-world applicability. They inadequately represent the coupled effects of heat, mass transfer, and magnetic fields, which significantly influence fluid dynamics. There is a need for more comprehensive models to improve accuracy and understanding of fluid behavior in complex porous structures.

Objectives of the Study

The study aims to develop and validate mathematical models for viscous fluid flow in

Impact of Geometrical Configurations on Fluid Dynamics

The geometrical configuration of a porous medium influence's fluid flow dynamics. Two-dimensional models simplify analysis but may overlook complexities, while three-dimensional models offer more accuracy at the cost of higher computational demands. Thermal Conductivity and Mass Transfer in Porous Media: Thermal conductivity and mass transfer significantly impact fluid flow in porous media. Thermal dispersion, driven by temperature gradients, influences fluid movement in applications like geothermal energy and chemical reactors. Mass transfer, including diffusion and convection, complicates flow dynamics and is essential for accurately predicting fluid behavior in processes like contaminant transport and catalytic reactions.

Role of Magnetic Fields in Fluid Dynamics

Magnetic fields can significantly influence fluid flow in porous media, especially in MHD applications, altering flow patterns depending on the field's strength and orientation. This is important for processes like electromagnetic filtration and enhanced oil recovery with magnetic nanoparticles. The interaction is described by MHD equations, combining Maxwell's electromagnetism equations with Navier-Stokes fluid dynamics, offering a more comprehensive understanding of fluid behavior under magnetic influences.

Summary and Identification of Research Gaps

Despite advancements, existing models of viscous fluid flow in porous media often assume homogeneity and isotropy, limiting their applicability. They also inadequately account for thermal conductivity, mass transfer, and magnetic field effects, leading to less accurate predictions. Future research should focus on developing comprehensive models that integrate these factors and use advanced numerical techniques, with experimental validation to improve reliability. Addressing these gaps will enhance predictive capabilities and support more efficient, sustainable industrial processes.

Mathematical Formulation of the Problem

The study investigates viscous fluid flow in a porous medium under a transverse magnetic field, focusing on fully developed plane Poiseuille flow of an incompressible, electrically conducting fluid through a channel with variable permeability, with its graphical representation in Figure 1(b). The study models fluid flow in a channel with stationary plates, driven by a constant pressure gradient. The channel has variable permeability, k=k(y)k = k(y)k=k(y), and is subject to a uniform transverse magnetic field, assuming negligible induced current and a small Magnetic Reynolds number.

Where 𝑢𝑢 is the velocity vector, 𝜇𝜇𝑒𝑒𝑓𝑓𝑓𝑓 is the effective viscosity, (𝑦𝑦) is the permeability, 𝐽𝐽 is the current density, and 𝐵𝐵 is the magnetic induction vector. Assuming the absence of an external electric field and negligible internal effects such as charge separation or polarization, the current density is given by:

Where 𝜎𝜎 is the electrical conductivity of the fluid. Consequently, the Lorentz force 𝐹𝐹𝐿𝐿𝑜𝑜𝑟𝑟𝑒𝑒𝑛𝑛𝑡𝑡𝑧𝑧 and velocity 𝑢𝑢 are collinear and opposite in direction, hence:

Thus, the governing equation simplifies to:

Boundary Conditions

For the plane Poiseuille flow, the boundary conditions are: (0) = 0 and (1) = 0 Numerical methods, like the Galerkin method, are used to solve the differential equation with boundary conditions, yielding the velocity distribution across the channel. This helps analyze how permeability, magnetic field strength, and other parameters affect fluid flow in porous media. Where 𝑝𝑝 is the pressure, 𝜇𝜇 is the dynamic viscosity, 𝑔𝑔 is the gravitational acceleration, and 𝐹𝐹𝑒𝑒𝑥𝑥𝑡𝑡 represents external forces, such as magnetic forces in magneto-hydrodynamic (MHD) flows. where 𝑘𝑘 is the permeability of the porous medium.

Incorporation of Heat and Mass Transfer

Heat and mass transfer affect fluid flow in porous media. Figure 3 shows the temperature gradient as a cosine function, with the energy equation describing heat transfer, is given by: a) Energy Equation where 𝑐𝑐𝑝𝑝 is the specific heat at constant pressure, 𝑇𝑇 is the temperature, 𝑘𝑘 is the thermal conductivity, and 𝑄𝑄 represents heat sources or sinks.

b) Mass Transfer Equation

Mass transfer, involving species concentration C, is described by the convection-diffusion equation:

where 𝐷𝐷 is the diffusion coefficient and 𝑅𝑅 represents reaction terms. Coupling these equations with the momentum equation allows for the analysis of how thermal and concentration gradients affect fluid flow .

Influence of Magnetic Fields on Fluid Flow:

Magnetic fields introduce additional forces in fluid flow, described by MHD equations. Figure 4 shows how different magnetic field strengths (B = 0, 0.5, 1.0) dampen the velocity profile, with stronger fields leading to greater damping. The Lorentz force F𝐿𝐿𝑜𝑜𝑟𝑟𝑒𝑒𝑛𝑛𝑡𝑡𝑧𝑧 is responsible for this effect is given by:

where 𝐽𝐽 is the current density and 𝐵𝐵 is the magnetic flux density. where 𝜎𝜎 is the electrical conductivity and 𝐸𝐸 is the electric field. Incorporating these into the momentum equation, the MHD momentum equation becomes:

Similarity Transformation Techniques:

Similarity transformation techniques simplify partial differential equations (PDEs) by combining independent variables (e.g., x and t) into a single similarity variable η.

Similarity transformations can reduce PDEs to ODEs. For example, the momentum equation might simplify to 𝑓𝑓′′(𝜂𝜂) + 12𝜂𝜂𝑓𝑓′(𝜂𝜂) = 0; where 𝑓𝑓(𝜂𝜂) is the transformed function.

Conversion of Partial Differential Equations to Ordinary Differential Equations

Converting PDEs to ODEs simplifies the solution process, making it more suitable for numerical analysis, often by assuming a specific solution form and applying boundary conditions. For example, if we assume a self-similar solution for the velocity field (𝑥𝑥, 𝑡𝑡), we might write:

where 𝛼𝛼 is a constant and 𝐹𝐹(𝜂𝜂) is a function of the similarity variable. Substituting this form into the original PDE, we obtain an ODE for (𝜂𝜂):

where 𝛽𝛽 and 𝛾𝛾 are constants derived from the original PDE coefficients. By solving this ODE, we can obtain insights into the behavior of the original fluid flow problem, greatly simplifying the analysis and computational effort required.

Methodology Model Development

Developing a comprehensive model for viscous fluid flow in porous media requires integrating factors like thermal conductivity, mass transfer, and magnetic fields. The model development begins by defining governing equations, including the continuity, momentum (Navier-Stokes), energy, and convection-diffusion equations.

Boundary conditions, such as inlet velocities, pressures, and temperature gradients, are crucial for defining fluid behavior at porous medium interfaces. The model also accounts for permeability variations, which significantly affect flow patterns.

Validation of Mathematical Models

Validating mathematical models involves comparing predictions with experimental data and numerical simulations. Controlled experiments measuring flow rates, pressure drops, and temperature distributions provide empirical data to verify the model's accuracy. Numerical simulations using CFD software validate the models by analyzing fluid flow under various conditions. Discrepancies between model predictions and experimental or numerical results are addressed through iterative model refinement.

Geometrical Configurations and Their Effects

The geometrical configuration of the porous medium affects fluid flow dynamics. Two-dimensional flows are simpler but less accurate, while three-dimensional flows offer better accuracy but require advanced numerical techniques and significant computational resources. Radial flow configurations, important in applications like groundwater flow and petroleum engineering, require consideration of radial symmetry and varying distances. Numerical simulations and experiments help study how geometry affects fluid dynamics in porous media.

Numerical Simulation Setup

Numerical simulations are set up by discretizing governing equations over a computational grid and applying boundary conditions. The finite volume method is used, dividing the domain into small control volumes and converting the equations into algebraic ones for numerical solution.

The simulation setup involves defining initial conditions, fluid properties, and porous medium characteristics. Advanced CFD software like ANSYS Fluent or OpenFOAM iteratively solves the discretized equations, with validation against experimental data to ensure accuracy and reliability.

Experimental Procedures

Experimental procedures involve selecting porous media with varying properties and using fluids of different viscosities to study flow behavior. The setup includes a fluid reservoir, pump, flow meters, pressure transducers, and temperature sensors to monitor flow rates, pressure drops, and temperature variations. Experiments are conducted under controlled conditions, varying parameters like flow rates and temperature gradients. Data collected at multiple points ensures accuracy, providing a foundation for validating and refining theoretical models.

Test Conditions: Table 1 outlines the experimental conditions, including porous medium type, fluid viscosity, flow rate, inlet temperature, and magnetic field strength.

Data Collection and Analysis

Data from experiments and simulations, including flow rates, pressure drops, temperature distributions, and velocity profiles, are analyzed to identify trends and validate mathematical models. Statistical analysis assesses data reliability, and discrepancies are addressed through model refinement. Validated models are used to predict fluid behavior for practical applications.

Results and Discussion

The comparison between numerical simulations and experimental data shows strong correlation, confirming the accuracy of the developed models for viscous fluid flow in porous media. Experimental velocity profiles, pressure drops, and temperature distributions align well with predictions, validating the governing equations and boundary conditions. This confirms the robustness of the mathematical framework. Table 2 summarizes experimental results, including pressure drop, outlet temperature, average velocity, and heat transfer rate. Table 3 provides detailed numerical data, covering initial and final pressures, temperatures, flow rates, viscosities, magnetic field strengths, velocities at the inlet and outlet, heat transfer rates, and porosity percentages. Figure 5(a) Pressure Drop vs. Flow Rate, plot shows that higher flow rates lead to greater pressure drops due to increased resistance in the porous medium. illustrates the effect of flow rate on temperature change, with higher flow rates leading to greater temperature variations. Figure 5(e) demonstrates how magnetic field strength impacts pressure drop, highlighting its influence on flow resistance and pressure dynamics. The validated models have significant implications for industries like chemical engineering, geothermal energy, and environmental remediation. They enable optimized system design and operation by accurately predicting fluid behavior in porous media. Incorporating thermal and mass transfer, along with magnetic field effects, offers a comprehensive approach to addressing fluid dynamics challenges and fostering innovative solutions across various industries.

Hypothesis Testing

The hypothesis that a comprehensive model incorporating heterogeneity, anisotropy, heat and mass transfer, and magnetic fields can accurately predict fluid behavior in porous media was validated. The results showed that including heterogeneity and anisotropy improves prediction accuracy compared to conventional models assuming homogeneity and isotropy.

The coupled heat and mass transfer analysis highlighted their critical role in fluid dynamics. Heat transfer, influenced by thermal conductivity, affects temperature distribution and flow patterns, while mass transfer, driven by concentration gradients, impacts overall flow behavior. Integrating these effects enhances model predictive capabilities, especially for thermal and chemical processes. Experimental validation confirmed the reliability of the theoretical models. Key parameters such as pressure drops, temperature variations, and flow rates were measured, and the empirical data closely aligned with the model predictions, ensuring the models' practical applicability. Numerical simulations using CFD confirmed the experimental results, providing detailed insights into the interactions between thermal effects, fluid viscosity, and porous medium properties. The alignment between numerical and experimental data validated the hypothesis that the comprehensive model accurately predicts fluid behavior.

Conclusion

This study presents a validated model for viscous fluid flow through porous media, incorporating factors like heterogeneity, anisotropy, heat and mass transfer, and magnetic fields. The model improves prediction accuracy for complex fluid behavior in real-world conditions. Experimental and numerical results confirm its reliability. Applications include chemical engineering, geothermal energy, and environmental remediation. Future research should refine the models and expand validation for broader applicability

References

1. Srivastava BG, Deo S. Effect of magnetic field on the viscous fluid flow in a channel filled with porous medium of variable permeability. Appl Math Comput. 2013;219:8959-64. doi:10.1016/j.amc.2013.03.065.
2. Beg OA, Zaman A, Ali N, Gaffar SA, Bég TA. Numerical computation of non-linear oscillatory two-immiscible magnetohydrodynamic flow in dual porous media system: FTCS and FEM study. Heat Transf Asian Res. 2019;48(7):293147.
3. Khedr MEM, Chamkha AJ, Bayomi M. MHD flow of a micropolar fluid past a stretched permeable surface with heat generation or absorption. Nonlinear Anal Model Control. 2009;14(1):27-40.
4. Chamkha AJ. MHD-free convection from a vertical plate embedded in a thermally stratified porous medium with Hall effects. Appl Math Model. 1997;21(10):603-9.
5. Magyari E, Chamkha AJ. Exact analytical results for the thermosolutal MHD Marangoni boundary layers. Int J Therm Sci. 2008;47(7):848-57.
6. Shahnazian H, Odenbach S. Rheological investigations of ferrofluids with a shear stress-controlled rheometer. J Phys Condens Matter. 2008;20:204137.
7. Hamdan MH, Kamel KT, Siyyam HI. A permeability function for Brinkman's equation. In: Proceedings of 11th WSEAS International Conference on Mathematics Methods and Computational Techniques in Intelligent Systems. 2009. p. 198-205.
8. Deen NG, Peters E, Padding JT, Kuipers J. Review of direct numerical simulation of fluid-particle mass, momentum, and heat transfer in dense gas-solid flows. Chem Eng Sci. 2014;116:710-24.
9. Verma VK, Datta S. Analytical solution of slow flow past a heterogeneous porous sphere with radial variation of permeability using Brinkman model. J Porous Media. 2012;15(4):355-9.
10. Wood BD. A primer on inertial flows in porous media. AGU Fall Meet Abstr. 2009;2009:H43E-1069. International Journal of Statistics and Applied Mathematics https://www.mathsjournal.com
11. Turkyilmazoglu M. MHD natural convection in saturated porous media with heat generation/absorption and thermal radiation: closed-form solutions. Arch Mech. 2019;71(1):49-64.
12. Raptis A, Perdikis C. Magnetohydrodynamics effects on mass transfer flow through porous medium. Astrophys Space Sci. 1985;113:53-8.
13. Turkyilmazoglu M. Velocity slip and entropy generation phenomena in thermal transport through metallic porous channel. J Non Equilib Thermodyn. 2020;45(4):341-55.
14. Yadav R, Agarwal M, Joseph J, Saxena P. Numerical analysis of MHD flow of viscous fluid between parallel porous bounding wall. Int J Adv Res. 2016;4(3):1448-57.