Horizontal Barrier Research Articles

Numerical Study on the Effects of Horizontal Barriers on Fire Behavior in Logistics Facility Racks

In this study, a numerical analysis was conducted to derive the optimal installation conditions for flame-spread barriers within racks, to reduce the damage caused by recurrent fire accidents in logistics facilities. The effectiveness of the fire extinguishing system according to domestic regulations was reviewed to ensure the feasibility of installation of the horizontal barrier. This review revealed that sprinkler systems were ineffective in delaying or preventing the initial spread of fire under certain fire conditions, confirming the need for horizontal barriers. Furthermore, although the installation of horizontal barriers was effective in delaying the initial fire spread, it led to an increase in fire size by blocking the extinguishing water and altering the fire spread path. It was confirmed that deriving the optimal installation conditions for barriers requires the consideration of the interactions between the in-rack sprinkler heads and the surrounding fire extinguishing systems.

Performance Improvement of an Integrated Collector Storage Solar Water Heater by Using an Additional Glass Cover: A Numerical Study

This paper describes the analysis of a solar storage water heating system, emphasizing its cost-effectiveness and integrated design. The study employs numerical analysis to assess and optimize the performance of a rectangular solar collector with an added glass cover and an internal horizontal barrier. Using COMSOL software, a three-dimensional modeling study was conducted in the Kufa-Najaf climate on May 9th. The analysis focused on fundamental convection equations governing mass, momentum, and energy conservation. Temperature and velocity distributions were obtained at 13:30 PM, revealing a peak storage water temperature of 56.3 °C. The additional glass cover-barrier storage collector, set at a 60° angle, exhibited the highest stored energy at 2065.6 W, surpassing the conventional storage collector at 1007.06 W. Furthermore, the model demonstrated an increased instantaneous efficiency of 58.3%, compared to 28.4% for the conventional collector. The simulation results align well with previous findings, showing a maximum mean absolute percent error (MAE) of 4.5%.

Integrated operation of sand dam and groundwater dam to increase water supply capacity in mountainous areas of Chuncheon, South Korea

Study regionA valley area in South Korea Study focusSand dams are cost-effective water-storage facilities for water-stressed areas in arid and semi-arid climates, with the advantages of reducing evaporation losses and filtering pollutants. Generally, sand dams are constructed by placing a horizontal barrier across the river bed, creating a naturally sand-filled structure upstream. Unlike this typical construction, a sand dam in Korea in the monsoon climate region was built near a mountainous valley, not in a stream, in Mulori, Chuncheon City, to ensure a sustainable water supply throughout the year, and especially to prevent freezing in the winter. Furthermore, a small-scale groundwater dam was additionally constructed downstream of the sand dam to increase the water supply capacity in preparation for severe droughts. In this study, the hydrological reservoir-routing simulations were performed to evaluate the effects of the combined use of the two dams on increasing water supply. New hydrological insights for the regionDaily discharge rates for the period from March 16, 2020 to July 6, 2022 were simulated for several hybrid operations that supplies water primarily through the sand dam under the normal conditions, but combines both water sources to manage severe droughts. The simulated results showed that, compared with operating the sand dam alone, the hybrid operations could satisfy the village’s water demand by increasing the minimum water supply rates by 24.6–50.0 m³/day for different simulation conditions. The combined use of sand and groundwater dams presented in this study will help provide sustainable water supply in water-stressed areas, especially mountainous uplands vulnerable to droughts and freezes.

Analytical model for the mitigation of VOC vapor with horizontal permeable reactive barrier in the contaminated site considering non-uniform source

Volatile organic compounds (VOCs) contamination at the groundwater may cause vapor intrusion and pose significant threats to human health. As a novel low-carbon mitigation technology, a horizontal permeable reactive barrier (HPRB) is proposed to remove the VOC vapor in the vadose zone and mitigate the vapor intrusion risk. To estimate the performance of HPRB in the contaminated site with a non-uniform source, a transient two-dimensional analytical model is developed in this study to simulate the VOC vapor migration and oxidation processes in the layered soil. The analytical model is verified against the experimental results and numerical simulation first and the parameter study is then conducted. The HPRB has good performance for the contaminated sites involving factors including deep source and local soil with low effective diffusivity. To consider the vertical heterogeneity of the local soil, the traditional equivalent homogeneity method has limitations in considering the horizontal migration of VOC vapor and is not suitable for the two-dimensional model. On the contrary, the artificial layered method based on the proposed analytical model has better accuracy and is recommended to be adopted in practice. Leading to the exponential decrease in the VOC vapor concentration at the ground surface, increasing the thickness of HPRB is an effective measure to enhance the performance of HPRB. The fitting exponential function can be applied to determine the minimum design value of the thickness of HPRB in practice.

