Flow Path Research Articles

Tip air injection and plasma DBD actuator as active flow control in axial compressors - Review

Abstract Nowadays, flow control is extensively studied in fluid dynamic domain, and the principal aim is to manipulate the flow path that passes over a surface to reduce losses generated by turbulent flow, generated by detachment. The Flow control involves passive or active methods, and for approximately thirty years, active flow control methods have been studied and developed in axial compressors more intensively, active flow control methods are those in which energy is added to the fluid through an external source. The flow control reduces losses generated by secondary flow effects, such as flow detachment at rotor tip, flow detachment at casing and flow detachment at stator vanes hub. The present paper provides a summary of investigation of active flow control methods tip air injection and plasma actuators.

Density spatial distribution and anisotropy of full-scale bentonite-sand blocks

Constructing buffer barriers using high-quality bentonite-based blocks is crucial for the secure disposal of high-level radioactive waste (HLW). However, the density distribution and anisotropy of full-size bentonite-sand blocks are not yet fully understood, which may affect their quality and performance as buffer barriers. This study aimed to address this critical issue by conducting computed tomography (CT) scanning tests on full-size bentonite-sand blocks in the radial, angular, and vertical directions to obtain CT numbers and density. The analysis of the density distribution and anisotropy in these three directions followed. Subsequently, thermal conductivity and hydraulic conductivity tests were conducted in these directions to explore their relationship with density anisotropy. The findings revealed that the block density gradually decreased with increasing radius in the radial direction and first decreased slightly and then increased slightly in the angular and vertical directions. However, the density of the bentonite-sand mixture inside the blocks ranged from 1.96 g/cm³ to 2.09 g/cm³, with minimal differences. Additionally, there was no significant difference in density between the radial and angular directions, with an average density of 2.06 g/cm³. The density in the vertical direction was slightly lower, with an average density of 2.04 g/cm³. Scanning electron microscope images validated the presence of anisotropy in block density caused by horizontally oriented bentonite aggregates resulting from axial compression. However, the presence of quartz sand and incompressible impurities in the blocks created preferential flow paths, mitigating the impact of block density anisotropy on thermal and hydraulic conductivity. Consequently, thermal conductivity, averaging 1.60 W/m·K, and hydraulic conductivity, averaging 3.12 × 10−10 cm/s, demonstrated isotropy in all three directions. These block parameters meet the design requirements for an effective buffer barrier. This study systematically reveals the density distribution of full-size bentonite-sand blocks and its impact on the performance, providing important references for the safe design of HLW repositories.

Conceptualizing Controlling Factors for PFAS Salting Out in Groundwater Discharge Zones Along Sandy Beaches.

Understanding fate and transport processes for per- and poly-fluoroalkyl substances (PFAS) is critical for managing impacted sites. "PFAS Salting Out" in groundwater, defined herein, is an understudied process where PFAS in fresh groundwater mixes with saline groundwater near marine shorelines, which increases sorption of PFAS to aquifer solids. While sorption reduces PFAS mass discharge to marine surface water, the fraction that sorbs to beach sediments may be mobilized under future salinity changes. The objective of this study was to conceptually explore the potential for PFAS Salting Out in sandy beach environments and to perform a preliminary broad-scale characterization of sandy shoreline areas in the continental U.S. While no site-specific PFAS data were collected, our conceptual approach involved developing a multivariate regression model that assessed how tidal amplitude and freshwater submarine groundwater discharge affect the mixing of fresh and saline groundwater in sandy coastal aquifers. We then applied this model to 143 U.S. shoreline areas with sandy beaches (21% of total beaches in the USA), indirectly mapping potential salinity increases in shallow freshwater PFAS plumes as low (<10 ppt), medium (10-20 ppt), or high (>20 ppt) along groundwater flow paths before reaching the ocean. Higher potential salinity increases were observed in West Coast bays and the North Atlantic coastline, due to the combination of moderate to large tides and large fresh groundwater discharge rates, while lower increases occurred along the Gulf of Mexico and the southern Florida Atlantic coast. The salinity increases were used to estimate potential perfluorooctane sulfonic acid (PFOS) sorption in groundwater due to salting out processes. Low-category shorelines may see a 1- to 2.5-fold increase in sorption of PFOS, medium-category a 2.0- to 6.4-fold increase, and high-category a 3.8- to 25-fold increase in PFOS sorption. The analysis presented provides a first critical step in developing a large-scale approach to classify the PFAS Salting Out potential along shorelines and the limitations of the approach adopted highlights important areas for further research.

