Water Flow Model Research Articles

Effects of geometry on the thermal-fluid dynamics of a transition water flow in a square heated asymmetric cavity channel

Cavity channels candidate as a promising passive cooling configuration in thermal management. The superior performance and the limited scientific investigation regarding cavity channels dictate the need for better understanding of the thermal-fluid dynamics. To the best knowledge of the authors, there are few works investigating the effect of the geometry on the thermal-fluid dynamics of cavity channels with no focus on thermal data and quantification of the recirculation effect which can depend on geometrical parameters. In this work, numerical RANS simulations were performed in ANSYS Fluent with the kT-kL-ω turbulence model for a transitional water flow in an asymmetric square cavity channel with different attachment angles in the range of 20°–90°. The optimal heat transfer performance in terms of uniform and high Nusselt numbers is provided by the configuration with an angle in the range of approximately 30°–40°. Recirculation plays an important role in enhancing heat transfer and its strength is associated with the penetration capacity in the boundary layer. Two regions of strong and weak recirculation are identified along the channel; when the recirculation is weak, other mechanisms such as the thinning of the thermal boundary layer play a more important role. By increasing the attachment angle it is shown the local thinning of the thermal boundary layer, an augmentation of the magnitude of the negative turbulent heat flux peak, and a similarity of the conductive sublayer thickness with the recirculation strength.

Macropore flow in relation to the geometry and topology of soil macropore networks: Re-visiting the kinematic wave equation

The rapid flow of water through soil macropores significantly affects the partitioning of precipitation between surface runoff and infiltration and also the rate of solute transport in soil, both of which have an impact on the risk of contamination of surface water and groundwater. The kinematic wave equation is often employed as a model of gravity-driven water flow through soil macropores. The exponent in this simple model influences the pore water velocity attained in the macropores at any given input rate and is usually estimated by inverse modelling against measured flow rates or water contents. In theory, the exponent in the kinematic wave equation should depend on the geometry and topology of the conducting macropore networks, although these relationships have not so far been investigated. In this study, we related metrics of soil structure derived from X-ray images to values of the kinematic exponent estimated from drainage experiments on twenty-two columns sampled at three different field sites under two contrasting land uses and at three different depths.We found that smaller values of the exponent in the kinematic wave equation, which would equate to more rapid flow of water through soil macropores, were found in plough pan and subsoil columns of smaller macroporosity, for which biopores comprised a significant fraction. The macroporosity in these columns was more vertically oriented and poorly inter-connected, though still continuous across the sample. In contrast, topsoil columns from both arable land and grassland had better connected, denser and more isotropically-distributed macropore networks and larger values of the kinematic exponent. Our results suggest that for predictive modelling at large scales, it may be feasible to estimate the kinematic exponent using class pedotransfer functions based on pedological information such as land use and horizon type.

Study of the efficiency of machine learning algorithms based on data of various rocks

Background: Absolute permeability plays an important role in studying the fluids flow in porous media during the development of oil and gas reservoirs, the injection of CO2 into reservoirs for storage, the monitoring of pollutants migration in underground aquifers, and the modeling of catalytic systems. Therefore, an accurate and fast evaluation of its values is an important task.
 Aim: The purpose of this article is to study the applicability of machine learning methods for predicting the absolute permeability of carbonate samples, as well as ways to improve the prediction of permeability.
 Materials and methods: The input data is 408 small volumes extracted from four cylindrical carbonate samples composed almost entirely of calcite. Input data includes total and connected porosity, specific surface area, radii of all and only connected pores, coordination number, throat radius and length, tortuosity, and absolute permeability. Permeability prediction is carried out using regression machine learning methods such as random forest, extremely random trees and extended gradient boosting. Parameters (data) of small volumes were determined using pore-scale modeling of water flow in their pore space applying a specialized Avizo software.
 Results: Data of small volumes extracted from fractured and non-fractured samples were analyzed, and the results showed that there are good relationships between many parameters of small volumes. For example, the connected and total porosity have a second-order polynomial relationship with a high correlation coefficient. Using the above-mentioned regression machine learning methods, absolute permeability values were predicted when input data divided into training and testing data in a ratio of 80/20 and 70/30.
 Conclusion: Using the logarithm of permeability instead of permeability itself and considering fractured and non-fractured samples separately, can increase the accuracy of absolute permeability prediction using the above-mentioned machine learning methods up to 90%. The extremely random trees method is the most accurate among the three machine learning methods considered for our task.

A fully coupled dynamic water-mooring line system: Numerical implementation and applications

