Simulation Of Solute Transport Research Articles

A finite particle method (FPM) for Lagrangian simulation of conservative solute transport in heterogeneous porous media

When a smoothed particle hydrodynamics (SPH) method, a Lagrangian and meshfree numerical scheme, is used to solve the advection-dispersion equation (ADE), SPH solutions may not be accurate when particles are irregularly distributed, i.e., the distance between two neighbor particles varies irregularly in a simulation domain. Particle irregularity may be caused by nonuniform groundwater flow in a heterogeneous field of hydraulic conductivity. This study explores for the first time whether the Finite Particle Method (FPM) can provide more accurate ADE solutions than SPH does for irregularly distributed particles. FPM is similar to SPH in theory, but uses a modified kernel gradient to construct a SPH approximation of solute concentration gradients. Performance of SPH and FPM with irregularly distributed particles is evaluated by using two numerical cases. The first case considers only diffusive transport, and has an analytical solution for the evaluation. The second case considers both advection and dispersion, and uses a numerical solution as a reference for the evaluation. For each of the two cases, several numerical experiments are conducted using multiple sets of irregularly distributed particles with different levels of particle irregularity due to different levels of heterogeneity of hydraulic conductivity. Numerical results indicate that, for the numerical experiments of this study, FPM outperforms SPH to yield more accurate ADE solutions. However, FPM solutions are still subject to numerical errors, and the errors increase when the level of heterogeneity of hydraulic conductivity increases. Further improvement of FPM is warranted in a future study.

Effects of impurities H2S and N2 on CO2 migration and dissolution in sedimentary geothermal reservoirs

CO2 geological storage is a key technology for reducing greenhouse gas emissions. But merely injecting CO2 into sedimentary strata has high costs and no direct economic benefit. A more cost-effective and environment-friendly way is storing and utilizing CO2 in deep sedimentary reservoirs for heat extraction. However, the effects of ubiquitous impurities on the sequestration safety and heat extraction performance of CO2 in high-temperature reservoirs, are unclear. In this study, we have performed numerical simulations of multiphase flow, solute transport, and heat conduction at reservoir temperature 50–250 °C, to reveal the effects of two common impurities (i.e., H2S and N2) on the migration, dissolution and thermal absorption of CO2 in sedimentary geothermal reservoirs. It is found that, for higher reservoir temperature, the impurity N2 leads to smaller promotion on the migration of gas phase, but retains the acceleration effect on the cooling in reservoir. Above 150 °C, the deviations on the gas phase migration velocity and the CO2 dissolution efficiency caused by 10 mol% N2 are less than 12%. Meanwhile, the 10 mol% N2 can increase the total heat absorbed from reservoir by about 40% at least. It means that N2 has the potential to benefit a CO2 plume geothermal system significantly without distinct risks. Besides, it is also found that the CO2 dissolution efficiency can be increased and decreased by H2S and N2 through partitioning phenomenon, respectively.

Using Complex Probability Amplitudes to Simulate Solute Transport in Composite Porous Media

AbstractProbability amplitudes are fundamental to quantum mechanics and offer robust descriptions of complicated systems, which have allowed physicists to explain behaviors inaccessible to classical physics. This article ponders how some of the same conceptual underpinnings of the mathematics used for modeling quantum systems might be applied to subsurface water resources problems and speculates how these tools could facilitate applications on quantum computers. A probability amplitude‐based model for describing advective‐dispersive transport in porous media using linear operators is investigated. The proposed complex valued model decomposes spreading into two “sub‐continuum partial dispersion” coefficients, and this recovers classical spreading when the sum of these coefficients is the Fickian dispersion coefficient. However, the probability amplitudes have a many‐to‐one relationship to a probability distribution, so it embeds a level of heterogeneity into seemingly equivalent functions. Two propagators with different sub‐continuum coefficients may have the same macroscopic behavior when either is considered in isolation, but when they act on the other the system’s behavior changes. Differences in the amplitudes cause a reduction in spreading as velocity correlations are disrupted, despite both propagators having identical dispersion coefficients, and this cannot be achieved using classical methods without changing the dispersion coefficient. The main point is that these amplitude‐based models offer a way to embed information about the system into the propagators, instead of just “averaging it out” when making an upscaled model.

Simulation of solute transport behaviors in saturated karst aquifer system

