Mathematical Model Of Flow Research Articles

Statistical and mathematical modeling of viscoelastic 3D Casson nanofluid flow with activation energy and motile microbes induced by exponential sheet bounding with Darcy-Forchheimer porous medium

Statistical and mathematical modeling of viscoelastic 3D Casson nanofluid flow with activation energy and motile microbes induced by exponential sheet bounding with Darcy-Forchheimer porous medium

An Elastohydrodynamic model of the slot-die coating process

Abstract The slot-die coating process plays an important role in the industry, as it is employed in many different fields. The characteristics of the final application are determined by the flow between the die and the roller. This research paper aims to develop a mathematical model of such flow that takes into account the roller deformations caused by high pressure values reached by the coating fluid. This Elastohydrodynamic model is made up of a coupling between the mathematical model of the flow and the mathematical model of the roller deformations. Model resolution is undertaken numerically by deforming the flow domain according to the roller deformations using computational fluid dynamics (CFD) and computational solid mechanics (CSM) techniques. For its part, the finite volume method (FVM) is used to perform the flow model analysis while the finite element method (FEM) is employed to deal with roller deformations. The results obtained from this model give information on the flow pressure distribution, coating gaps, meniscus position, extent of roller deformations in the coating flow and the influence of different operating conditions. The information obtained from this study is valuable for industrial applications, as it gives insights into the coating process that can help manufacturers to define a suitable combination of operating parameters in order to obtain coating applications that meet quality and performance requirements.

Upper convected Maxwell fluid in an isothermal flow with thermal radiation

In this research paper, a mathematical model of the isothermal flow of upper convected Maxwell (UCM) fluid with thermal radiation through a saturated porous media was studied in detail, with different physical parameters being considered. Maple 18 software is used to implement this strategy. The analysis of the results showed that the flow system was significantly influenced by the Schmidt, Darcy, and Deborah numbers, the unsteadiness parameter, the Maxwell parameter, the magnetic field parameter, and other physical characteristics. The influence of several physical characteristics on the flow system are discussed and shown graphically.

Shape factor and thermal radiation effects on convective MHD flow of Casson blood ternary hybrid nano‐fluid through a stretched rotating rigid disk

AbstractThe mathematical modeling of blood flow with the incorporation of ternary hybrid nanoparticles (i.e., titania, silica, and alumina nanoparticles) is the focus of this work. The flow of blood is simulated using the Casson fluid model. The impacts of ternary hybrid nanoparticles, shape factor, and Geometry of solid nanoparticles are visualized and investigated. The momentum and thermal characteristics of a flowing liquid are determined by considering magnetization, porosity, and nonlinearized radiating thermal flux. The basic partial differential equations resulting from mathematical modeling are turned into non‐linear ordinary differential equations using suitable velocity transforms. The derived nonlinear equations are then numerically calculated using the 4th‐5th order Runge‐Kutta‐Fehlberg method with the shooting technique and analytically solved by the Adomian decomposition method (ADM). The effects of the factors involved on the dimensionless profiles produced are fully addressed. It is found that the Skin friction factor and heat transfer rate in the radial direction upsurge with the augment of both rotation parameter and nanoparticles volume fraction; however, the Skin friction factor drops in the azimuthal direction. Also, results obtained reveal an enhancement in the local Nusselt number with the upsurge in the magnitude of radiation parameter, Rd, solid nanoparticles concentration, , shape factor value, s, and disk temperature, . For validation, the outcomes of this inquiry were compared to the outcomes of the HAM‐based Mathematica software. In addition, the acquired analytical DRA data are compared to numerical RKF45 values and those given in the literature.

Pore-scale modeling of wettability alteration coupled two-phase flow in carbonate porous media

The Middle East marine carbonate reservoirs are the primary targets for global crude oil exploration and production. In contrast to the clastic rock reservoirs, the marine carbonate reservoirs have undergone more complicated deposition, diagenesis and transformation, and the rock surface tends to be oil-wet or mixed-wet. One of the most effective methods for enhancing oil recovery in marine carbonate reservoirs is the use of an ion matched (IM)-surfactant (S) system to alter the wettability of rock surfaces. In this paper, we propose a microscale mathematical model of wettability alteration-coupled fluid flow by considering oil–water two-phase flow, wettability alteration, solute advection- diffusion and surfactant adsorption. The proposed microscale mathematical model is solved using the multi-component lattice Boltzmann model, and its accuracy is further verified. Finally, the effects of displacing fluid, ion-matched water concentration and surfactant concentration on oil mobilization at pore-scale are studied. Results show that, ion-matched water has a positive effect on wettability alteration, and typically shifts the cross-point of oil–water relative permeability curve to the right. The concentration of ion-matched water primarily affects fluid morphology of the wetting-phase and has little effect on phase permeability. After adding surfactant to the ion-matched water, the Euler coefficient of oil-phase gradually increases as the displacement continues, a large amount of remaining oil is gradually dispersed into scattered oil droplets and the oil-phase relative permeability is further improved.