Several numerical challenges exist in the analysis of water-mooring line systems which require robust, yet practical, methods to address this type of fully coupled nonlinear dynamic problems. The present study proposes a novel class of numerical techniques for the formulation and implementation of a fully coupled dynamic system which involves water flows and catenary mooring line system. In particular, the three-dimensional water flow model is replaced by a simplified multilayer shallow water system with mass exchange terms between the layers including frictional forces at the bed topography and wind-driven forces at the water free-surface. Coupling conditions between the multilayer shallow water model and the mooring line system are also investigated in the current work. As numerical solvers we implement a fast finite volume method for the multilayer shallow water equations and a nonlinear dynamic analysis for the mooring line based on elastic catenary cable elements. Efficient calculations of the interaction forces between the shallow water flow and the submerged mooring line system and associated numerical implementations are also discussed. The accuracy and computational advantages of the proposed fully coupled system are verified using a series of well-established benchmark problems and wind-driven flows over both flat and non-flat beds. The computational results obtained show high performance the developed model and demonstrate the ability of the method to simulate fully coupled dynamic water-mooring line systems.

A Numerical Simulation of Moisture Reduction in Fine Soil Subgrade with Wicking Geotextiles.

A new wicking geotextile is proposed to control the water content of fine-grained soil subgrade. By comparing the spatial distribution of volumetric water content and matric suction before and after the installation of the wicking geotextile, the effectiveness of the geotextile in controlling the subgrade humidity is evaluated. Firstly, the hydraulic parameters of the wicking geotextile are obtained through laboratory tests using a pressure plate apparatus. Then, a numerical model for water flow in the subgrade is established using COMSOL to obtain the spatial distribution characteristics of humidity in the subgrade under different groundwater levels (2~8 m). The results show the wicking geotextile exhibits strong hydrophilicity, low water retention, and high horizontal permeability. Compared to the subgrade without geotextile, the water content of the soil above the geotextile decreases significantly by 7.6~9.6% at groundwater levels of 4~8m, while the saturation decreases by 18.3~23.0%, and the matric suction increases by 2~2.3 times. The wicking fabric functions as an effective drainage material to serve as a capillary barrier in the cross-plane direction and an effective drainage tunnel to transport water in the in-plane direction. The dynamic resilient modulus of the subgrade increases by 23.2~43.6%. The wicking geotextile effectively absorbs and drains weakly bound water in unsaturated soil due to the matric suction difference and its horizontal drainage capacity to improve the bearing capacity of the subgrade. It suggests that using wicking geotextile for drainage and reinforcement in fine-grained soil subgrades with groundwater levels ranging from 4 to 8 m is beneficial.

Assessment of the Impacts of Rainfall Characteristics and Land Use Pattern on Runoff Accumulation in the Hulu River Basin, China

Climate change causes the river basin water cycle disorders, and rainfall characteristics frequently result in flood disasters. This study aims to simulate and assess the response behavior of basin floods under the influence of rainfall characteristics and land use changes in the Hulu River basin using a 2D hydrological and hydraulic GAST (GPU Accelerated Surface Water Flow and Transport Model). The peak flow rate and water depth during floods were examined by simulating the evolution process of basin floods and related hydraulic elements under the independent effects of various rainfall characteristics or land use and further simulating the response results of basin floods under the combined effects of rainfall characteristics and land use. The seven scenarios were set to quantify the degree of influence that land use and rainfall characteristics have on the basin flood process based on examining changes in land use and rainfall characteristics in the research area. The results from different rainfall characteristics scenarios depicted that as the rainfall return period is shorter, the peak flow rate is higher, and the peak flow rate is lower as the return period is prolonged. Under different rainfall characteristics, the peak flow rate in scenario R8 is 41.30%, 40.00%, and 34.51% higher than the uniform distribution of rainfall, while water depth is decreased by 0.55%, increased by 4.96% and 2.92% as compared to the uniform distribution of rainfall. While under different land use scenarios, it is observed that the change in land use has increased 2.7% in cultivated land and 1.1% in woodland. In addition, the interactive effect of different rainfall characteristics and land use it can be seen that the scenario with the greatest reduction in flood risk due to rainfall characteristics and land use is RL2-4, representing a 12.55% decrease in peak flow and a 37.69% decrease in peak water depth. In this scenario, the rainfall is heavier in the southeast and northwest regions and lighter in the northeast and southwest regions. The land use type is characterized by reforestation and the return of cultivated land to forests. The changes in rainfall distribution and the increase in grassland contribute to the decrease in flood threat. Future research in the erodible parts of the Hulu River basin, planning for water resources, and soil and water conservation can all benefit from the study’s conclusions.