Ventilation performance evaluation of an operating room with temperature-controlled airflow system in contaminant control: A numerical study

This article investigates the efficacy of temperature-controlled airflow systems in modern operating rooms for contaminant control, a critical factor in preventing surgical site infections. We have conducted experimental measurements in an operating room equipped with temperature-controlled ventilation to map the airflow field and contaminant dispersion (airborne particles with diameters ranging from 0.5 to 1 μm). The results were used to validate the computational fluid dynamics code, which was then employed to simulate and examine different conditions, including contaminant release locations and air supply rates. Realizable k-epsilon and passive scalar models were utilized to simulate airflow and airborne particle phases. We assessed the airflow distribution and contaminant dispersion, utilizing indices such as ventilation and air change efficiency scales. The analysis provided quantitative insights into the distribution and removal of contaminants, as well as the speed at which the room air was replaced. Contamination was found to be effectively reduced when contaminants were released near exhaust outlets or under central unidirectional inlets. The presence of the operating table caused a big distortion of the central downward airflow, forming a horizontal air barrier at the periphery. Under this unique interior configuration, an appropriate air supply ratio between central and periphery zones was required to achieve optimal overall ventilation performance.

Numerical Simulation of Fracture Propagation during Temporary Plugging Staged Fracturing in Tight-Oil Horizontal Wells.

Fracture propagation with temporary plugging hydraulic fracturing in tight-oil reservoirs is simulated in this study. The research considers dynamic fluid redistribution, with stress differences among multiple fractures. The fracture morphology during temporary plugging staged fracturing (TPSF) is investigated by using a user-defined perforation element combined with a pore-pressure finite-element model. The precision of the integrated model is verified by using the standard finite-element approach. Then, case studies are presented to investigate the influence of cluster spacing, horizontal stress difference coefficients (SDC), injection rates, and barrier tensile strengths. The simulation results show that central cluster fractures are hampered by side-cluster fractures, while TPSF can alleviate the effect and lead to a more uniform propagation of all fractures. Stress interference weakens as cluster spacing increases, and propagation patterns are minimally influenced once spacing reaches 40 m. Higher injection rates can improve the injection pressure, enlarging fracture width, and potentially increasing the risk of fracture penetration. Barrier tensile strength and horizontal SDC can modify fracture geometries and determine the penetration behavior of multiple fractures.

Dispersion effects in porous medium gravity currents experiencing local drainage

We develop a theoretical model to study (dense) two-dimensional gravity current flow in a laterally extensive porous medium experiencing leakage through a discrete fissure situated along this boundary at some finite distance from the injection point. Our model, which derives from the depth-averaged mass and buoyancy equations in conjunction with Darcy's law, considers dispersive mixing between the gravity current and the surrounding ambient by allowing two different gravity current phases. Thus do we define a bulk phase consisting of fluid whose density is close to that of the source fluid and a dispersed phase consisting of fluid whose density is close to that of the ambient. We characterize the degree of dispersion by estimating, as a function of time, the buoyancy of the dispersed phase and the separation distance between the bulk nose and the dispersed nose. On this basis, it can be shown that the amount of dispersion depends on the flow conditions upstream of the fissure, the fissure permeability and the vertical and horizontal extents of the fissure. We also show that dispersion is larger when the gravity current propagates along an inclined barrier rather than along a horizontal barrier. Model predictions are fitted against numerical simulations. The simulations in question are performed using COMSOL and consider different inclination angles and fissure and upstream flow conditions. Our study is motivated by processes related to underground $\mathrm {H}_2$ storage e.g. an irrecoverable loss of $\mathrm {H}_2$ when it is injected into the cushion gas saturating an otherwise depleted natural gas reservoir.

METHOD FOR PREDICTION OF DAMAGE TO BUILDING STRUCTURES OF BUILDINGS AND STRUCTURES UNDER VIBRODYNAMIC IMPACTS TRANSMITTED THROUGH THE GROUND ENVIRONMENT, AND WAYS TO ELIMINATE THEM

