Hydraulic Parameters Research Articles

The allocation and erosive effects of on-wall flow on loess gully heads

On-wall flow and jet flow converted from upstream flow at the gully brink play important but differential roles in gully head erosion. However, on-wall flow allocation at gully brinks and its erosive effects on vegetation-covered loess gully heads remain unclear. Simulated flow scouring experiments were conducted on grass-covered gully heads under different inflow rates (3.0–7.2 m3/h) and upstream slope gradients (1–7°) to investigate the on-wall flow proportion and relationships with the hydraulic parameters of upstream flow, as well as erosion and sediment yield characteristics. The results showed that the on-wall flow rates linearly increased with increasing inflow rate (P < 0.01) but decreased as the slope gradient increased, whereas the jet flow rate increased as the inflow rate or slope gradient increased. The on-wall flow proportions ranged from 24.6 % to 58.7 % and decreased with increasing inflow rate or slope gradient, whereas the jet flow proportions exhibited the reverse trend. The on-wall and jet flow proportions had the most significant correlations with the Froude number, following a negative power function and positive logarithmic function (P < 0.01), respectively. The breadth and depth of the scour hole that developed on the gully headwall increased linearly with the on-wall flow rate (P < 0.01). Jet flow had a limited effect of accumulation on scour hole development, with average increases of 13.5 % and 8.6 % in the breadth and depth, respectively. On-wall flow generally contributed 17.6 %–48.4 % of the sediment driven by on-wall and jet flows but played a dominant role under 1° upstream slope gradient and 3.6 m3/h inflow rate in which the contribution reached 64.9 %. On-wall flow might have dominated the runoff volume and played a dominant role in the sediment yield under low upstream slope gradients and inflow rates and was essentially responsible for scour hole development, which could destabilize and retreat the gully head.

Exploring the Effects of Fissures on Hydraulic Parameters in Subsurface Flows from the Perspective of Energy Changes

Reynolds number (Re), pore water pressure (P), and water flow shear force (τ) are primary indicators reflecting the characteristics of subsurface flow. Exploring the calculation of these parameters will facilitate the understanding of the hydrodynamic characteristics in different subsurface flows and quantify their differences. Hence, we conducted a study to monitor soil water content, matrix potential, and pore water pressure in two typical soil profiles (with and without fissures). The distribution of Re, P, and τ in both matrix flow (MF) and preferential flow (PF) were calculated with an improved calculation method, focusing on their energy changes. Results showed that these hydrologic parameters are quite different between MF and PF. Re values in MF remained below 0.1, indicating lower water flow velocities, while the Re values ranged from 0.8 to 2 in PF, indicating higher flow velocities. The P values in PF was tens to hundreds of times higher than that in MF, which is mainly due to the rapid accumulation and leakage of water within soil fissures. Additionally, the larger hydraulic radius and gradient in PF also resulted in higher τ values in PF (2~6 N m−2) than in MF (0~1.5 N m−2). In PF, the pressure potential was the significant factor for τ, while τ in MF was dominated by the matrix potential and varies with the magnitude of the matrix potential gradient. This study suggests that Re, P, and τ could be considered as the major indexes to reflect dynamic characteristics of subsurface flow.

Improved partitioning between matrix and macropore flow: Novel bimodal lognormal functions for water retention and hydraulic conductivity in pumice and non-pumice soils

Dual-porosity models have been shown to improve models of soil–water movement and enhance the water balance of structured soils. In this study we introduce novel, continuous, closed-form, bimodal, lognormal functions for soil–water retention, θ(ψ), and unsaturated hydraulic conductivity, K(ψ), enhancing traditional models and significantly improving predictions, particularly for pumice soils. These functions incorporate the thresholds for (a) water pressure, (b) soil water content, and (c) unsaturated conductivity, which accurately differentiates macropores from matrix pores. Validation using 313 observation points from laboratory data shows an increase in the Nash–Sutcliffe efficiency coefficient of θ(ψ) from 0.94 to 0.97, and for K(ψ) from 0.83 to 0.95. The model also requires six constant, semi-empirical tortuosity parameters and provides a physically constrained approach for the hydraulic parameters that reduces non-uniqueness risks. The derived functions yield improved predictions and enable the computation of macropore soil water content and flow contributions, with potential applications for the advancement of preferential flow modelling.