Root distributions predict shrub–steppe responses to precipitation intensity

Abstract. Precipitation events are becoming more intense around the world, changing the way water moves through soils and plants. Plant rooting strategies that sustain water uptake under these conditions are likely to become more abundant (e.g., shrub encroachment). Yet, it remains difficult to predict species responses to climate change because we typically do not know where active roots are located or how much water they absorb. Here, we applied a water tracer experiment to describe forb, grass, and shrub root distributions. These measurements were made in 8 m by 8 m field shelters with low or high precipitation intensity. We used tracer uptake data in a soil water flow model to estimate how much water respective plant root tissues absorb over time. In low-precipitation-intensity plots, deep shrub roots were estimated to absorb the most water (93 mm yr−1) and shrubs had the greatest aboveground cover (27 %). Grass root distributions were estimated to absorb an intermediate amount of water (80 mm yr−1) and grasses had intermediate aboveground cover (18 %). Forb root distributions were estimated to absorb the least water (79 mm yr−1) and had the least aboveground cover (12 %). In high-precipitation-intensity plots, shrub and forb root distributions changed in ways that increased their water uptake relative to grasses, predicting the increased aboveground growth of shrubs and forbs in these plots. In short, water uptake caused by different rooting distributions predicted plant aboveground cover. Our results suggest that detailed descriptions of active plant root distributions can predict plant growth responses to climate change in arid and semi-arid ecosystems.

Using specularity content to evaluate eight geothermal heat flow maps of Totten Glacier

Abstract. Geothermal heat flow (GHF) is the dominant factor affecting the basal thermal regime of ice sheet dynamics. But it is poorly defined for the Antarctic ice sheet. We compare the basal thermal state of the Totten Glacier catchment as simulated by eight different GHF datasets. We use a basal energy and water flow model coupled with a 3D full-Stokes ice dynamics model to estimate the basal temperature, basal friction heat and basal melting rate. In addition to the location of subglacial lakes, we use specularity content of the airborne radar returns as a two-sided constraint to discriminate between local wet or dry basal conditions and compare the returns with the basal state simulations with different GHFs. Two medium magnitude GHF distribution maps derived from seismic modelling rank well at simulating both cold- and warm-bed regions, the GHFs from Shen et al. (2020) and Shapiro and Ritzwoller (2004). The best-fit simulated result shows that most of the inland bed area is frozen. Only the central inland subglacial canyon, co-located with high specularity content, reaches the pressure melting point consistently in all the eight GHFs. Modelled basal melting rates in the slow-flowing region are generally 0–5 mm yr−1 but with local maxima of 10 mm yr−1 at the central inland subglacial canyon. The fast-flowing grounded glaciers close to the Totten ice shelf are lubricating their bases with meltwater at rates of 10–400 mm yr−1.

Numerical modelling of flow and heat transport in closed mines. Case study Walsum drainage province in the Ruhr coal-mining area

High geothermal potential and multiple mine-water-based geothermal installations in Germany and other countries improve the relevance of detailed studies and modeling of promising sites. In this context, we developed a numerical model of water flow and heat transport in the Walsum mine drainage province in the west of the Ruhr coalmining area using the available data on geology, mining, water levels, pumping, and the temperatures of deep rocks and mine water. The model was validated by varying the parameters of groundwater recharge and hydraulic conductivity to achieve sufficient consistency with measured inflows and pumping rates from the central pumping facility located in the Walsum 2 shaft. The calculated mine water temperature of 30.3 ºC is close to the average of the measured temperature varied within the range of 29 – 33 ºC during the last years of mine maintenance. Using the numerical model, we evaluated the expected thermal capacity of a hypothetical open-loop circulation system and two closed-loop geothermal systems within the study area. The installation and operation of these systems would enable the generation of a thermal capacity from a few dozen kW to 1 MW sufficient for small-size to mid-size heat consumers with insignificant impact on the high thermal energy potential of the Walsum mine drainage province.

A Mechanistic Model for Simulation of Carbendazim and Chlorothalonil Transport through a Two-Stage Vertical Flow Constructed Wetland