A method for determining the possibility of damage to the building structures of buildings and structures in the case of transmission of vibrodynamic effects on them through the foundation from external vibration sources is proposed. The essence of the technique is to determine the maximum vibration speed of the foundation and the soil in front of it by the finite element method, followed by a comparison of the obtained value with the limited vibration speed. Based on this comparison, 4 possible options are considered: there are no damages, minor, significant or critical defects appear. In the last 3 cases, measures are taken to protect building structures from vibration. 3 variants of vibration isolation are considered, which allow to reduce the rate of vertical oscillations by more than 80%: installation of a horizontal inertial barrier in the form of a deepened concrete slab; installation of a vertical barrier in the form of adjoining boreholes filled with easily compressible materials or cased with gas-filled cylinders under pressure; a combination of these methods. It is shown that the use of these methods of vibration isolation makes it possible to avoid damage to buildings in case of exceeding the limiting speed of vibrations.

Evaluation of Effectiveness of Application of Horizontal Inertial Barriers to Reduce Vibrations Propagating in Ground Environment

The paper presents calculations of vibrations of the soil mass when a horizontal inertial barrier in the form of a rectangular concrete slab buried in the ground is placed on the propagation path of vibrations. The damping effect of a surface wave upon its contact with an inertial plate is associated with its reflection, refraction and partial absorption. Theoretical studies have been carried out using the finite element method. The ground medium has been considered as an elastic inertial array bounded by non-reflecting boundaries. Various variants of the geometry of the inertial plate and its spatial arrangement on the surface of the soil medium between the vibration source and the considered point behind the barrier are modeled. The effectiveness of each variant of vibration isolation has been quantified by the value showing how many times the speed of vertical oscillations of the ground behind the barrier decreases compared to the free propagation of surface waves. It is shown that an intensive decrease in vertical displacements occurs starting from the side face of the inertial plate. Here, the amplitude of vertical oscillations decreases by 9,8 times for a concrete slab 15 m long, and by 4,2 times for a 3-meter slab. At a distance of 22 m from the point of application of the dynamic load, the amplitudes decrease by factors of 5.48 and 2.95, respectively, for slabs of 15 and 3 m width. This method of reducing vibro-dynamic effects has a simple design and can be used in cramped conditions of urban development to protect existing and planned buildings and structures.

Immune cell extravasation in an organ-on-chip to model lung inflammation

Acute respiratory distress syndrome (ARDS) is a severe lung condition with high mortality and various causes, including lung infection. No specific treatment is currently available and more research aimed at better understanding the pathophysiology of ARDS is needed. Most lung-on-chip models that aim at mimicking the air-blood barrier are designed with a horizontal barrier through which immune cells can migrate vertically, making it challenging to visualize and investigate their migration. In addition, these models often lack a barrier of natural protein-derived extracellular matrix (ECM) suitable for live cell imaging to investigate ECM-dependent migration of immune cells as seen in ARDS. This study reports a novel inflammation-on-chip model with live cell imaging of immune cell extravasation and migration during lung inflammation. The three-channel perfusable inflammation-on-chip system mimics the lung endothelial barrier, the ECM environment and the (inflamed) lung epithelial barrier. A chemotactic gradient was established across the ECM hydrogel, leading to the migration of immune cells through the endothelial barrier. We found that immune cell extravasation depends on the presence of an endothelial barrier, on the ECM density and stiffness, and on the flow profile. In particular, bidirectional flow, broadly used in association with rocking platforms, was found to significantly delay extravasation of immune cells in contrast to unidirectional flow. Extravasation was increased in the presence of lung epithelial tissue. This model is currently used to study inflammation-induced immune cell migration but can be used to study infection-induced immune cell migration under different conditions, such as ECM composition, density and stiffness, type of infectious agents used, and the presence of organ-specific cell types.

Insights into catalyst-free, highly effective degradation of pharmaceutical contaminant in aqueous solution by a dielectric barrier discharge system