Comprehensive analysis of carbon emission reduction technologies (CRTs) in China's coal-fired power sector: A bottom-up approach

Carbon dioxide (CO2) reduction technologies (CRTs) in the coal-fired power sector play an imperative role in the mitigation of environmental challenges and reducing CO2 emissions to help achieve the 2 °C target. However, a compelling necessity persists for a unified framework that can effectively and accurately estimate the costs and potentials associated with CRTs, arising from the diversity of technologies and unit types. Therefore, this study employs a bottom-up approach to analyze the costs and potential of 25 advanced CRTs in the coal-fired power sector, excluding CO2 capture, utilization and storage (CCUS), across 19 types of coal-fired power units. This analysis amalgamates the CO2 reduction supply curve (CRSC) approach with levelized cost of CO2 reduction (LCOC) and a broken-even analysis which could reflect the sector's average carbon reduction cost. The outcomes reveal the following key insights: (1) These technologies collectively harbor a substantial CO2 conservation potential amounting to 925.61 Mt CO2, with a broken-even price of CNY 410.1/tCO2. This reduction could potentially lower 20%–25 % of CO2 emissions in China's power system in 2022, highlighting the crucial role of CRTs in achieving a low-carbon transition and national climate targets; (2)Steam turbine flow path retrofit (T8) not only offers a relatively high CO2 reduction potential (>60 Mt CO2), but also features a low cost (<CNY 150/tCO2) and high profit (>CNY 0.85/MWh). Prioritizing the development of this technology can significantly accelerate the low-carbon transition of coal power; (3) Carbon price have a paramount influence on the cost-effectiveness of CRTs and driving their adoption.

Tailings particle size effects on pollution and ecological remediation: A case study of an iron tailings reservoir

The particle size distribution in tailings notably influences their physical properties and behavior. Despite this, our understanding of how the distribution of tailings particle sizes impacts in situ pollution and ecological remediation in in-situ environment remains limited. In this study, an iron tailings reservoir was sampled along a particle flow path to compare the pollution characteristic and microbial communities across regions with different particle sizes. The results revealed a gradual reduction in tailings particle size along the flow direction. The predominant mineral composition shifts from minerals such as albite and quartz to layered minerals. Total nitrogen, total organic carbon, and total metal concentrations increased, whereas the acid-generating potential decreased. The region with the finest tailings particle size exhibited the highest microbial diversity, featuring metal-resistant microorganisms such as KD4–96, Micrococcaceae, and Acidimicrobiia. Significant discrepancies were observed in tailings pollution and ecological risks across different particle sizes. Consequently, it is necessary to assess tailings reservoirs pollution in the early stages of remediation before determining appropriate remediation methods. These findings underscore that tailings particle distribution is a critical factor in shaping geochemical characteristics. The responsive nature of the microbial community further validated these outcomes and offered novel insights into the ecological remediation of tailings.

Rice response to fluridone following topsoil removal on a precision‐leveled field

AbstractFurrow‐irrigated rice (Oryza sativa L.) production requires additional management options without the cultural strategy of a flood. In 2023, fluridone (Brake) was registered for use in rice production. Precision‐leveling land is an important aspect of maintaining an effective irrigation flow path; however, the removal of topsoil can severely impact crop response from residual herbicides. The objective of this study was to evaluate the impact of fluridone on rice response when applied to a precision‐leveled field following topsoil removal. An on‐farm field study was conducted in 2023 near Osceola, AR, with a Sharkey‐Steele clay complex soil texture. The study consisted of six herbicide treatments applied at the rice three‐leaf stage. These treatments included clomazone, quinclorac, and fluridone at three rates (0.5×, 1×, and 2× of a label rate) and a nontreated control. Across evaluation times, maximum visual rice injury was greater than 25% and 65% for the fluridone 1× and 2× treatments, respectively. Rice canopy coverage was reduced in the fluridone 1× and 2× treatments compared to all other treatments from 6 to 10 weeks after application (WAA). Rice heading in the fluridone 2x treatment was delayed by more than 60 and 30 percentage points at 11 and 12 WAA, respectively, compared to all other treatments. Rough rice yield in the fluridone 2× treatment was reduced by 21% compared to all other treatments. Overall, fluridone 1× and 2× rates caused substantial rice injury and would not be recommended on a precision‐leveled field.