A mechanistic model was developed to simulate one-dimensional pesticide transport in two-stage vertical flow constructed wetland. The two pesticides taken under study were carbendazim and chlorothalonil. The water flow patterns within the constructed wetland were simulated using the Richards equation. Water content and vertical flux, which are the outputs of the substrate water flow model, were used to calculate the substrate moisture-related parameters and advection term in the solute transport model. The governing solute transport equation took into account a total of six processes: advection, molecular diffusion, dispersion, adsorption to the solid surface, degradation and volatilization. A total of 14 simulation cases, corresponding with available experimental data, were used to calibrate the model, followed by further simulations with standardized influent pesticide concentrations. The simulations indicated that the constructed wetland reached a steady state of pesticide removal after 7 days of operation. Two distinct water flow patterns emerged under saturated and unsaturated conditions. The patterns observed while varying the hydraulic loading rates were similar for each individual saturation condition. Two-factor ANOVA of the simulated data further revealed that the carbendazim and chlorothalonil removal was dependent on the hydraulic loading rates, but it was independent of the influent pesticide concentration. Analysis of the simulated pesticide removal showed that degradation emerged as the predominant removal process over time for both the pesticides. The model developed in this study can be an important tool for the design and construction of treatment wetlands for pesticide removal from wastewater.

From global glacier modeling to catchment hydrology: bridging the gap with the WaSiM-OGGM coupling scheme

Coupled glacio-hydrological models have recently become a valuable method for predicting the hydrological response of catchments in mountainous regions under a changing climate. While hydrological models focus mostly on processes of the non-glacierized part of the catchment with a relatively simple glacier representation, the latest generation of standalone (global) glacier models tend to describe glacier processes more accurately by using new global datasets and explicitly modeling ice-flow dynamics. Yet, to the authors' knowledge, existing catchment-scale coupled glacio-hydrological models either do not include these most recent advances in glacier modeling or are simply not available to other users. By making use of the capabilities of the free, distributed, physically-based Water Flow and Balance Simulation Model (WaSiM) and the Open Global Glacier Model (OGGM), a coupling scheme is developed to bridge the gap between global glacier representation and local catchment hydrology. The WaSiM-OGGM coupling scheme is used to further assess the impacts under future climates on the glaciological and hydrological processes in the Gepatschalm catchment (Austria), by considering a combination of three climate projections under the Representative Concentration Pathways (RCP) 2.6, 4.5, and 8.5. Additionally, the results are compared to the original WaSiM model with the integrated Volume-Area (VA) scaling approach for modeling glaciers. Although both models (WaSiM with VA scaling and WaSiM-OGGM coupling scheme) perform very similar during the historical simulations (1971–2010), large discrepancies arise when looking into the future (2011–2100). In terms of runoff, the VA scaling model suggests a reduction of the mean monthly peak between 10–19%, whereas a reduction of 26–41% is computed by the coupling scheme. Similarly, results suggest that glaciers will continuously retreat until 2100. By the end of the century, between 20–43% of the 2010 glacier area will remain according to the VA scaling model, but only 1–23% is expected to remain with the coupling scheme. The results from the WaSiM-OGGM coupling scheme raises awareness of including more sophisticated glacier evolution models when performing hydrological simulations at the catchment scale in the future. As the WaSiM-OGGM coupling scheme is released as open-source software, it is accessible to any interested modeler with limited or even no glacier knowledge.

Three-dimensional X-ray imaging of macropore flow

Macropores are known to be important pathways for the rapid transport of water, solutes and colloids in soil. Nevertheless, we still know very little about how the topology and geometry of macropore networks govern water flow configurations and velocities in natural soil. In this study, we aimed at gaining more insight into macropore flow by using X-ray tomography to quantify air–water distributions in the macropore networks of undisturbed topsoil and subsoil columns of a clay soil at varying steady-state flow rates. We observed that while large fractions of the macropore network remained air-filled, the air phase only became entrapped when the irrigation rate was very close to the saturated hydraulic conductivity of the soil. The data enabled us to parameterize a kinematic wave model for water flow following the approach proposed in Jarvis et al. (2017a). Follow-up experiments would be required to evaluate whether these kinematic wave parameters derived by X-ray imaging match with those obtained from outflow measurements. We found that quantitative X-ray imaging of macropore flow through soils still remains a challenging task. We recommend that future experiments are conducted on smaller soil samples to improve image resolution and minimize experimental time spans as well as X-ray image noise and illumination bias. Such experiments could also include 3-D tracer imaging to identify the imaged macropore networks transporting most of the water (i.e. the backbone) at varying steady irrigation rates.

Modeling Water Flow in Variably Saturated Porous Soils and Alluvial Sediments