The emerging pharmaceutical wastewater including non-steroidal anti-inflammatory drugs has posed a serious threat to the ecological environment and human health, and therefore effective treatment technologies are being urgently needed. Herein, a novel horizontal continuous-flow dielectric barrier discharge (DBD) system without catalyst was developed for the highly effective degradation of diclofenac in aqueous solution. The degradation efficiency of 10 mg/L diclofenac reached about 90% in 14 min, and the corresponding energy efficiency was 0.181 g/kWh. Furthermore, the treated water quality was evaluated in terms of TOC, BOD5/COD and toxicity. The biodegradability rapidly improved with the increasing of DBD treatment time as be expected, although the mineralization efficiency of diclofenac was still low. To gain insight into the degradation mechanism, short-lived transient and long-lived species have been identified during DBD process, in which transient species (·OH, ·O2–) were the major contribution to the oxidation of diclofenac. Furthermore, based on the intermediates identified by LC-MS, a possible degradation pathway of diclofenac by DBD system in aqueous solution was proposed. Therefore, the novel DBD reactor with horizontal continuous-flow appears to be more promising for the treatment of pharmaceutical wastewater in practical applications.

Laboratory Experiment on Formation of MICP Horizontal Seepage-Reducing Body in Confined Aquifer for Deep Excavation

Using a soil improvement biotechnology based on microbial-induced carbonate precipitation (MICP) to form a horizontal water barrier in foundation pit dewatering can not only achieve a good dewatering control effect, but also control the cost and obtain good economic and environmental benefits. In this study, taking the Shanghai layer ⑨ sand as the research object, the solution injected into the layer ⑨ sand by MICP grouting technology includes bacterial solution (Sporosarcina pasteurii) and cementation solution (urea and CaCl2 solution). The migration mechanism of the bacterial solution and cementation solution in groundwater was analyzed through batch experiment and column test, and the mathematical model of migration was established. The dynamic adsorption characteristics of the S. pasteurii, Ca2+, and urea were obtained by dynamic adsorption test. Through isothermal adsorption test, the Freundlich isotherm was used to describe the isothermal adsorption of S. pasteurii onto fine sand, and the partition coefficient of the S. pasteurii was found to be 1.5 × 10−7 L/mg. The Langmuir isotherm can be used truly reflect the isothermal adsorption of Ca2+ onto fine sand. The distribution coefficient of the Ca2+ was 7 × 10−8 L/mg. The maximum adsorption capacity of the solute was 2404.8 mg/kg. The adsorption capacity of the urea was very small. According to the sand column test, the vertical dispersion was 11.5 cm. Due to the obvious size effect of dispersion, the sensitivity of dispersion in actual working conditions should be analyzed in a subsequent numerical simulation. The determined solute transport parameters provide support for the further study of the numerical simulation of the solute transport process of the HSRB formed by MICP grouting technology and provide a basis for further field application.

Breakdown path control by barrier effect for improving withstand characteristics of air‐insulated gaps under alternating voltages

A polymer barrier inserted in the interelectrode gap can significantly improve the insulation performance of air-insulated high-voltage electrical equipment. With the development of power equipment towards integration and miniaturisation, multiple barriers can be used to further increase the insulation level. This paper investigated the withstand characteristics of the air-insulated rod-rod gap with multiple barriers by changing the barrier numbers, the barrier position and the barrier insertion direction. Four control modes were set up based on the assumed breakdown paths under alternating voltages. It was found that the control modes with barriers inserted from the opposite direction had a more significant improvement effect than that from the same direction. When the horizontal barrier position was closer to the electrode and the vertical barrier position was farther to the electrode axis, the breakdown voltage was higher. The vertical barrier position was a key factor affecting the probability of breakdown path type, and the breakdown voltages along different types of breakdown paths were almost the same when the barriers were at the same position. The increase of the breakdown voltage was mainly resulted from the change of electric field distribution in the insulated gap and the extension of the breakdown path length.

Submerged breakwater of a flexible porous membrane with a vertical flexible porous wall over variable bottom topography

An analytical model of oblique wave interaction with a breakwater consisting of a submerged horizontal flexible porous membrane near a vertical porous barrier over variable bottom topography is developed based on linear water wave theory. Using Green's second identity, a mild-slope equation for oblique waves propagating over a submerged porous membrane and a finite variable bottom topography is derived. The analytical solution of the referred problem is obtained by using the eigenfunction expansion method along with the mild-slope equation and Darcy's law is used for wave past horizontal and vertical porous structure. The present results are compared with existing published analytical results and experimental data and observed that they have a good level of agreement. The effect of different design parameters on the wave quantities and wave forces on the porous wall in variable bottom depth is analyzed. It is found that full-wave reflection is drastically reduced and the wave energy is dissipated by the submerged horizontal flexible porous membrane and porous barrier. The analysis of the reflection and dissipation coefficients as well as wave forces on the vertical wall in terms of wave energy dissipation and wave loads by the proposed submerged flexible porous breakwater will be helpful for better understanding in the design of various types of protecting infrastructures for marine operations.