The karst development makes aquifer have strong anisotropy and heterogeneity. In order to reveal the characteristics of solute transport in the karst fissure–conduit aquifer system, this study presents a physical model of fissure–conduit in laboratory experiments to carry out the solute transport simulation. In this paper, the tracer tests of fissure–conduit combination, fissure, and conduit solute transport process in saturated flow are designed. We found that different aquifer structures and tracer injection points have an influence on the shape of the breakthrough curve. Besides, the two-dimensional dispersion model of tracer injection of the instantaneous point was used to calculate the dispersion parameters of each group of experiments. Then, the dynamic responses of the linear distance (x) between the injection point and the receiving point, initial time (t0), peak time (tm), peak concentration (cm), average tracer transport velocity (V), and porosity (p) of aqueous media to the longitudinal dispersion coefficient are discussed. In addition, according to the measured data, Gaussian multi-peak fitting can be used to reflect the overall shape and change trend of the multi-peak BTC. These results demonstrate the solute transport behaviors in the saturated karst aquifer system, which have important reference significance for solving the engineering environmental problems in the karst area.

Efficient simulation of groundwater solute transport using the multipoint flux approximation method with arbitrary polygon grids

Groundwater solute transport models are critical for simulating subsurface contaminant processes, guiding policy making for groundwater management, and pollution remediation. In this study, we develop an improved groundwater solute transport model which uses arbitrary-shaped polygons for spatial discretization. This allows for the efficient simulation of solutes in complex boundary aquifer systems. The governing equation is discretized in space using the finite volume method to ensure flux conservation both locally and globally. The dispersion term is discretized by the multipoint flux approximation (MPFA) method which is a finite volume method that is suitable on arbitrary polygon grids. A total variation diminishing (TVD) scheme is used to control numerical oscillation in the advection term. The functionality and performance of the new model are demonstrated by comparing the simulated results between our model and MT3DMS model in five case studies. The comparisons show that the new model can well represent the solute transport processes. In a real-world watershed with complex-shaped boundaries, the new model outperforms the conventional groundwater solute transport model with smaller computed root-mean-square errors. These modeling results suggest that our model can provide robust modeling solutions for simulating groundwater solute transport processes, especially in aquifer systems with complex shaped boundaries. Furthermore, our model can provide a more flexible discretization solution for coupling surface water models and groundwater solute transport models.

Numerical Simulation of Solute Transport in Saturated Porous Media with Bounded Domains.

This work presents the first attempt to develop unconditionally stable, implicit finite difference solutions of one-sided spatial fractional advection-dispersion equation (s-FADE) by imposing the nonzero Dirichlet boundary condition (ND BC) or the nonzero fractional Robin boundary condition (NFR BC) at inlet boundary and the zero fractional Neumann boundary condition (ZFN BC) at outlet boundary. The results of the numerical studies performed using artificial solute transport parameters demonstrated that the numerical solution with the NFR BC as the inlet boundary produced much more realistic concentration values. The numerical solution with the NFR BC at the inlet boundary was capable of correctly describing the Fickian and non-Fickian behaviors of the solute transport at different α values, and it had the relatively same accuracy at different numbers of the spatial nodes. Also, the practical application of the numerical solution with the NFR BC as the inlet boundary was investigated by conducting tracer experiments in homogeneous and heterogeneous soil columns. According to the obtained results, this numerical solution described well solute transport in the homogenous and heterogeneous soils. The α values of the homogeneous and heterogeneous soils were obtained in the ranges of 1.849-1.999 and 1.248-1.570, respectively, which were in excellent agreement with the physical properties of the soils. In a nutshell, the numerical solution of the s-FADE with the NFR BC as the inlet boundary can be successfully applied to describe the solute transport in the homogeneous and heterogeneous soils with bounded spatial domains.

Assessing the Impact of Lining Polluted Streams on Groundwater Quality: A Case Study of the Eastern Nile Delta Aquifer, Egypt