River quality management: Integrating uncertainty, failure probability, and assimilation capacity

Managing river water quality is challenging due to uncertainties in hydraulic and hydrologic parameters. This study integrates the symmetric exponential function (SEF) approach for solving the advection-dispersion equation with the Monte Carlo method in MATLAB. This combination allows us to explore the river's assimilation capacity and the failure probability (Pf) of maintaining desired water quality standards. Here, Pf represents the likelihood of pollutant concentrations exceeding acceptable limits under varying river conditions. A key contribution of this study is the introduction of a novel equation, developed using the Genetic Programming (GP) soft computing tool, to calculate assimilation capacity considering the failure probability of water quality provision. This equation provides a valuable tool for risk assessment in water resource management by quantifying pollutant assimilation dynamics. Its robustness is validated through high Coefficient of Determination (R2) and Overall Index (OI) values near 1, along with low Root Mean Square Error (RMSE) and Mean Absolute Error (MAE). The study identifies critical river characteristics, such as flow velocity and pollutant load, significantly influencing the reliability index (β). By outlining how adjustments in these parameters can achieve a target reliability index (β=4.526), our study offers a practical approach to safeguarding river ecosystems. For example, increasing flow velocity by 76 % can shift the river from a safe state (Pf=3×10−5) to a hazardous state (Pf=1), while a 44 % decrease in velocity allows for 57 % more pollutant assimilation. These findings highlight the importance of flow control as a cost-effective strategy for mitigating high pollutant concentrations and ensuring sustainable water quality management. By integrating numerical approaches with reliability sampling methods and soft computing techniques, this study enhances understanding of river system dynamics and supports informed decision-making for protecting water resources.

Integrating calibrated PTFs and modified OpenKarHydro framework to map the responses of ecohydrological processes to climate change across the Loess Plateau

Climate change is exacerbating the risk of soil water stress, soil erosion and ecological degradation in the Loess Plateau (CLP). However, the spatiotemporal dynamics of ecohydrological processes, including soil water balance (SWB), soil erosion (SE) and vegetation net primary productivity (NPP), in response to climate change remain isolated and ambiguous. Additionally, existing studies often rely on limited measured data from specific watersheds with particular land use/cover types to establish the pedotransfer functions (PTFs) for single soil hydraulic parameter (SHP). There is a notable lack of comprehensive studies that explore PTFs encompassing all SHPs applicable to the entire CLP, predict the spatiotemporal response of ecohydrological processes to climate change, and identify future risk periods and regions for SWB, SE, and NPP. To address these gaps, we first screened the PTFs applicable for the whole CLP, and then projected the spatiotemporal dynamics of ecohydrological processes from 2020 to 2030 under CMIP6 scenarios by integrating the PTFs with the OpenKarHydro, RUSLE, and CASA model, and finally classified the soil water content, SE and NPP to identify their risk periods and regions. The results showed that the PTFs achieved satisfactory accuracy, with RMSEs for Ks, θs, θr, α, n, θfc and θw being 2.682, 0.109, 0.016, 0.111, 0.897, 0.060 and 0.058, BIASs of 1.365, 0.182, 0.012, 0.031, 0.214, 0.057 and 0.048, R2s of 0.445, 0.500, 0.430, 0.400, 0.694, 0.453 and 0.453, and NSEs of 0.645, 0.737, 0.874, 0.349, 0.567, 0.756 and 0.458, respectively. SE decreased significantly from 2020 to 2030, with an average annual rate of −6.18 %. Soil water storage decreased significantly between 2020 and 2030, and declining from southeast to northwest. The proportion of area in the wet zone decreased by 4 %, while the proportion of area in the dry zone increased by 11.9 % from 2020 to 2030. The average NPP in 2020–2030 is 320.07 gC·m−2·a−1, with the largest in summer and smallest in winter, and decreasing from southeast to northwest. The NPP increased significantly in 2020–2030, with average annual value and rate of 19.61 gC·m−2·a−1 and 70.84 %, respectively. The NPP increased significantly in spring, summer and autumn, but remained stable in winter. This investigation bridges the gap between existing soil properties, missing SHPs, the strong regional applicability of PTFs, and isolated ecohydrological processes. It is promises to provide valuable insights into the response to climate change in the CLP and another water-limited regions.