Analysis of thermophoretic and Brownian diffusions in hydromagnetic curvilinear flow of Carreau nanofluid with activation energy and heat generation

This study presents a theoretical analysis of the thermophoretic and Brownian diffusion effects on the magnetized flow of Carreau nanofluid over a stretchable curved surface. The investigation includes the impact of heat generation and activation energy, which are incorporated into the energy and concentration equations, respectively. The mathematical model of the boundary layer flow is formulated as a set of nonlinear partial differential equations using curvilinear coordinates. These equations are transformed into ordinary differential equations through the application of appropriate similarity variables. The resulting system is solved numerically by using the shooting method along with the Runge-Kutta technique. The accuracy of the numerical results is further verified by the Keller-box method, a finite difference technique. The obtained results are also checked by giving a comparison table with the published data in the literature. Furthermore, the convergence of the solutions is also checked by a convergent series method known as the homotopy analysis method (HAM). The study examines the influence of several key parameters, including the radius of curvature, Weissenberg number, Prandtl number, Lewis number, Brownian and thermophoretic diffusions parameters, Hartmann number, heat generation parameter, activation energy parameter, temperature difference parameter, chemical reaction parameter, and power law index. The effects of these parameters on momentum, temperature, concentration distributions, mass and heat transfer rates, and as well as the skin friction coefficient are examined through graphs and tables. Graphical observations declare that the velocity field diminishes against the rising values of the Hartmann number and Weissenberg number. However, the velocity profile rises with improved values of the radius of curvature parameter and power law index parameter. It is also noticed that the magnitude of the concentration field rises with growing values of the Brownian parameter, activation energy parameter and Lewis number, whereas it shows a decreasing manner against the enhanced values of the thermophoresis parameter and temperature difference parameter.

Investigations of the Cold Box Core Curing Stage Using an Augmented Simulation Approach

AbstractFoundry sand cores are used to produce complex metal castings. They are manufactured on universal coremaking machines in two distinct stages, shooting and curing. Various simulation tools exist to optimize the coremaking process. Numerical modeling of the two distinct process stages is challenging as it involves different tooling, raw materials, and process parameters. Additionally, each stage requires appropriate mathematical modeling of the fluid flow since the relevant physical phenomena are different. The state-of-the-art combined simulation approach (CA) can simulate these two stages in sequence using a Core shooting & curing module. The gassing system (GS) is also important to consider when calculating the complete physical mass flow during the curing stage in CA. An augmented simulation approach (AA) is proposed, which can include GS in the coremaking simulations. This approach enables the execution of two-stage combined simulations (CA) using Magma C + M software (SW1) and incorporates GS data using FLOW-3D CAST software (SW2). In order to analyze the influence of including selected factors such as the hopper during the shooting stage and GS during the curing stage, the simulation approaches CA and AA were compared. The simulation results were correlated with experiment results (EX) to analyze sand density distribution and different curing times, focusing on the final core quality. Quantity analysis has been done among EX samples to observe the cured core trend. The results obtained from AA exhibit a significantly better correlation with EX than with CA. The proposed augmented approach offers significant potential for the improved analysis, optimization, and accurate prediction of complex coremaking processes within polyurethane (PU) Cold Box Systems.

Mathematical model of molten iron and slag flow in the hearth

The stable hearth discharge is crucial for the operation of a blast furnace (BF). This study established a mathematical model for hearth discharge based on Bernoulli's principle, which considers the blast pressure, kinetic energy of molten iron and slag, gravitational potential energy, the pressure loss of molten iron and slag flowing through the deadman and taphole. It is found that the sharp increase of molten iron and slag flow is probably due to the fracture of the mud, while the sharp decrease of molten iron and slag flow is probably due to the blockage of the taphole by coke. This paper also studies the effect of operating parameters on the hearth discharge. With the increase of taphole plug time, the production rate and slag viscosity, the tap-end time will increase. With the increase of deadman voidage, the tap-end time will be shortened. The correctness of the mathematical model was verified by the actual flow rate of slag and molten iron.