Groundwater is considered to be an important water supply for domestics, industry, and irrigation in many areas of the world. Renewable groundwater is recharged by rainfall and seepage from canals and open drain networks. Agricultural and industrial drainage, as well as domestic drainage, represent the main discharges into open drains. Therefore, these drains are considered to be a source of recharge as well as a source of pollution. In this study, we aim to evaluate the impact of the Bahr El Baqar drain system on groundwater quality in the Eastern Nile Delta, Egypt. MODFLOW was used to create a numerical model to simulate groundwater flow in an aquifer and MT3DS was used to simulate solute transport from the open contaminated Bahr El Baqar drain to the groundwater. Two approaches were developed in the study area. The first approach was applied to investigate the impact of increasing the abstraction rates on the contaminant transport into the aquifer, the second approach was developed to identify the effect of lining the drain using different materials on contaminant extension in the aquifer to protect groundwater quality in the east Nile Delta Aquifer. The results showed that the TDS values increased by 18.23%, 23.29%, and 19.24% with increased abstraction rates of 15%, 34%, and 70%, resulting from population increases in 2010, 2025, and 2040, respectively; however, the TDS in the aquifer decreased from 0.6%, to 6.36%, 88.35%, and 90.47% by using lining materials.

Simultaneous identification of groundwater pollution source spatial–temporal characteristics and hydraulic parameters based on deep regularization neural network-hybrid heuristic algorithm

The prerequisites for groundwater pollution remediation include providing accurate identification of hydraulic parameters of the contaminated site and retrospecting pollution source release history from sparse and discrete observation data. However, expensive time running and high nonlinearity of groundwater solute transport simulation model hinder precise and fast identification. To overcome expensive running cost, a deep regularization neural network (DRNN) was introduced to emulate the highly nonlinear simulation model. A hybrid heuristic algorithm, integrating a particle swarm optimization (PSO) and a differential evolution markov chain (DEMC) method was employed to provide reliable identification results of pollution source characteristics and hydraulic parameters. In particular, DEMC conducted a global search and PSO performed a local search. To evaluate transferability for various groundwater pollution cases, the proposed method was utilized to provide an identification for a point source pollution case and a non-point source pollution case. The results indicated that DRNN was capable of emulating input–output relationship of a highly nonlinear and high time expense simulation model of two cases with a R2 of 0.997 and 0.984, respectively, accelerating identification process. Furthermore, when the PSO was combined with DEMC, the optimization was accelerated while escaping from the local optimum, substantially improving the efficiency and accuracy of the identification task.

WITHDRAWN: Efficient simulation of groundwater solute transport using the multipoint flux approximation method with arbitrary polygon grids

WITHDRAWN: Efficient simulation of groundwater solute transport using the multipoint flux approximation method with arbitrary polygon grids

Topographic features of nano-pores within the osteochondral interface and their effects on transport properties –a 3D imaging and modeling study

Recent insights suggest that the osteochondral interface plays a central role in maintaining healthy articulating joints. Uncovering the underlying transport mechanisms is key to the understanding of the cross-talk between articular cartilage and subchondral bone. Here, we describe the mechanisms that facilitate transport at the osteochondral interface. Using scanning electron microscopy (SEM), we found a continuous transition of mineralization architecture from the non-calcified cartilage towards the calcified cartilage. This refurbishes the classical picture of the so-called tidemark; a well-defined discontinuity at the osteochondral interface.Using focused-ion-beam SEM (FIB-SEM) on one osteochondral plug derived from a human cadaveric knee, we elucidated that the pore structure gradually varies from the calcified cartilage towards the subchondral bone plate. We identified nano-pores with radius of 10.71 ± 6.45 nm in calcified cartilage to 39.1 ± 26.17 nm in the subchondral bone plate.The extracted pore sizes were used to construct 3D pore-scale numerical models to explore the effect of pore sizes and connectivity among different pores. Results indicated that connectivity of nano-pores in calcified cartilage is highly compromised compared to the subchondral bone plate. Flow simulations showed a permeability decrease by about 2000-fold and solute transport simulations using a tracer (iodixanol, 1.5 kDa with a free diffusivity of 2.5 × 10−10 m2/s) showed diffusivity decrease by a factor of 1.5. Taken together, architecture of the nano-pores and the complex mineralization pattern in the osteochondral interface considerably impacts the cross-talk between cartilage and bone.

Simulation of reactive solute transport in the critical zone: a Lagrangian model for transient flow and preferential transport