Water wave scattering by a submerged horizontal semi-circular barrier in a two-layer fluid

This paper studies water wave scattering by a submerged horizontal bottom-mounted semi-circular barrier in a two-layer fluid based on the linear potential flow theory. Two cases, i.e., the barrier completely located in the lower fluid and the barrier intersecting the internal interface of two fluids, are considered independently. In the former case, an analytical solution for the wave scattering problem is derived by using the method of multipole expansions, whereas in the latter case, a numerical solution is developed by employing the multi-domain boundary element method (BEM). The conservation relations of wave energy for incident waves in surface and internal modes are derived, and the reflection and transmission energy coefficients as well as the wave forces acting on the barrier are calculated. Comparisons of calculated results between the analytical solution and the BEM numerical solution are conducted. Calculation examples are presented to examine the effects of various physical parameters (angular frequency, interface position, density ratio of two fluids and incident angle) on wave forces acting on the barrier and the energy transformation between two wave modes. It is found that there is a noticeable difference in the energy transformation between the barrier intersecting the fluid interface or not.

Hydrodynamic performance of a submerged horizontal porous wave barrier

The hydrodynamic performance of a submerged horizontal porous wave barrier has been investigated using the analytical, numerical, and experimental approaches in regular waves. First, the analytical model is formulated under the assumption of potential flow and is solved using the matched eigenfunction expansion method, where the presence of porous barrier is recognized by the equivalent linearized quadratic energy dissipation model. Second, the numerical simulations based on a 3D Reynolds Averaged Navier–Stokes equations with standard low-Re k-ε turbulent closure and volume of fluid approach are accomplished. Analytical and numerical results are compared with the experimental results conducted in a 2D wave tank. The estimation of reflection, transmission, energy loss coefficients, and vertical force on the submerged horizontal barrier is shown to be in satisfactory agreement. From the systematic parametric study, it is found that there exists an optimal porosity and submergence depth near P=0.1,d/h=0.05−0.1 within the range of 2.0<kh<8.0. The analytical formulations using quadratic energy dissipation model predicts well the physical flow behavior up to H/λ=0.02. The present methodology imparts a reciprocative approach by using analytical and numerical models for a better design of the submerged horizontal porous wave barrier.

Restoring groundwater levels after tunneling: a numerical simulation approach to tunnel sealing decision-making

Tunneling is often unpopular with local residents and environmentalists, and can cause aquifer damage. Tunnel sealing is sometimes used to avoid groundwater leakage into the tunnel, thereby mitigating the damage. Due to the high cost of sealing operations, a detailed hydrogeological investigation should be conducted as part of the tunneling project to determine the impact of sealing, and groundwater modeling is an accurate method that can aid decision-making. Groundwater-level drawdown induced by the construction of the Headrace water-conveyance tunnel in Sri Lanka dried up 456 wells. Due to resulting socio-environmental problems, tunnel sealing was decided as a remedy solution. However, due to the expectation of significant delays and high costs of sealing, and because the water pressure in the tunnel may prevent groundwater seepage into the tunnel during operation, there was another (counter) decision that the tunnel could remain unsealed. This paper describes groundwater modeling carried out using MODFLOW to determine which option—sealed or unsealed tunnel—is more effective in groundwater level recovery. The Horizontal Flow Barrier and River packages of MODFLOW were used to simulate sealed and unsealed tunnels, respectively. The simulation results showed that only through tunnel sealing can the groundwater level be raised to preexisting levels after 18 years throughout the study area. If the tunnel remains unsealed, about 1 million m3/year of water conveyed by the tunnel will seep into the aquifer, reducing the operational capacity of the tunnel as a transport scheme. In conclusion, partial tunnel sealing in high-impact sections is recommended.

Experimental and numerical investigation of a surface-fixed horizontal porous wave barrier