Effects of porosity and porosity distribution in gas diffusion layer on the performances of proton exchange membrane fuel cell

The transmission capability of the gas diffusion layer directly impacts electrochemical reaction and water drainage, consequently, the comprehensive performance and lifetime of the proton exchange membrane fuel cell. The gas diffusion layer's porosity and distribution must be designed effectively as a porous component. This study focuses on experimentally validated models to investigate the effects of porosity and distribution in the gas diffusion layer using three-dimensional numerical simulation. A porosity of 0.6, with a unitary distribution, provides comprehensive improvements in power density while minimizing pressure drops. Additionally, four different distribution patterns of porosity (28 cases) are studied while maintaining an overall porosity of 0.6. The linear porosity distribution (along the flow path) with positive slopes outperforms the alternant distribution due to the cumulative effects on the micro-subsection of the gas diffusion layer. The stepped and/or the sinusoidal distribution can also improve the performances, but just when the step gradient and/or the amplitude are sufficiently small. The adverse effects on the uniformity of current density cause the sinusoidal distribution to be inferior to linear and stepped distribution patterns. The transmission of oxygen significantly affects the dynamic performances, with the distribution patterns of porosity critically influencing sensitivity.

Estimation of Effective Fracture Aperture in Glacial Tills by Analysis of Dye Tracer Penetration.

This study advances a methodology to estimate effective apertures of fractures in glacial tills based on dye tracer infiltration tests and numerical simulations. The approach uses the visible penetration depth of the dye tracer along fracture flow paths as primary information to calculate effective fracture apertures. Further data used in the calculation are the dye tracer input concentration and retardation, the duration of the tracer injection, and the hydraulic gradient applied to control the infiltrating water fluxes. The method does not require measurement of hydraulic conductivity for the fractured till and enables direct observation of flow and transport patterns within the fractures (e.g., uniform flow and dye tracer distribution, channeling due to aperture variability, and presence of biogenic macropores in fractures). The approach was successfully verified by using the estimated effective fracture aperture values in Large Undisturbed Columns (LUCs) to consistently simulate both the observed LUC effluent breakthrough of a conservative bromide tracer and the water fluxes with the hydraulic gradient applied in the experiments. Sensitivity analyses revealed that estimation of small effective fracture apertures (<10 μm) required accurate determination of the dye tracer retardation factor. By contrast, in the case of larger effective apertures (>20 μm), the sensitivity of the estimated effective fracture aperture to variations in the porous material and solute transport parameters was low compared to the dominant sensitivity to the water flow through the fractures (cubic relation between flow and aperture). The proposed approach may be extended beyond laboratory applications and assist in characterizing field-scale fracture networks.

Seasonal variations in lithium and magnesium isotopic composition of the Nujiang River water at the Daojieba Hydrological Station: Implications for chemical weathering