Abstract. We present a method to simulate fluid flow with reactive solute transport in structured, partially saturated soils using a Lagrangian perspective. In this context, we extend the scope of the Lagrangian Soil Water and Solute Transport Model (LAST) (Sternagel et al., 2019) by implementing vertically variable, non-linear sorption and first-order degradation processes during transport of reactive substances through a partially saturated soil matrix and macropores. For sorption, we develop an explicit mass transfer approach based on Freundlich isotherms because the common method of using a retardation factor is not applicable in the particle-based approach of LAST. The reactive transport method is tested against data of plot- and field-scale irrigation experiments with the herbicides isoproturon and flufenacet at different flow conditions over various periods. Simulations with HYDRUS 1-D serve as an additional benchmark. At the plot scale, both models show equal performance at a matrix-flow-dominated site, but LAST better matches indicators of preferential flow at a macropore-flow-dominated site. Furthermore, LAST successfully simulates the effects of adsorption and degradation on the breakthrough behaviour of flufenacet with preferential leaching and remobilization. The results demonstrate the feasibility of the method to simulate reactive solute transport in a Lagrangian framework and highlight the advantage of the particle-based approach and the structural macropore domain to simulate solute transport as well as to cope with preferential bypassing of topsoil and subsequent re-infiltration into the subsoil matrix.

Adaptation and validation of the ParSWMS numerical code for simulation of water flow and solute transport in soilless greenhouse substrates

Numerical simulation of three-dimensional water flow and solute transport in containerized variably saturated soilless substrates with complex hydraulic properties and boundary conditions necessitates high-resolution discretization of the spatial and temporal domains, which commonly leads to several million nodes requiring numerical evaluation. Even today’s computing prowess of workstations is not adequate to tackle such problems within a reasonable timeframe, especially when numerous realizations are required to optimize the geometry, substrate properties, and irrigation and fertigation management of soilless plant growth modules. Hence, the parallelization of the numerical code and utilization of high performance computing (HPC) are essential. Here, we adapted and applied the ParSWMS parallelized code that is amenable to solving the 3D Richards equation for water flow and the convection-dispersion equation for solute transport subject to linear solute adsorption. The code was modified to allow for nonlinear equilibrium solute adsorption with new boundary conditions and applied to simulate water flow and nitrogen and phosphorus transport in containerized soilless substrates. Multi-solute transport simulations with the modified Linux ParSWMS code were first performed on a workstation and referenced to the Windows-based HYDRUS (2D/3D) numerical code. After confirming the agreement between the modified ParSWMS code and HYDRUS (2D/3D), various preconditioners and iterative solvers were evaluated to find the computationally most efficient combinations. The performance of the modified ParSWMS code and its stability were compared to HYDRUS (2D/3D) simulations for three soilless substrates consisting of horticultural perlite, volcanic tuff, and a volcanic tuff/coconut coir mixture. Considering the solute mass balance error as a stability measure, ParSWMS outperformed HYDRUS (2D/3D). Moreover, simulations with the modified ParSWMS code were about 22% faster than simulations with HYDRUS (2D/3D) on the workstation. Tests of the modified ParSWMS on two HPC clusters with 28 and 94 cores revealed a potential computational speedup of 94% relative to the HYDRUS (2D/3D) simulations performed on the workstation.

A Spatially Distributed, Physically-Based Modeling Approach for Estimating Agricultural Nitrate Leaching to Groundwater

Nitrogen-nitrate, while being fundamental for crop production, is of particular concern in the agricultural sector, as it can easily leach to the water table, worsening groundwater quality. Numerical models and Geographic Information System may support the estimation of nitrate leaching rates in space and time, to support sustainable agricultural management practices. In this paper, we present a module for the simulation of the processes involved in the nitrogen cycle in the unsaturated zone, including nitrate leaching. This module was developed taking steps from the ANIMO and EPIC model frameworks and coupled to the hydrological models integrated within the FREEWAT platform. As such, the nitrogen cycle module was then included in the FREEWAT platform. The developed module and the coupling approach were tested using a simple synthetic application, where we simulated nitrate leaching through the unsaturated zone for a sunflower crop irrigated district during a dry year. The results of the simulation allow the estimation of daily nitrate concentration values at the water table. These spatially distributed values may then be further used as input concentration in models for simulating solute transport in aquifers.

New Finite Volume–Multiscale Finite-Element Model for Solving Solute Transport Problems in Porous Media

AbstractThis paper proposes a finite volume–Yeh’s finite-element–multiscale finite-element model (FVYMSFEM) for solute transport simulations, which is extended from its elliptic groundwater-flow eq...

How to relevantly characterize hydraulic properties of saline and sodic soils for water and solute transport simulations