The sustainable exploitation of groundwater resources is a multifaceted and complex problem, which is controlled, among many other factors and processes, by water flow in porous soils and sediments. Modeling water flow in unsaturated, non-deformable porous media is commonly based on a partial differential equation, which translates the mass conservation principle into mathematical terms. Such an equation assumes that the variation of the volumetric water content (θ) in the medium is balanced by the net flux of water flow, i.e., the divergence of specific discharge, if source/sink terms are negligible. Specific discharge is in turn related to the matric potential (h), through the non-linear Darcy–Buckingham law. The resulting equation can be rewritten in different ways, in order to express it as a partial differential equation where a single physical quantity is considered to be a dependent variable. Namely, the most common instances are the Fokker–Planck Equation (for θ), and the Richards Equation (for h). The other two forms can be given for generalized matric flux potential (Φ) and for hydraulic conductivity (K). The latter two cases are shown to limit the non-linearity to multiplicative terms for an exponential K-to-h relationship. Different types of boundary conditions are examined for the four different formalisms. Moreover, remarks given on the physico-mathematical properties of the relationships between K, h, and θ could be useful for further theoretical and practical studies.

HAM-based analysis of nonlinear convection of two-layer water and air-TiO2 flow model in a vertical channel

This article investigates the nonlinear convection transport of heat transfer in a nanoparticle laden two-layer flow in a vertical channel under the additional effects of magneto hydrodynamic (MHD) and radiation. Water and compressed air are chosen as the fluids in two-layer flow. This type of flow is commonly seen in the transportation of oil and gas in pipelines, therefore, correct understanding of flow dynamics including convective transport phenomenon, layers interaction, and graphical interpretation of flow variables is important for better designing of equipments. Titanium dioxide ( Ti O 2 ) nanoparticles are added into the air to magnify the transport mechanism. The mathematical model incorporating buoyancy, nonlinear convection, MHD, viscous dissipation, and radiation effects is presented in the form of nonlinear Ordinary Differential Equation (ODE’s). Then utilizing homotopy analysis method, these interconnected differential equations are resolved. A visual illustration shows how different physical factors affect fluid rate and temperature. Tabulated values of Nusselt number at both walls are presented to discuss the numerical data of convection and found that it decreases toward the left boundary but rises toward the right wall by increasing Eckert number and Prandtl number. The results of this two-layer flow will help in understanding this complex flow which has practical applications in various engineering processes and further enhance the scientific knowledge related to multilayer problems.

A 1D-2D dynamic bidirectional coupling model for high-resolution simulation of urban water environments based on GPU acceleration techniques

Urban water pollution is an important issue of global concern. An integrated model that can be employed to efficiently and accurately simulate the whole system of urban surfaces, underground drainage systems and river and lake water environments is an important tool for comprehensively treating urban water pollution. Therefore, we developed a 1D-2D dynamic bidirectional coupling model for high-resolution simulation of the entire urban water environment (GAST-GCSE) comprising the graphics processing unit (GPU)-accelerated surface water flow and transport (GAST) model and the storm water management model (SWMM). The SWMM and GAST model were used for coupled simulation of hydrological and hydrodynamic quantities and water quality in 1D and 2D regions, respectively. The time steps of the two models are the same. Compute Unified Device Architecture (CUDA) parallel architecture programming and advanced model algorithms were adopted to improve the simulation accuracy and efficiency of the GAST model. The two models were dynamically coupled by calling the dynamic link library (DLL) in the SWMM and adding external functions. Two test cases were considered to verify the performance of the proposed model. The results showed that the proposed model provides high simulation accuracy. Compared with the results obtained by commercial software and theoretical values, the relative errors of the water quantity and water quality simulation results were less than 3%. The whole water environment of the main urban area of Changzhi city in North China was evaluated by the GAST-GCSE model, and the simulation results were compared with monitoring data and pure SWMM simulation results. The simulation errors at three waterlogging points were less than 14%, and the Nash efficiency coefficients (NSE) of the NH4+-N concentration in the three river sections were 0.77, 0.56 and 0.64, respectively. Notably, the model needed only 4.39 h to simulate rainfall production and confluence on a 2D surface, a river flood in complex terrain and its accompanying pollutant evolution process involving 3,849,428 cells over 4 h. This study provides new insights, ideas and tools for urban water environment simulation and change mechanism study.

Estimating cumulative evapotranspiration using superconducting gravimeter data: a study in Buenos Aires Province, Argentina