Simulation and analysis of proppant transportation and distribution rules under multi-cluster perforation in horizontal wellbore

Unconventional reservoirs require segmented multi-cluster fracturing technology for commercial development. Due to the settling and inertial effects of proppant particle flow, even if the fracturing fluid is evenly distributed among each perforation, the proppant may still be unevenly distributed. However, the uniform distribution of proppant in multi-fracture is a favorable condition for high-yield gas wells. Based on the proppant dynamics model in the wellbore, coupled with the non-uniform fracture propagation model, a mathematical model of proppant wellbore transportation flow under multi-fracture is established. Through numerical calculations, the proppant transport and distribution mechanism are studied. The study shows that due to gravity, distribution in the upstream is more uniform, and the proppant in the downstream gradually settles and most particles flow toward the bottom. Proppant distributed in each cluster is affected by the flow rate of slurry to a certain extent, and is also controlled by its own properties. The turning efficiency determines the inflow of proppant in perforations. Increasing the pump rate can alleviate gravity but reduce turning efficiency. The perforation azimuth primarily regulates the flow of proppant into downstream perforations. The research results can provide theoretical guidance for the design of proppant transport and cluster perforation in hydraulic fracturing.

Numerical simulation of two-phase oil–water flow in fractured-vuggy reservoirs based on the coefficient of porous medium proportion and coupled regions

Based on the flow characteristics of fluids in various reservoir media, fractured-vuggy oil reservoirs can be classified into seepage zones and conduit flow zones. An interface exists between these two regions, where the movement of formation fluid near this interface is characterized by a coupling or transitional phenomenon between seepage and conduit flow. However, the complexity of this coupling interface poses challenges for traditional numerical simulations in accurately representing the intricate fluid dynamics within fractured-vuggy oil reservoirs. This limitation impacts the development planning and production adjustment strategies for fractured-vuggy oil reservoirs. Consequently, achieving accurate characterization and numerical simulation of these systems remains a critical challenge that requires urgent attention. A new mathematical model for oil-water two-phase flow in fractured-vuggy oil reservoirs is presented, which developed based on a novel coupling method. The model introduces the concept of the proportion coefficient of porous media within unit grids and defines a coupling region. It employs an enhanced Stokes–Brinkman equation to address the coupling issue by incorporating the proportion coefficient of porous media, thereby facilitating a more accurate description of the coupling interface through the use of the coupling region. Additionally, this proportion coefficient characterizes the unfilled cave boundary, simplifying the representation of model boundary conditions. The secondary development on the open-source fluid dynamics software is conducted by using matrix & laboratory (MATLAB). The governing equations of the mathematical model are discretized utilizing finite volume methods and applying staggered grid techniques along with a semi-implicit calculation format for pressure coupling—the Semi-Implicit Method for Pressure Linked Equations algorithm—to solve for both pressure and velocity fields. Under identical mechanism models, comparisons between simulation results from this two-phase flow program and those obtained from Eclipse reveal that our program demonstrates superior performance in accurately depicting flow states within unfilled caves, thus validating its numerical simulation outcomes for two-phase flow in fractured cave reservoirs. Utilizing the S48 fault-dipole unit as a case study, this research conducted numerical simulations to investigate the water-in-place (WIP) behavior in fractured-vuggy oil reservoirs. The primary focus was on analyzing the upward trend of WIP and its influencing factors during production across various combinations of fractures and dipoles, thereby validating the feasibility of the numerical modeling approach in real-world reservoirs. The simulation results indicated that when multiple dissolution cavities at different locations communicated with the well bottom sequentially, the WIP in the production well exhibited a staircase-like increase. Furthermore, as the distance between bottom water and well bottom increased, its effect on water intrusion into the well diminished, leading to a slower variation in the WIP curve. These characteristics manifested as sudden influxes of water flooding, rapid increases in water levels, and gradual rises—all consistent with actual field production observations. The newly established numerical simulation method for fractured-vuggy oil reservoirs quantitatively describes two-phase flow dynamics within these systems, thus effectively predicting their production behaviors and providing guidance aimed at enhancing recovery rates typically observed in fractured-vuggy oil reservoirs.

Thermal recovery of heavy oil reservoirs: Modeling of flow and heat transfer characteristics of superheated steam in full-length concentric dual-tubing wells