Experimental and numerical investigations were conducted to study the performance of a surface-fixed horizontal porous wave barrier in regular waves. The characteristics of the reflection and transmission coefficients, energy dissipation, and vertical wave force were examined versus different porosities of the barrier. Numerical simulations based on 3D Reynolds Averaged Navier-Stokes equations with standard low-Re k-e turbulent closure and volume of fluid approach were accomplished and compared with the experimental results conducted in a 2D wave tank. Experimental measurements and numerical simulations were shown to be in satisfactory agreement. The qualitative wave behavior propagating over a horizontal porous barrier such as wave run-up, wave breaking, air entrapment, jet flow, and vortex generation was reproduced by CFD computation. Through the discrete harmonic decomposition of the vertical wave force on a wave barrier, the nonlinear characteristics were revealed quantitatively. It was concluded that the surface-fixed horizontal barrier is more effective in dissipating wave energy in the short wave period region and more energy conversion was observed from the first harmonic to higher harmonics with the increase of porosity. The present numerical approach will provide a predictive tool for an accurate and efficient design of the surface-fixed horizontal porous wave barrier.

Assessment of assumptions underlying models of prokaryotic pangenome evolution

BackgroundThe genomes of bacteria and archaea evolve by extensive loss and gain of genes which, for any group of related prokaryotic genomes, result in the formation of a pangenome with the universal, asymmetrical U-shaped distribution of gene commonality. However, the evolutionary factors that define the specific shape of this distribution are not thoroughly understood.ResultsWe investigate the fit of simple models of genome evolution to the empirically observed gene commonality distributions and genome intersections for 33 groups of closely related bacterial genomes. A model with an infinite external gene pool available for gene acquisition and constant genome size (IGP-CGS model), and two gene turnover rates, one for slow- and the other one for fast-evolving genes, allows two approaches to estimate the parameters for gene content dynamics. One is by fitting the model prediction to the distribution of the number of genes shared by precisely k genomes (gene commonality distribution) and another by analyzing the distribution of the number of genes common for k genome sets (k-cores). Both approaches produce a comparable overall quality of fit, although the former significantly overestimates the number of the universally conserved genes, while the latter overestimates the number of singletons. We further explore the effect of dropping each of the assumptions of the IGP-CGS model on the fit to the gene commonality distributions and show that models with either a finite gene pool or unequal rates of gene loss and gain (greater gene loss rate) eliminate the overestimate of the number of singletons or the core genome size.ConclusionsWe examine the assumptions that are usually adopted for modeling the evolution of the U-shaped gene commonality distributions in prokaryote genomes, namely, those of infinitely many genes and constant genome size. The combined analysis of genome intersections and gene commonality suggests that at least one of these assumptions is invalid. The violation of both these assumptions reflects the limited ability of prokaryotes to gain new genes. This limitation seems to stem, at least partly, from the horizontal gene transfer barrier, i.e., the cost of accommodation of foreign genes by prokaryotes. Further development of models taking into account the complexity of microbial evolution is necessary for an improved understanding of the evolution of prokaryotes.

Review of Modeling Approaches to Groundwater Flow in Deformed Carbonate Aquifers.

We discuss techniques to represent groundwater flow in carbonate aquifers using the three existing modeling approaches: equivalent porous medium, conduit network, and discrete fracture network. Fractures in faulted stratigraphic successions are characterized by dominant sets of sub-vertical joints. Grid rotation is recommended using the equivalent porous medium to match higher hydraulic conductivity with the dominant orientation of the joints. Modeling carbonate faults with throws greater than approximately 100 m is more challenging. Such faults are characterized by combined conduit-barrier behavior. The barrier behavior can be modeled using the Horizontal Flow Barrier Package with a low-permeability vertical barrier inserted to represent the impediment of horizontal flow in faults characterized by sharp drops of the piezometric surface. Cavities can occur parallel to the strike of normal faults generating channels for the groundwater. In this case, flow models need to account for turbulence using a conduit network approach. Channels need to be embedded in an equivalent porous medium due to cavities a few centimeters large, which are present in carbonate aquifers even in areas characterized by low hydraulic gradients. Discrete fracture network modeling enables representation of individual rock discontinuities in three dimensions. This approach is used in non-heavily karstified aquifers at industrial sites and was recently combined with the equivalent porous medium to simulate diffusivity in the matrix. Following this review, we recommend that the future research combines three practiced modeling approaches: equivalent porous medium, discrete fracture network, and conduit network, in order to capture structural and flow aspects in the modeling of groundwater in carbonate rocks.