The response of chemical weathering to climate could be examined in a way by investigating seasonal changes in river geochemistry. Here, we determined the major elements, δ13C of DIC, 87Sr/86Sr, δ7Li, and δ26Mg of monthly time-series water samples at the Daojieba Hydrological Station in lower reach of the Nujiang River over a hydrological year. During high-flow seasons, the dissolved 87Sr/86Sr, δ7Li, and δ26Mg exhibit comparatively high values. 87Sr/86Sr ratio, together with δ13C of DIC, Si/Ca ratio, and pH value of river water, indicate an enhanced silicate weathering contribution to river dissolved load, and river δ7Li is impacted by secondary mineral formation-related fractionation during these periods. Similar to the 87Sr/86Sr ratio, δ26Mg of river water is indicative of a dominant source control. The positive relationship between Mg/Na ratio and δ26Mg value suggests that the prior weathering of Mg-rich silicate minerals could be important for elevated δ26Mg and Mg/Na molar ratio during high-flow seasons, while the Mg isotope fractionation signal can be seen when runoff increases dramatically. During low-flow seasons, river 87Sr/86Sr and δ26Mg values evolve to a carbonate weathering signal, whereas river δ7Li reflects a more congruent silicate weathering, and these changes point to the obvious influence of geothermal water in these periods. The covariations of δ7Li, δ26Mg, and discharge are observed. Changes in runoff can significantly impact the isotopic composition of Li and Mg by changing the hydrological flow path, in which dominant reactions of mineral dissolution and secondary mineral formation vary. This study highlights the hydrological control on Li and Mg isotopic compositions of the Nujiang River water, which could be considered the dominant control of monsoon climate on weathering processes in high-relief river basins at a short-term scale.

Visualizing preferential flow paths using dye tracer and species diversity theory methods to explore their correlation to soil properties with random forest algorithm

The preferential flow path development is potentially the result of spatial variations in soil properties with soil depth. However, visualizing the evolution of the preferential flow path with soil depth remains a challenge. This paper presents dye tracer and species diversity theory methods for characterizing preferential flow paths. Field dye tracer experiments were performed at three sites (tree, bush, and grass) in the Yellow River Delta wetland and dye distribution diversity indices (Simpson index (Ds), Shannon-Wiener index (H), Margalef index (Dm), and Pielou index (E)) were applied to verify their availability for preferential flow assessment. The results showed that the uniformity of the shallow-infiltrated dye at the tree site, non-uniformity of the shallow-infiltrated dye at the bush site, and deep dye infiltration at the grass site were the three typical infiltration types. The average proportion of dye-stained areas (PDA) gradually decreased with increasing soil depth. The quantitative effects of soil properties on PDA changes were profound, indicating that soil clay content at 0–10 cm depth, soil sand content at 10–20 cm depth, soil drainage capacity at 20–30 cm depth, and soil bulk density at 30–40 cm depth were the most predictive factors controlling PDA changes. Our results also showed that dye-stained patches with extremely high and high dye concentrations were the most distributed; Ds, H, Dm, and E were the highest at the tree site and E was the diversity index with the greatest importance for PDA change. The findings reveal the soil properties controlling the formation of preferential flow paths, which will improve our understanding of water resource management in the vadose zones of coastal wetlands.

Design and thermal-hydraulic analysis of multi-target system with 100 MeV proton linear accelerator for the production of 67Cu and 68Ge radioisotopes

Radioisotopes are widely used in the fields of medicine, science, and industry. The growing demand for medical radioisotopes has driven research on alternative production methods. In particular, both isotopes of 67Cu and 68Ge play vital roles in the medical environment in many countries to be used in the radio-immunotherapy and the positron emission tomography imaging, respectively. This study designed a multi-target system consisting of two Zn and one Ga2O3 plates to enable simultaneous production of the medical radioisotopes 67Cu and 68Ge using 100 MeV proton beams. To understand the thermal effect on the multi-targets, we examined the distribution of energy absorbed in each solid plate target when exposed to an accelerated proton beam through the thermal-fluid analysis based on ANSYS simulation. For confirming thermal stability for two Zn targets and one Ga2O3 target, the modified water flow path inside the multi-target system was designed effectively with the controlled distribution of multiple sub-holes between main inlet and the following four channels. It was confirmed that the newly designed multi-target system of Zn and Ga2O3 solid plates shows higher thermal stability than the case of uniform distribution of water inlet, which means it could be exposed to a higher current beam of 7.57% to decrease the processing time.

A five field formulation for flow simulations in porous media with fractures and barriers via an optimization based domain decomposition method