Hydraulic and Hydrogeochemical Characterization of Carbonate Aquifers in Arid Regions: A Case from the Western Desert, Egypt

Using geochemical and pumping test data from 80 groundwater wells, the chemical, hydrologic, and hydraulic properties of the fractured Eocene carbonate aquifer located west of the Al-Minya district, the Western Desert, Egypt, have been characterized and determined to guarantee sustainable management of groundwater resources under large-scale desert reclamation projects. The hydrochemical data show that groundwater from the fractured Eocene carbonate aquifer has a high concentration of Na+ and Cl− and varies in salinity from 2176 to 2912 mg/L (brackish water). Water–rock interaction and ion exchange processes are the most dominant processes controlling groundwater composition. The carbonate aquifer exists under confined to semi-confined conditions, and the depth to groundwater increases eastward. From the potentiometric head data, deep-seated faults are the suggested pathways for gas-rich water ascending from the deep Nubian aquifer system into the overlying shallow carbonate aquifer. This mechanism enhances the dissolution and karstification of carbonate rocks, especially in the vicinity of faulted sites, and is supported by the significant loss of mud circulation during well drilling operations. The average estimated hydraulic parameters, based on the analysis of step-drawdown, long-duration pumping and recovery tests, indicate that the Eocene carbonate aquifer has a wide range of transmissivity (T) that is between 336.39 and 389,309.28 m2/d (average: 18,405.21 m2/d), hydraulic conductivity (K) between 1.31 and 1420.84 m/d (average: 70.29 m/d), and specific capacity (Sc) between 44.4 and 17,376.24 m2/d (average: 45.24 m2/d). On the other hand, the performance characteristics of drilled wells show that well efficiency ranges between 0.47 and 97.08%, and well losses range between 2.92 and 99.53%. In addition to variations in carbonate aquifer thickness and clay/shale content, the existence of strong karstification features, i.e., fissures, fractures or caverns, and solution cavities, in the Eocene carbonate aquifer are responsible for variability in the K and T values. The observed high well losses might be related to turbulent flow within and adjacent to the wells drilled in conductive fracture zones. The current approach can be further used to enhance local aquifer models and improve strategies for identifying the most productive zones in similar aquifer systems.

Conceptual Modflow Model Application with Boundary Condition Analysis

In this study, groundwater flow modeling of the aquifer in the Porsuk Basin, located in the Eskisehir province of Turkey, was conducted using the MODFLOW modeling method and the Groundwater Modeling System (GMS) software interface. The modeling was performed using the finite difference method under the assumption of steady flow, incorporating well data, aquifer boundaries, topographic elevations, riverbed data, and hydraulic parameters from the 1971 water year. This study provided insights into changes in groundwater over time and yielded a comprehensive water budget. By creating a three-dimensional solid model of the Porsuk Basin aquifer, general information about the region's hydrogeology was obtained. The results are expected to play a significant role in identifying potential drilling areas for groundwater extraction. Since the modeling utilizes objects from geographic information systems, it will enhance visualization of the regional structure for researchers, allowing them to observe the topographic features of the area in three dimensions. Additionally, it is anticipated that the groundwater conceptual model will inform drilling studies and serve as a foundation for research on groundwater transport and pollution.

Numerical and sensitivity analysis of hydraulic characteristics of triangular labyrinth side weir

Triangular labyrinth side weirs have significantly more discharge capacity than traditional nonlinear weirs, and the complex hydraulic parameter interaction mechanism has been the research focus. This study used Computational Fluid Dynamics (CFD) to analyze the side weir's flow characteristics. Then, the Bayesian optimization algorithm and Extreme Learning Machine (BELM) developed a prediction model for the side weir's discharge coefficient. Finally, Sobol's method performed a sensitivity analysis for hydraulic parameters. The results show that the main channel's streamline is evenly distributed and begins to shift when it is close to the side weir. The overflow front is increasing and secondary flow also increases. BELM's Mean Absolute Percentage Error and Root Mean Square Error are 8.793 % and 0.455 in the testing stage, respectively, declined by about 56.24 % and 32.29 % compared with ELM; Froude number Fr, weir crest angle θ and the ratio of overflow front length to weir head l/h1 are the most critical hydraulic parameters affecting the discharge coefficient, the global sensitivity coefficients are 0.4393, 0.4218 and 0.4152, respectively.