Proper sets of soil hydraulic parameters are important for modeling water flow and solute transport. Using field saline water and laboratory tap water, this study systematically analyzed and compared laboratory saturated hydraulic conductivity, Kls and soil water retention curve, SWRC of saline soil within three soil multi-layered profiles to find out the adequate initial hydraulic parameter set of a hydrological model. Falling and constant head methods were applied on three undisturbed samples per pedogenic horizon. SWRC were established on the same samples, using both saline and tap water. Results showed the spatial variability of hydraulic properties in vertical and horizontal dimensions when comparing soil profiles' properties. The constant head method with saline water resulted in the highest Kls values followed by the falling head method with saline water, the constant head method with fresh water, and the falling head method with fresh water. The differences in Kls were up to two orders of magnitude when comparing all datasets. Fresh water reduced soil water retention significantly. A significant correlation between the hydraulic properties and some physical and chemical properties suggests that the effects of soil texture, structure, calcium carbonate and gypsum content, and CEC (clay mineralogy) were larger than those of salt content. The change in hydraulic properties and the dependency on the procedure use may be attributed to aggregate breakdown and plugging of macropores, clay heterogeneity and dispersion, calcium carbonate and gypsum dissolution, all damaging soil structure. Simulation runs illustrated that hydraulic properties obtained with saline water resulted in better model performance than when using the conventional procedures with fresh/tap water. A proper set of initial hydraulic parameter values is an important step in model calibration and thus must be characterized correctly. Therefore, to hydraulically characterize saline soil, procedures should be standardized. Moreover, the approach and results of this study may provide a reference for characterization of soil hydraulic properties of saline and sodic soils effectively.

Soil drainage water and nutrient leaching in winter wheat field lysimeters under different management practices

The objective of this study was to investigate the effects of different irrigation strategies, planting methods and various nitrogen application rates on water and solutes transport under the root zone. Furthermore, Hydrus-1D model was used to simulate the drainage water and soil nutrient (N, K) leaching by the drainage water in the experimental treatments. Results indicated that variable alternate furrow irrigation (VAFI) decreased drainage water by a mean value of 30% compared with that obtained in full irrigation (OFI); whereas, precipitation played a significant role in increasing drainage water. About 60% and 80% of drainage water occurred due to rainfall in the full irrigation strategy in the first and second year, respectively, and most of the drainage water in VAFI strategy occurred due to precipitation as 90% for both years. Additionally, precipitation during the low growth stages of winter wheat (fall and winter) had noticeable influence on solute transport. On average, 75%, and 65% of leached nitrate, and potassium, respectively, occurred due to rainfall in OFI method; however, about 85% of leached nitrate, and potassium in VAFI strategy occurred due to rainfall events. Application of 150 kg N ha−1 decreased about 40% nitrate leaching below the root zone in comparison with that obtained in 300 kg N ha−1 without significant difference in crop yield. Hydrus-1D model reasonably simulated solute transport especially nitrate and potassium. Therefore, Hydrus-1D could be a successful tool for predicting N, and K transport and making proper management decisions to improving environmental problems.

Pore-to-Core Upscaling of Solute Transport Under Steady-State Two-Phase Flow Conditions Using Dynamic Pore Network Modeling Approach

We present a solute transport model, developed by employing a dynamic pore network modeling approach, to investigate dispersive solute transport behaviors in consolidated porous media. The model is capable of upscaling solute transport processes from pore to core, under two-phase fluid configurations. The governing equations of fluid flow, fluid displacement, and solute transport are solved at the pore level. A heavily parallelized computing scheme is utilized to simulate dynamic fluid displacements and transport processes in a core-scale pore network constructed from micro-computed tomography (micro-CT) images of a Berea sandstone sample. A series of solute transport simulations are conducted under the single-phase condition to validate the model by comparing the computed longitudinal dispersion coefficients against the experimental data over a wide range of Peclet numbers (Pe), i.e., $$3 \times 10^{-2}{\sim}3 \times 10^5$$ . The model is then used to simulate solute transport under two-phase fluid configurations to examine the effects of the non-wetting phase saturation on solute transport behaviors. More specifically, solute transport is studied at different Pe and different water saturations obtained at the end of imbibition processes. We find that in a two-phase system, the longitudinal dispersion coefficient substantially increases with enhanced convective mixing under the transport regime with high Pe but decreases with restricted diffusive mixing at low Pe. In addition, the results indicate a non-monotonic dispersion–saturation relation under the convective transport condition. We illustrate that a strong correlation exists between the extent of dispersive mixing and the heterogeneity of fluid saturation profile across a core-scale network.