The present work deals with the numerical resolution of coupled 3D–2D problems arising from the simulation of fluid flow in fractured porous media modeled via the Discrete Fracture and Matrix (DFM) model. According to the DFM model, fractures are represented as planar interfaces immersed in a 3D porous matrix and can behave as preferential flow paths, in the case of conductive fractures, or can actually be a barrier for the flow, when, instead, the permeability in the normal-to-fracture direction is small compared to the permeability of the matrix. Consequently, the pressure solution in a DFM can be discontinuous across a barrier, as a result of the geometrical dimensional reduction operated on the fracture. The present work is aimed at developing a numerical scheme suitable for the simulation of the flow in a DFM with fractures and barriers, using a mesh for the 3D matrix non conforming to the fractures and that is ready for domain decomposition. This is achieved starting from a PDE-constrained optimization method, currently available in literature only for conductive fractures in a DFM. First, a novel formulation of the optimization problem is defined to account for non permeable fractures. These are described by a filtration-like coupling at the interface with the surrounding porous matrix. Also the extended finite element method with discontinuous enrichment functions is used to reproduce the pressure solution in the matrix around a barrier. The method is presented here in its simplest form, for clarity of exposition, i.e. considering the case of a single fracture in a 3D domain, also providing a proof of the well posedness of the resulting discrete problem. Four validation examples are proposed to show the viability and the effectiveness of the method.

A distributed parameter model of refrigerant-cooled multi-channel evaporator for battery thermal management

Application of refrigerant-cooled multi-channel evaporators (MCEs) is a promising solution to cool electric vehicle batteries because the non-uniformly distributed heat from a battery can be efficiently dissipated by suitable arrangement of the multi-channels in an MCE. In order to provide a simulation tool for the design of an MCE, a distributed parameter model for quickly predicting the performance of an MCE at various channel arrangements and non-uniformly distributed heating conditions is useful, and thus it is developed in this study. In the model, heat and mass transfer in multiple flow paths in an MCE are separately described; the uneven mass flow distribution among these channels is predicted based on the pressure balance; each flow path is divided into multiple control volumes to reflect the non-uniformly distributed heating conditions; and the mass flow distribution among flow paths and the heat transfer in all the control volumes are jointly solved by a combined iteration algorithm of pressure drop and heat transfer. The proposed model is validated by experiments under the condition of battery heat flux ranging from 500 to 1500 W m−2. The results show that the prediction deviations of the average temperature and the highest temperature on the MCE surface from the experimental ones are 1.5 °C and 1.6 °C, respectively, which indicates that the proposed model can well predict the thermal performance of MCE under non-uniformly distributed battery heat.

Seasonal and spatial variations of water recharging in the Yangtze River using hydrogen and oxygen stable isotopes

Water and material cycles in rivers are controlled by diverse sources of recharge and flow paths, which exhibit spatial and temporal variations in the contributions of water sources and hydrologic processes. To assess the variation of water recharging, the stable isotope abundances of hydrogen and oxygen (δ2H and δ18O) of 152 water samples collected from the Yangtze River mainstream during May (the flat season) and September (the wet season) in 2021 were analyzed. Bayesian modelling was used to calculate the contribution rates of precipitation, glacier snow meltwater, and groundwater. The isotopic compositions of precipitation and the river exhibited comparable seasonal and spatial characteristics. The Bayesian model results showed that the mainstream was governed by precipitation, with a greater contribution of precipitation in the lower reaches during the wet season. Lakes had a larger contribution during the flat season, as evidenced by the greater increase in the magnitude of δ18O values in the mid-lower reaches during this period. The mid-lower reaches of the river had a decreased river evaporation line (REL) slope, which could be owe to lake replenishment. In addition to intensifying the seasonal variations in the contribution ratios of the lakes, the emplacement of the Three Gorges Dam (TGD) reduced the groundwater contribution of the reservoir. This study enhances our understanding of the hydrological and material cycles in the Yangtze River Basin and provides vital baseline knowledge on the contribution of different recharge sources in the mainstream.

A novel compact supercritical CO2 solar receiver based on impinging jet: Heat transfer performance and thermal stress analysis