ABSTRACT A new method for estimating cumulative evapotranspiration using superconducting gravimeter data is presented. Evapotranspiration is computed from the classical water balance equation using precipitation data, estimated runoff, and water storage assessed from temporal variations of gravity. The gravity-based methodology is applied to a study site located at the Argentinian–German Geodetic Observatory (AGGO) in the flatland area of Buenos Aires Province (Argentina), which houses the only superconducting gravimeter in South America. For the period from April 2017 to May 2018, the cumulative evapotranspiration is estimated at 710.7 ± 25.7 mm. Evapotranspiration values are successfully compared with estimates provided by three methods, based on the water balance equation using soil moisture data, the numerical modelling of water flow in the unsaturated zone and the MODIS product (MOD16A2). The analysis of evapotranspiration also allows us to identify a period of soil water stress which is validated by soil moisture measurements.

Validating coupled flow theory for bare‐soil evaporation under different boundary conditions

AbstractEvaporation from bare soil is an important hydrological process and part of the water and energy balance of terrestrial systems. Modeling bare‐soil evaporation is challenging, mainly due to nonlinear couplings among liquid water, water vapor, and heat fluxes. Model concepts of varying complexity have been proposed for predicting evaporative water and energy fluxes. Our aim was to test a standard model of coupled water, vapor, and heat flow in the soil using data from laboratory evaporation experiments under different boundary conditions. We conducted evaporation experiments with a sand and a silt loam soil and with three different atmospheric boundary conditions: (i) wind, (ii) wind and short‐wave radiation, and (iii) wind and intermittent short‐wave radiation. The packed soil columns were closed at the bottom (no water flux) and instrumented with temperature sensors, tensiometers, and relative humidity probes. We simulated the evaporation experiments with a coupled water, vapor, and heat flow model, which solves the surface energy balance and predicts the evaporation rate. The evaporation dynamics were predicted very well, in particular the onset of stage‐two evaporation and the evaporation rates during the stage. A continuous slow decrease of the measured evaporation rate during stage‐one could not be described with a constant aerodynamic resistance. Adding established soil resistance parametrizations to the model significantly degraded model performance. The use of a boundary‐layer resistance, which takes into account the effect of point sources of moisture, improved the prediction of evaporation rates for the sandy soil, but not for the silt loam.

Development of a water flow and nutrients concentration simulation model for cascade hydroponic systems

Development of a water flow and nutrients concentration simulation model for cascade hydroponic systems

Approximate Solutions for A Fractional Shallow Water Flow Model

This paper presents the solutions of fractional Drinfeld-Sokolov-Wilson (DSW) equations
 that occur in shallow water flow models using the residual power series method.
 The fractional derivatives and integrals are considered in the conformable sense. In
 addition, surface plots of the solutions are given. The solutions and results show that
 the present method is very efficient and effective due to the lack of a need for complex
 calculations and that the method also has a wide range of practicability in the resolution
 of partial differential fractional equations.

Development of mechanistic-artificial intelligence model for simulation of numerical data of water flow in porous materials

Fluid dynamics of water flow through porous metallic media is significant for cooling and heating applications. The prediction of the velocity of fluid flowing inside the porous media could provide useful data for pressure drop calculations. This prediction is usually performed by a mathematical modeling approach like computational fluid dynamics (CFD). The computational method is a powerful means of precision calculation, although it could take more time and expenditure for more complex geometries or more complex fluid flow regimes. Artificial intelligence (AI) algorithms could learn and map CFD data under several conditions. AI algorithms could continue the prediction of physical data, saving computing time and performance. The present study focuses on CFD modeling of water flow inside a pipe filled with copper porous media. For the first time, the fuzzy bee algorithm (BAFIS) maps the generated CFD data to inlet velocities of 0.5, 0.7, 1.1, and 1.3 m/s. The BAFIS intelligence assessment for accurate convective flow prediction in porous media is not available in the literature. Additionally, the reliability of this method for missing CFD data prediction has not been considered yet. The results showed that the maximum intelligence (regression=0.97) is for a bee number of 140. Using the most intelligent BAFIS, the outlet velocity can be predicted by the artificial intelligence method for further nodes and inlet velocities without CFD modeling. BAFIS could precisely predict the outlet velocity for missing data with an inlet velocity of 0.91 m/s based on previously mapped data. A comparison was made between BAFIS and the fuzzy neural network (ANFIS) for CFD data predictions. The mean-square error and the root-mean-square error of ANFIS were slightly more than BAFIS (i.e., 0.1% and 3%, respectively).