Accurately predicting the flow and heat transfer characteristics of superheated steam (SHS) in the wellbore is essential for efficiently developing heavy oil reservoirs. However, few studies have examined SHS flow in full-length concentric dual-tubing wells (CDTWs). In this paper, firstly, based on fluid dynamics and thermodynamics theories, a mathematical model for SHS flow in vertical and horizontal wellbores is proposed. Then, this model is solved using a finite difference method for spatial discretization and an iterative technique. Finally, model verification, type curve analysis, and sensitivity analysis are conducted sequentially. The results indicate that SHS temperature and pressure are independent variables, and the effect of friction energy on temperature is greater than its impact on pressure. The injection rate has a critical value, and the critical injection rates of the main controlling factors affecting the pressure drop and temperature drop of SHS in the vertical tubing are different. When the injection rate is below the critical value, the gravitational force of the SHS helps maintain high enthalpy. The SHS should be transported to the well bottom as quickly as possible. To enhance the uniform heating effect of the reservoir and improve heat utilization efficiency, a relatively small injection rate, high injection pressure, and low injection temperature are recommended. The uniformity of SHS pressure and temperature distribution in the horizontal annulus positively correlates with the uniformity of the reservoir's heat absorption rate. Achieving more uniform SHS pressure and temperature profiles in the horizontal annulus benefits the uniform heating of the reservoir.

On the mechanics of conducting micropolar fluids with magnetic particles: Vorticity–microrotation difference