The compact solar receiver can withstand high temperature and high pressure, which provides promising solutions for direct supercritical CO2 (S-CO2) solar receiver applications. However, it still suffers from local high temperature and large thermal stress due to the non-uniform heating boundary condition. The impinging jet method can significantly enhance the local heat transfer coefficient, holding significant potential in reducing the local high temperature of the compact S-CO2 receiver. Inspired by this, a novel design for a compact S-CO2 receiver using the impinging jet was proposed in this work. The heat transfer performance and stress behavior of the jet channel within the receiver were numerically investigated. The results showed that the jet channel has a significantly higher heat transfer performance than the rectangular channel, and the shorter flow path of the jet channel also effectively reduces the pressure loss. Meanwhile, the jet pattern significantly enhances the temperature uniformity, thereby reducing the thermal stress. In addition, increasing the channel aspect ratio improves the thermo-hydraulic performance, but increase the stress at the same time. The increase in the number of jet holes and the jet hole diameters not only enhance the thermo-hydraulic performance, but also reduce the stress of the jet channel.

3D‐Printed Porous Organic Cages for Gas Filtration: Fabrication and Flow Simulations

AbstractPorous organic cages (POCs) are a class of emerging porous materials with high porosity and selectivity for gas adsorption and separation. As a representative of first‐generation POCs, the cage molecule CC3, has huge potential for separating noble gases and volatile organic compounds (VOCs). 3D printing CC3 using direct ink writing (DIW) offers significant advantages to build complex structures. Through rational design, formulations that are both printable and functional are achieved. Optimised formulations have the elasto‐visco‐plastic behavior required while retaining the functional properties of the cages. The characterization of printed parts evidence that the CC3 structure and adsorption capacity are well preserved. BET surface areas vary with CC3 concentration and can reach up to ≈249 m2 g−1 for 70 wt% CC3. However, high CC3 content leads to a decrease in shrinkage, higher porosity, and low compressive strength ≈0.7 MPa. Computational fluid dynamics (CFD) is used to investigate the gas flow behavior through grid‐type structures with tailored geometric parameters. The results reveal that changing the offset distance and porosity manipulates the flow path, velocity distribution, and pressure drop. This work provides a pathway to design and fabricate structures with multi‐scale porosity that can be widely used for other functional materials and applications.

Contributions of soil organic carbon-induced root- and soil properties complexity to water flow in eastern China

Water flow within the soils affects the efficiency of materials transfer mediated by water. Soil organic carbon (SOC) as an important role in active water flow events can drive the complexity of root-soil synthesis by improving root and soil properties. However, contributions of SOC-induced root- and soil properties complexity to water flow are not well understood. In this study, dye tracing experiments at the three forest stands (oak, pine, and bamboo forests) were conducted to explore water flow patterns, i.e., preferential flow paths (PFP), stream buffer zones (SBZ), and water flow zones (WFZ). X-ray microtomography (CT) scanning was performed to reconstruct the root architecture. The partial least squares path model was applied to quantitatively explore the effects of root- and soil properties on water flow. The results showed that the index of water flow connectivity (IWFC) in the PFP and WFZ patterns decreased with increasing soil depth, while IWFC in the SBZ pattern increased at first and then decreased. In the PFP pattern, soil physical properties had the larger total effects (TE = 0.624) on IWFC change compared with root properties (TE = 0.257). In the SBZ pattern, the total effects of root properties controlling IWFC change (TE = 0.510) were greater than soil physical properties (TE = −0.386). Both of them can equally affect the IWFC in the WFZ pattern. In conclusion, the influences of SOC by driving the changes of soil properties on gravity-driven convective flow process were dramatically stronger than root properties, while SOC could primarily drive the changes of root properties and thereby affect capillary-driven convective flow process. The present results can provide a scientific basis for sustainable forestry management and also a better understanding of the forestry hydrology.

Real-time adaptive scheduling optimization for inter-satellite contact window resources in dynamic satellite networks