Stochastic assessment of the effect of land-use change on nonpoint source-driven groundwater quality using an efficient scaling approach

Assessing the impact of a land-use change (LUC) or change in land-use management on nonpoint source-driven groundwater quality in heterogeneous aquifers requires complex analysis. Stochastic methods have been used to account for prediction uncertainty but at high computational cost, which significantly limits the application of these approaches. As an efficient alternative, this study evaluates the application of a meta-analytical solution for evaluating the change in contaminant breakthrough curves at extraction wells in response to LUC. The solution uses the concentration percentiles from a reference stochastic simulation of water flow and solute transport in a groundwater system, assuming a reference land-use distribution pattern. Reference land-use controls the spatially variable rates of both, recharge and contaminant mass loading. The effect of a LUC is evaluated by scaling the ratio between the reference and the new (post-LUC) average input concentrations. The validity of the proposed meta-analysis tool is tested by comparing the results of the meta-analytical solution with those from a full stochastic simulation of the post-LUC scenario. Simulation results show that the accuracy of the meta-analytical solution is best when the regional average recharge rates for both pre- and post-LUC remain approximately unchanged, for any change in contaminant mass loading. Results also indicate that changes in spatial variability and pattern of the recharge rate do not significantly impact the flow field, travel times, and resulting concentrations, if the magnitude of local recharge remains about the same. Lastly, the results show large variability among wells of (and—for an individual well—uncertainty about) the time lag between the time of LUC and the time of consequential effective change in concentrations across wells in the affected region, captured here using statistical metrics.

A distributed-order time fractional derivative model for simulating bimodal sub-diffusion in heterogeneous media

Bimodal transport due to dual-mobile advection or mass exchange between mobile and immobile zones has been widely observed for pollutants moving in heterogeneous media. Existing nonlocal transport models, however, cannot capture either heterogeneity in mobile domains or bi-peak concentration phenomena for breakthrough curves (BTCs) of solute transport in complex media. Therefore, this study proposed a dual heterogeneous domain model (DHDM) framework, conceptualized from the distributed-order time fractional derivative that incorporates a broader spectrum of particle movement in the dual-domain system than the classical single-domain model. A Lagrangian scheme was developed to solve the DHDM using the Bernoulli trial with the transition probability simulating mass exchange between different domains. Phase transition probabilities calculated based on the zeroth spatial moments of dual domain transport equations were used to simulate the particle partitioning between different domains. Applications to transport experiments in silt–clay columns conducted in our laboratory showed that the DHDM yielded significant improvement for simulating solute transport behavior compared to the single-domain anomalous transport model. Parameter analysis of the DHDM further revealed that the velocity disparity of the two domains controlled the BTC rising time and the shape of the second peak in the BTC. When the velocity of the slow domain is extremely small, the second BTC peak develops a ‘shoulder’ characteristic during the late-time tail of the BTC. This study improved our understanding for bimodal transport in a dual domain system and provides a feasible tool for capturing a wide range of bimodal transport in complex media.

Numerical analysis of sensor-based flood-floor ebb-and-flow subirrigation system with saline water

ABSTRACT Ebb-and-flow is a well-known and common greenhouse irrigation system. In this paper, a specific numerical model for a sensor-based flood-floor ebb-and-flow irrigation system has been developed using one-dimensional Richards’ equation for water flow and two-site sorption advection-dispersion model for simulation of solute transport. The model was validated against analytical and numerical results proposed in the literature, and acceptable agreement was obtained in all cases. Maximum RMSE (MAE) was 0.053 m (0.273 m) for water flow and 0.021 (0.065) for solute transport. The model was then used for a parametric study of how different factors (such as salinity of irrigation water, triggering and cutoff thresholds, sensor position and flooding depth) would affect transpiration of sensor-equipped solitary pot under ebb-and-flow irrigation system for two different crops. It was revealed that sensor position and water salinity were significantly influential on transpiration of the pots. Finally, the synchronized sensor equipped-sensorless module of the numerical model was used to optimize the flood-floor ebb-and-flow irrigation system considering different factors using Taguchi approach. It was concluded that installing a water content sensor, whose triggering and cutoff thresholds were at −10 and −2 m, at the upper-half of the center-positioned pot would result in higher transpiration and lower unevenness of water distribution.