Fitting the van Genuchten model to the measured hydraulic parameters in soils of different genesis and texture at the regional scale

Fitting the van Genuchten model to the measured hydraulic parameters in soils of different genesis and texture at the regional scale

Soil hydraulic property maps for the contiguous United States at 100‐m resolution and seven depths: Code design and preliminary results

AbstractEstimates of the van Genuchten (1980, abbreviated as VG) parameters and saturated hydraulic conductivity (Ks) were made for the contiguous United States at a resolution of 100 m and seven soil depths by combining the SoilGrids+ (SG+) soil property maps of Ramcharan et al. with the R3H3 member of the Rosetta3 hierarchical pedotransfer functions (PTFs) of Zhang et al. To this end, we developed multi‐threaded code that significantly speeds up computation (up to a factor 25) depending on the level of parallelism. We verified estimates first by calculating simple summary statistics of estimated basic properties of SG+ with actual measured soil properties for 14,113 pedons in the National Cooperative Soil Survey (NCSS) (2023) labsample database. Next, we computed summary statistics of PTF‐estimated moisture contents for NCSS and SG+ data. The results show estimation errors are dominated by intrinsic errors of the PTF, and that (potentially correctable) systematic errors are present in SG+ soil properties and PTF estimates. The resulting hydraulic property maps contain well over 750 million points for each of the seven layers and show considerable horizontal and depth variation for each VG parameter and Ks, except the VG “n” parameter, which is dominated by values between 1.25 and 1.6. The hydraulic property maps are 99.9% complete, and we demonstrate that plausible profiles and uncertainty information can be generated for virtually each point. The maps are available as two multi‐channel GeoTIFF maps per SG+ layer: one with the five hydraulic parameters and one with the corresponding covariances.

Увеличение продолжительности работы установки электроцентробежного насоса для повышения безопасности и сокращения рисков при спускоподъемных операциях

The study analyzes basic complications associated with generating mechanical impurities; numeric modeling to optimize the electrical submersible pump operation has been performed. Modeling the fluid flow in combination with erosion models has been used to determine the key factors facilitating the erosion of the electrical submersible pump action wheel impeller during oil production. In order to obtain a complete view of the flow and erosion properties, modeling has been performed for the entire fluid flow route. Modeling of the flow and the action wheel component erosion has been performed to determine the impact of the flow imbalance as well as to ensure the boundary conditions at the input for the modeling. The modeling results have shown: the rotational flow introduced in the mixing point at the action wheel input significantly impacts the strength properties of the modeled object; the area most exposed to erosion is the area of the action wheel impeller. The values of loss of material have been determined via the time sequence of geometry changes in the compounds in combination with the actual hydraulic operational parameters. The flow trajectories proposed for modeling show that the recirculating flow is the dominant. It is widely adopted that recirculation is one of the main causes of increased erosion in the compounds of the electrical submersible pump as it leads to the high frequency of high-speed impacts by sand particles. In order to reduce the erosion risk, the analysis aiming to identify the impeller geometrical parameters has revealed that the erosion rate in the compound can be reduced by decreasing the inclination of the action wheel impeller. The study has demonstrated that numerical modeling can provide a detailed and precise evaluation of erosion potential in cases associated with the high velocity of fluid flow with particles of mechanical impurities.

Statistical characteristics of aquitard hydraulic conductivity, specific storage and porosity