Real-time scheduling of contact window resources is a fundamental technology for the future construction of Space–Terrestrial Integrated Networks (STIN). The spatio-temporal dynamics of satellite networks result in continuous variations in satellite relative positions and inter-satellite visibility window resource status. Furthermore, frequent changes in task priority make real-time resource scheduling extremely challenging. Centralized online scheduling methods often lead to a significant amount of waiting time, while distributed online scheduling methods struggle to adapt to the high dynamics of satellite networks. Therefore, the real-time scheduling technology for contact window resources in satellite networks is still in the exploration phase. We propose an inter-satellite online resource allocation method to solve the latency issue caused by the seizure of contact window resources. First, we determine the shortest and maximum flow paths set in the satellite network according to task requirements and construct them as the core network. Then, we comprehensively analyze the supply and demand situation of the task and the core network at a specific moment and verify the reliability of our minimum cost flow-solving method based on the core network. Moreover, we consider factors such as changes in task priority and edge modifications and validate the correctness of our method to update the minimum cost flow at the next moment based on local network changes. Finally, through extensive comparative simulations, we evaluate the absolute advantage of our inter-satellite online scheduling method and validate a significant improvement in task scheduling timeliness.

Is order-2 proportionality good enough for approximating the most likely path flow in user equilibrium traffic assignment?

The proportionality condition is a standard approach to dealing with the non-uniqueness issue in the user equilibrium (UE) traffic assignment problems (TAP). Although the proportionality condition can reduce the degree of arbitrariness, it remains unclear how much arbitrariness remains and whether it can meaningfully affect model outcomes and relevant decisions that depend on them. The answers to these questions are impeded by the lack of an efficient algorithm that can find the exact maximum entropy UE path flow solution for networks of practical size. In this paper, we fill this gap by developing a high-performance augmented Lagrangian algorithm that effectively exploits the special problem structure. Our numerical results reveal that there are a considerable number of links with non-negligible arbitrariness in the solution generated by the proportionality condition, and that this problem becomes worse if the level of congestion increases in the network. Since about a decade ago, many practitioners have relied on state-of-the-art traffic assignment tools based on the proportionality condition to perform select link analysis, among other applications. The results reported herein are a reminder that their toolbox may need reevaluation and perhaps an upgrade.

Hydrogeochemical evolution and isotopic characteristics of groundwater-surface water interaction processes within the Kawere catchment and its hydrogeotechnical implication on pit slope stability at the Nsuta mine, Ghana

In geo-hazard prone heterogenous mining terrain with multiple geological structures, surface water may interact with its underlying groundwater system conveying large quantities of water to deeper aquifers, which could undoubtedly lead to pit-inflows and slope stability challenges for any close-by mining operation. A total of forty-three (43) water samples were collected from the atmosphere, surface, subsurface and the pit-water during dry and wet seasons, respectively to identify pit-inflow origin using hydrogeochemical evolution paths and end member mixing analysis of the major ions at the only manganese mine in Nsuta, Ghana. The major ions order was Na2+>Ca2+>Mg2+>K+ for cations and anions' order was HCO3− > SO42− > Cl−. The identified hydrogeochemical facies and proportions are Ca–HCO3- (68.75%), Ca–Cl2 (18.75%) and mixed Ca–Mg–Cl (12.5%) in the dry season and mixed Ca–Na–HCO3 (54%) followed by Ca–HCO3 (46%) in the wet season. Na+ is the preeminent cation whilst HCO3− is the preeminent anion. Saturation Index (SI) of carbonate and evaporate minerals, particularly calcite, identified groundwater-surface water (GW-SW) exchange flux pathways in the study area, and thus establishing geochemical modeling as a robust tool for identifying GW-SW exchange flux pathways in riverine and aquifer systems. The isotopic composition of the various water samples is depleted having strong mixing relationship of R2 = 0.96%, suggesting the waters have same origin with low d-excess indicating rainfall under lower humidity. End member mixing analysis suggested GW-SW interaction with the surface water contributing an average of 87% to groundwater in the mine pit during the dry season. The generally low hydrogeochemical and isotopic evidence of the various water sources indicate homogeneity pointing to the existence of GW-SW interaction. The fast-movement of the exchange flux has widened flow path leading to pit slope stability issues and huge pit-inflows. Therefore, the mine managers, stakeholders, and decision makers are strongly advised to promptly integrate hydrogeotechnical assessments into their operations as a sustainable remedial measure. This will enable them to continuously reduce the potentiometric head of the groundwater system and ensure cost-efficient mining activities. It is essential to take these actions due to the mine's location within two river stream catchments.