In this paper, mechanisms that differentiate microrotation (i.e., particles' internal rotation) and vorticity in magnetohydrodynamic micropolar flows with magnetic particles are examined. Micropolar fluids are characterized by the asynchronization of the microrotation and the vorticity of the liquid carrier. When the microrotation is equal to the fluid's vorticity, the flow is purely Newtonian. In this context, two classical examples of micropolar magnetohydrodynamic fluids with magnetic particles are used, namely, blood and ferrofluid. The effect of all dimensionless parameters associated with the mathematical model is examined in terms of differentiating microrotation and vorticity. These dimensionless parameters are specifically calculated using the physical properties of the ferrofluid and blood. It was shown that higher values of the rotational viscosity (which are associated with smaller values of the spin relaxation time), higher values of channel's height, and higher values of the microrotation wall parameter tend to equalize microrotation and vorticity. On the other hand, the spin viscosity and the micromagnetorotation (magnetic torque) are mechanisms that differentiate microrotation and vorticity. Lorentz force does not seem to have any noticeable effect on the microrotation–vorticity difference. It is anticipated that this study will reveal the cases where a magnetohydrodynamic micropolar fluid with magnetic particles, such as blood and ferrofluid, can be simplified to a Newtonian one, which brings many benefits associated with a simplified mathematical flow model (such as smaller computational cost and less time).

Mathematical modeling and analysis for radiative MHD peristaltic flow of Bingham nanofluid

Mathematical modeling and analysis for radiative MHD peristaltic flow of Bingham nanofluid

Mathematical Model of Drilling Mud Movement along a Shale Shaker Screen

Introduction. When drilling oil or gas wells, rock from the bottomhole is brought to the surface by drilling mud, which is cleaned of sludge by shale shakers at the first stage. Shale shakers are equipped with such a screen and create such trajectories of vibration of the frame in order to solve a dual problem: on the one hand, to provide the cleaning of the drilling mud coming from the well mouth, and on the other hand — to maintain the proper degree of cleaning. To correctly select the nomenclature of screens, it is necessary to reliably determine the throughput capacity of a shale shaker in the real-time setting. This will allow you to pre-order the required number of screens with the required cell size. Previously, studies were conducted by sieving a mixture of calibrated granules on a standard shaker or by straining mineral oil through a fixed screen. However, this does not fully correspond to the actual conditions of drilling mud screening. The objective of this article was to simulate the throughput of a drilling shale shaker under real conditions by calculating the movement of drilling mud along a vibrating screen using its previously obtained specific throughput capacity.Materials and Methods. When creating a mathematical model of the full throughput capacity of a shale shaker, a previously known mathematical model of fluid flow in an open channel and a finite-dimensional approximation in the form of a sequence of sections of concentrated containers connected by hydraulic conductivities (by analogy with the approximation of long electrical lines and extended gas pipelines by the finite element method) were used. The fluid flow rate over the screen was determined according to Chézy law. In this case, the cross-sectional area of the flow above the screen for a specific section was specified as the arithmetic mean between the initial and final values of the drilling mud height in the corresponding section.Results. A chain calculation scheme for sieving drilling mud was compiled. Based on the material balance, mathematical models were developed for determining the throughput and height of drilling mud above the screen: for the first section of the screen and the first nodal point of the screen; for subsequent screen sections and subsequent screen nodal points; for the last section of the screen and the last nodal point of the screen.A mathematical model of the drilling mud movement along the shale shaker screen has been obtained, which, due to the finite-dimensional representation of this flow along n sections of the screen length, is a system of n integral equations describing concentrated tanks, and n – 1 algebraic equations describing the flow of drilling mud between tanks.Discussion and Conclusion. The movement of drilling mud along the length of the shale shaker screen corresponds to the movement of a steady uneven flow in an open channel. However, in the first case, the volume flow along the length is variable (decreasing to zero), and in the second case, it is constant along the length. Therefore, the use of only one differential equation for steady uneven flow in an open channel is not sufficient in this case. The resulting mathematical model of the drilling mud movement along the shale shaker screen is a mathematical basis for the development of a modeling program to determine the throughput of the shale shaker for real operating conditions, i.e., for a specific drilling mud, a given screen or a given drive system of the shale shaker.

Impact of thermal conductivity of a Williamson’s nanofluid flow due to a stretching sheet

This paper intends to build a mathematical model of a two-dimensional Williamson nanofluid flow due to the stretching of a sheet linearly. The flow geometry is influenced by the magnetic field that is applied externally to the system. The induced magnetic field is infinitesimal and hence neglected. However, the flow structure is incorporated with the Brownian motion and thermophoresis effects as it alter the physics of the flow. The equations that model the flow comprise the highly nonlinear coupled simultaneous system. Thus, a numerical technique namely finite difference accompanied with the Thomas algorithm is adopted to approach the flow system. The motion, temperature, and concentration of the Williamson nanofluid flow are studied for the different flow controlling parameters. The co-efficient of skin-friction, heat, and mass transfer rates is also computed. The streamlines, isotherms, and iso-concentration are plotted to picturise the flow phenomena in the complete domain. The study reveals that the temperature of the flow raises with the Brownian motion of the nanoparticles and the trend is opposite with the concentration. The strength of the streamlines, isotherms, and the concentration contour is identified to be high for the least magnitude of Weissenberg number. The Brownian motion raises the heat transfer rate and slows down the mass transfer rate.

Numerical study of Casson dusty hybrid nanofluid flow with volume fraction over a Darcy–Forchheimer porous medium

AbstractThis study develops a mathematical model for dusty Casson hybrid nanofluid flow over a vertical stretching sheet, considering the influence of volume fractions for the dust particles and nanoparticles, magnetic field and Darcy–Forchheimer penetrating medium. The nanoparticles used are TC4 (Ti‐6Al‐4V) and Nichrome (80% Ni and 20% Cr), suspended in unused engine oil (EO). The governing partial differential equations are transformed into ordinary ones using suitable transforms. Solutions are obtained numerically via the Matlab built‐in program Bvp4c technique and represent graphically the velocity, temperature, and concentration profiles for both the fluid and dust phases. Additionally, tabular forms present wall shear stress and heat transfer rate variations under different parameter values. It is found that the enhancing volume fraction parameter helps to boost the magnitude of shear stress, heat, and mass diffusivity.

Flow and heat transfer insights into a chemically reactive micropolar Williamson ternary hybrid nanofluid with cross-diffusion theory

Abstract The need for efficient nanotechnology has led to unexpected developments. Conserving continuous thermal propagation is essential in many industrial and thermal systems because it improves the efficiency of thermal engineering engines and machinery. Therefore, a promising platform to increase thermal power energy is the hybridization of magnetic nanoparticles in a heat-supporting, non-Newtonian fluid. In light of the above applications, a mathematical model is established to analyze the variable fluid features of the thermally radiative and chemically reactive flow of a micropolar Williamson ternary hybrid nanofluid with electromagnetohydrodynamic and electroosmosis forces on a porous stretching surface. Stratification boundary conditions and variable fluid properties were used to analyze the thermal and solutal behavior of the fluid flow. Furthermore, to measure the disorder of the flow system, entropy generation was considered by the impact of Joule heating and viscous dissipation. To develop the numerical scheme BVP4C in MATLAB, we first converted the mathematical flow model into two ordinary differential equations using a suitable transformation. The graphical and numerical results were determined against several parameters of a ternary hybrid nanofluid ( MWCNT , A l 2 O 3 , SiC ) ({\rm{MWCNT}},\hspace{0.25em}{\rm{A}}{{\rm{l}}}_{2}{{\rm{O}}}_{3},\hspace{0.25em}{\rm{SiC}}) and unary nanofluid ( A l 2 O 3 ) ({\rm{A}}{{\rm{l}}}_{2}{{\rm{O}}}_{3}) . The results indicate that the heat transfer rate is more prominent in the ternary hybrid nanofluid than in the unary nanofluid because the addition of nanofluids to the base fluid is used to improve the heat transport rate. It can be seen from the figures that a greater estimation of the magnetic and electric field parameters improves the fluid velocity because, for low values of M ≤ 1 M\le 1 , the aiding force is dominant compared to the retarding force, which results in an improvement in the velocity profile. Furthermore, the entropy generation rate increases for higher values of the Brinkman number and temperature ratio parameter because more heat is produced due to the greater values of Br {\rm{Br}} .

Stability analysis of dual solutions for nonlinear radiative magnetohydrodynamic flow of [formula omitted] hybrid nanofluid over a nonlinearly shrinking surface

The present work explores the existence and temporal stability of dual solutions with the consequence of nonlinear thermal radiation on hybrid nanofluid flow (including Silver and Titanium dioxide nanoparticles in water) past a permeable shrinking sheet. The analysis also considers the united outcomes of suction/injection, Joule heating, viscous dissipation, and magnetohydrodynamics. Employing realistic assumptions and suitable similarity transformations, the governing nonlinear partial differential equations are formulated and altered into a nonlinear ordinary differential equations system. These equations are then solved numerically using the Lobatto IIIA technique. It is observed that dual solutions can occur within a precise range of the shrinking surface parameter. Analyzing the temporal stability of these dual solutions under small disturbances shows that the upper solution branch is stable and represents a physically realistic solution to the problem. In contrast, the lower solution branch is unstable and not physically realizable. No solution exists beyond the critical λc=−1.3915 for the shrinking parameter. The critical values of shrinking parameter λc=−1.2565,−1.3188, and −1.4161 correspond to volume fractions of hybrid nanofluid with φAg−TiO2 of 2%,4%, and 8%, respectively. No solution exists beyond this critical point. The graphical representation and quantitative explanation demonstrate how various emergent characteristics affect momentum and thermal boundary layer profiles, Nusselt number, and skin friction. The velocity and temperature profile of hybrid nanofluid can be improved by adding a small volume of nanoparticle volume fraction. Thermal boundary layer thickness is amplified by an accumulation in the shrinking parameter, power index of velocity component, thermal radiation parameter, and temperature difference ratio. Sheering stress and heat transfer coefficient improve with adding more nanoparticles in the base fluid for the first solution braches, but opposite effects are found in the second solutions. The contributions of this research are significant both theoretically, in terms of advancing the mathematical modeling of hybrid nanofluid flow with heat transfer in engineering systems, and practically, in terms of applications to engineering cooling systems.

Mathematical modelling and computational reduction of molten glass fluid flow in a furnace melting basin

Mathematical modelling and computational reduction of molten glass fluid flow in a furnace melting basin

A Mathematical Model of Interstitial Fluid Flow and Retinal Tissue Deformation in Macular Edema.

This study aims to develop a mathematical model to elucidate fluid circulation in the retina, focusing on the movement of interstitial fluid (comprising water and albumin) to understand the mechanisms underlying exudative macular edema (EME). The model integrates physiological factors such as retinal pigment epithelium (RPE) pumping, osmotic pressure gradients, and tissue deformation. It accounts for spatial variability in hydraulic conductivity (HC) across the retina and incorporates the structural role of Müller cells (MCs) in maintaining retinal stability. The model predicts that tissue deformation is maximal at the center of the fovea despite fluid exudation from blood capillaries occurring elsewhere, aligning with clinical observations. Additionally, the model suggests that spatial variability in HC across the thickness of the retina plays a protective role against fluid accumulation in the fovea. Despite inherent simplifications and uncertainties in parameter values, this study represents a step toward understanding the pathophysiology of EME. The findings provide insights into the mechanisms underlying fluid dynamics in the retina and fluid accumulation in the foveal region, showing that the specific conformation of Müller cells is likely to play a key role.