Aquitard hydraulic conductivity (K), specific storage (Ss) and porosity (ϕ) are critical to the modeling of flow and mass and heat transport in groundwater systems, yet comprehensive insights into their statistical characteristics remain limited. This study compiles 456 K values from 47 studies, 202 Ss values from 49 studies and 239 ϕ values from 18 studies, followed by rigorous statistical analysis of aquitard K, Ss and ϕ with dependence on the features of test method, lithology and burial depth. Additionally, the interrelationships among K, Ss and ϕ are analyzed. A predictive model for aquitard K and Ss is developed utilizing the random forest (RF) approach. The statistical findings unearth that laboratory experiments tend to yield lower K while higher Ss, particularly when burial depths are less than 50 m. Clay, silt clay, and glacial till exhibit comparable K values, in contrast to significantly lower K values observed in shale formations. Aquitard K, Ss, and ϕ exhibit depth-decaying characteristics. Between two adjacent depth intervals, 0 ∼ 50 m, 50 ∼ 150 m and >150 m, there is nearly an order of magnitude difference in either K or Ss. Aquitard K tends to increase with ϕ when ϕ is below 0.25 and decreases vice versa. This phenomenon is probably attributed to the significant amount of immobile-bound water adhering to clayey particles within the aquitard with a high ϕ value, which reduces the pore space available for water flow. Elastic Ss, which indicates recoverable aquitard groundwater depletion, is about one order of magnitude lower than inelastic Ss, which indicates permanent aquitard groundwater depletion, implying that the majority of groundwater depletion from aquitards under long-term groundwater exploitation is non-recoverable. The RF-based predictive model indicates that ϕ is a highly important feature for predicting K and outperforms lithology, test method and burial depth in predicting Ss. The statistical findings of this study present a valuable basis for aquitard hydraulic parameters.

Vineyard cover crop management strategies and their effect on soil properties across Europe

AbstractVineyard soils are often of inherently poor quality with low organic carbon content. Management can improve soil properties and thus soil fertility. In European wine‐growing regions, a broad range of inter‐row management strategies evolved based on specific local site conditions and the varying effects of management intensities on soil, water balance, yield and grape quality. Accordingly, there is a need to investigate the effects of locally common cover crop management strategies and tillage intensity on soil organic carbon content and soil physical parameters. In this study, we investigated the impact of the most common inter‐row management practices in Austria, France, Romania and Spain. In all countries, we compared paired sites. Each site with cover crops and inter‐row management of low intensity was compared with one site with (temporarily) bare soil and high management intensity. All studied sites with cover crops and low management intensity, except those in Spain, had higher organic carbon contents than the paired more intensively managed vineyards. However, the highly water‐limited Spanish vineyards with temporary cover crops had lower organic carbon contents than the paired sites with bare soil. Sites with more organic carbon had better results for bulk density, percolation stability (PS), hydraulic conductivity and available soil water, with soil hydraulic parameters being less pronounced than others. Country comparison of inter‐row weed control systems showed that PS was particularly low in sampled vineyards in Romania and Spain, where weed control is based on intensive mechanical tillage. Alternating management systems with tillage every second inter‐row showed a decrease in soil structure compared with permanent green cover. Thus, inter‐row management with cover crops and reduced tillage increases soil organic carbon content and improves soil structure compared with bare soil management. If local constraints, such as water scarcity, do not allow year‐round planting, alternating inter‐row management with several years of alternating periods may be an option to mitigate those adverse effects. However, negative impact on the soil structure occurs with the very first tillage operation, whereas negative effects on the carbon balance only appear after long‐term use of tillage.

Comparison of novel physics-guided machine learning models with empirical equations for predicting longitudinal dispersion coefficient in diverse natural river systems

Accurate calculation of the longitudinal dispersion coefficient (LDC) is crucial for assessing river fate and transport processes. Traditional methods, including experimental, analytical, and empirical equations, each have their challenges, such as complexity and susceptibility to errors. This study introduces and evaluates advanced physics-guided machine learning (ML) models, including hybrid hydrodynamic-particle simulation (HHPS), physics-enhanced ML (PEML), and physics-regularized regression trees (PRRT). Our dataset includes hydraulic and geometric parameters collected from 50 diverse rivers across the US and the UK, encompassing variables such as river width (W), flow depth (H), average flow velocity (U), flow shear velocity (U*), concentration of the solute (μgL) and LDC (ε). Input parameters were dimensioned using the UU∗ and WHratios, while the dimensionless LDC (εHU∗) served as the model output. By embedding core hydrodynamic principles such as continuity, mass and momentum conservation, and Navier-Stokes equations into the ML algorithms, these models achieved high predictive accuracy. Specifically, the HHPS model reached a coefficient of determination of 0.997, and both PEML and PRRT demonstrated Nash-Sutcliffe efficiencies exceeding 0.950, significantly outperforming traditional empirical equations. Additionally, SHapley Additive exPlanations (SHAP) sensitivity analysis revealed the substantial influence of channel width and flow conditions on LDC predictions. Our study highlights the superior performance of physics-guided ML models in providing accurate and reliable tools for river water quality management, demonstrating their substantial improvements over traditional methods and their potential broader applications in environmental studies.

Inverse analysis of soil hydraulic parameters of layered soil profiles using physics‐informed neural networks with unsaturated water flow models

AbstractInformation about the spatial distribution of soil hydraulic parameters is necessary for the accurate prediction of soil water flow and the coupled movement of chemicals and heat at the field scale using a process‐based model. Physics‐informed neural networks (PINNs), which can provide physical constraints in deep learning to obtain a mesh‐free solution, can be used to inversely estimate soil hydraulic parameters from less and noisy training data. Previous studies using PINNs have successfully estimated soil hydraulic parameters for homogeneous soil but estimating such parameters of layered soil profiles where the interface depth and the parameters are unknown still has some difficulties. The objective of this study was to develop PINNs to inversely estimate the distribution of soil hydraulic parameters, such as saturated hydraulic conductivity and α and n of the Mualem–van Genuchten model directly within layered soil profiles by predicting changes in the pressure head from training data based on simulation results at given depths during infiltration. The impact of factors affecting PINNs performance, such as the weights assigned to each component of the loss function, time range used in error computations, and number of samples used to assess the physical constraint, was investigated. By assigning a larger weight to the physical constraint and excluding the earlier stage of infiltration in the loss function, the changes in the pressure head and the three soil hydraulic parameter distributions within the layered soil were successfully estimated. The developed PINNs can be further applied to more complex soils and can be improved.

Numerical Modeling of Hydrological Mechanisms and Instability for Multi-Layered Slopes

The process of rainwater infiltration into unsaturated multi-layered slopes is complex, making it extremely difficult to accurately predict slope behaviors. The hydrological mechanisms in multi-layered slopes could be significantly influenced by the varying hydraulic characteristics of different soils, thus influencing slope stability. A numerical model based on Hydrus 2D was constructed to investigate the hydrological mechanisms of multi-layered slopes under different slope inclinations and rainfall intensities. The results revealed hydraulic processes in response to rainfall in unsaturated multi-layered slopes, in which layered soils retard the advance of wetting fronts and affect seepage paths in the slope. The results also showed the characteristics of hydraulic parameters, including pore water pressure and moisture content, under different conditions, and explained the crucial factors at play in maintaining slope stability.

Application of surface geophysical investigations and pumping test data analysis for better characterization of aquifer hydraulic parameters, Upper Awash Sub-Basin, Central Ethiopia

Study regionUpper Awash River Sub-basin in central Ethiopia Study focusThe research aimed to estimate aquifer parameters and conduct resistivity modeling in the study area. Data were collected from 890 Schlumberger Vertical Electrical Soundings (VES) points strategically distributed across six selected zones. The VES data were analyzed using computer modeling and curve-matching techniques to interpret the resistivity responses of the geological formations. This analysis provided insights into the subsurface characteristics and identified potential aquifer zones. New insights for the regionThe interpreted curves from the VES data and Root Mean Square (RMS%) values ranging from 0.36 % to 1.67 %, indicating a good fit between observed and modeled data. Resistivity modeling and geoelectrical sections were produced along specific lines, visualizing subsurface structures and confirming the resistivity responses of the geological formations. To determine aquifer parameters, the Dar-Zarrouk parameters were calculated from the interpreted curves, yielding crucial information about hydraulic properties such as hydraulic conductivity (K) and transmissivity (T). Hydraulic conductivity values ranged from 1.1 m/day to 21.7 m/day, while transmissivity values ranged from 144.3 m²/day to 2763.3 m²/day. Validation against borehole pumping test data showed strong correlations: R² = 0.9873 for transmissivity and transverse resistance, and R² = 0.9352 for hydraulic conductivity and resistivity. These findings enhance the understanding of the complex volcanic aquifer characteristics in the Upper Awash River Sub-basin, aiding sustainable groundwater management and development in the area.

Experimental study on flow and heat transfer of High-Pressure helium flow in compact helical tube heat exchanger in HTGR

In this study, experiments were conducted on helium flow within a compact helical tube heat exchanger. High-temperature helium and deionized cooling water, flowing in opposite directions, were allocated to the shell side and tube side, respectively. The Reynolds number of helium ranged from 3 × 103 to 1.6 × 104 under an operational pressure of 2.1 MPa. Measurements were taken for helium temperature, tube outer-surface temperature distribution, pressure, and mass flow rate to illustrate the flow and heat transfer characteristics. Sensitivity analysis of thermal–hydraulic parameters, including heat transfer power and efficiency under various operational conditions, was performed. The convective heat transfer coefficient and pressure drop were calculated. The transition point (Re = 9000) between laminar and turbulent flow of helium on the shell side was identified, and empirical correlations for the Nusselt number and friction coefficient were proposed. The experimental results indicate that the heat transfer performance of helium in the compact helical tube heat exchanger exceeds that of water in large-scale helical tube heat exchanger by approximately 20.75 % under fully developed turbulence, while the flow resistance characteristics remain essentially consistent.

Modeling soil water dynamics of an intensively cultivated histosol

AbstractUnderstanding soil–water dynamics in cultivated organic soils (histosols) is crucial for sustainable agriculture, ecosystem preservation, and climate change mitigation. Data on the soil water retention curves (SWRCs) of these histosols in Canada are not readily available in the literature. The Hydrus‐1D model was used to predict SWRCs for a cultivated organic soil in Quebec. The model was validated with matric potential measured in a soil column at 10, 26, and 48 cm. The optimized hydraulic parameters in the Hydrus model resulted in field capacity moisture contents of 0.6, 0.696, and 0.49 cm3 cm−3 at the three depths, respectively. Wilting point moisture contents of 0.165, 0.107, and 0.14 cm3 cm−3 were obtained at the same depths. Model accuracy was confirmed with a Nash–Sutcliffe efficiency of 0.76, 0.91, and 0.78 for the three depths, respectively. These findings can inform decisions regarding water management practices such as irrigation and drainage requirements for intensively cultivated organic soils.

Dynamic and empirical methods for field capacity estimation in fine textured soils with a coarse interlayer

Field capacity (FC) is an important soil hydraulic concept in soil science and irrigation management. It is generally determined from soil water content in a soil layer when soil profile reaches a steady pressure head or negligible drainage flux from an initially saturated soil. However, the proposed criteria are mainly tested for uniform soils and vary with soil textures. To quantify FC in layered soils, a Richards equation-based model was used to describe water flow in fine-textured soils with a coarse interlayer. With calibrated soil hydraulic parameters for loam and sand from infiltration measurements, drainage from saturation was simulated in the loam with a sand interlayer. A relative drainage rate (δ) was defined as a function of water storage and drainage flux to analyze soil water status at FC. Soil water content in the upper loam layer of layered profiles was improved compared with that in the uniform loam, which was negatively correlated with buried depth but positively correlated with thickness of the sand layer for a specified δ. Under different buried depths and thicknesses, soil water content decreased with the decline of δ and decreased rapidly as δ reduced to 1 % d−1. The drainage flux at δ = 1 % d−1 changed within a range of 0.056–0.26 cm d−1, and soil water content reached to 0.278–0.346 cm3 cm−3, which accounted for 70–87 % of the saturated water content of loam. Although the FC in the upper fine-textured soil layer varied for different buried depths and thicknesses of coarse interlayer, the proposed dynamic method is reliable and universal to estimate the FC in the above layered soils at δ = 1 % d−1. An empirical equation was also developed to calculate the FC in fine-textured soils with different buried depths and thicknesses of a coarse interlayer based on the critical δ value.