Transport Model Research Articles

Progress in Researching the Ability of Geopolymer Concrete to Withstand a Penetration of Chloride Ions

Geopolymer concrete is an innovative environmentally friendly construction material, and the transportation of chloride ions plays a crucial role in determining its durability. This study provides a summary of the characteristics and limitations of the test techniques used to measure the resistance of geopolymer concrete to the permeability of chloride ions, based on the introduction of the chloride ion transport mechanism in geopolymer concrete. This text provides an overview of the features and constraints of the test techniques used to assess the resistance of geopolymer concrete to chloride ion permeability. It also explores the connections between the mechanism of chloride ion transport and the resistance of geopolymer concrete to chloride ion permeability. This paper provides a concise overview of the properties and constraints of the test methods used to measure the resistance of geopolymer concrete to chloride ion permeability. It also discusses the factors that can affect the chloride ion permeability resistance of geopolymer concrete and presents a comparison between different methods. The article continues by highlighting that the chloride transport model of geopolymer concrete is complex. The essay continues by highlighting the chloride transport model of geopolymer concrete, specifically focusing on the impact of individual parameters such as high temperature, freezing-thaw cycles, and the resistance of geopolymer concrete to chloride ion permeability. The study investigates the impact of freeze-thaw cycles, alkali admixture, and water glass modulus on the resistance of geopolymer concrete to chloride penetration. The infiltration of chloride, as well as the precision of determining the concentration border of chloride ions for colour rendering, require further in-depth investigation.

Analytical model for two-dimensional contaminant transport in a cut-off wall and aquitard system

The cut-off wall is an effective countermeasure to prevent contaminant transport from a contaminated site. However, the aquitard beneath the cut-off wall may serve as a preferential path, leading to an earlier breakthrough of contaminants in the cut-off wall system. In this study, we present a 2-D analytical model for characterizing the contaminant transport in the cut-off wall and aquitard system. The Green function method and discretization method are used to investigate contaminant migration in a semi-infinite domain. We show that the one-dimensional cut-off wall model may be insufficient to predict contaminant transport when the aquitard hydraulic conductivity (k2) exceeds 5 × 10-10 m/s. The contaminant mass outflow from the aquitard for the case with k2 = 2 × 10-9 m/s may account for 53 % of the total contaminant mass outflow. Increasing k2 from 5 × 10-10 m/s to 5 × 10-9 m/s leads to a 2.9-fold increase of cumulative mass outflow. An approximate fitting formula is proposed to calculate the required minimum thickness of cut-off wall under different type of aquitards. The thicker cut-off wall is required to ensure the cumulative mass outflow will not exceed the allowable value when k2 is relatively high (e.g., 1 × 10-9 m/s).

Quantifying the Contribution of Multiple Processes to the Dust Decreasing Trend in the Guliya Ice Core Over the Past 50 Years

AbstractDust records extracted from ice cores can facilitate the reconstruction of historical atmospheric dust levels and climate change. However, interpreting dust variations in ice cores is intricate because of the compounded influence of emission, transport, and deposition processes. This study investigated dust records retrieved from the Guliya ice cap drilled in 2015 on the West Tibetan Plateau using a mean trajectory transport and deposition model. Results showed that the Guliya dust concentration has exhibited a declining trend since the 1960s (−751 μg kg−1 yr−1). Applying an attribution approach, we discovered that low dust emission (80.3%) was the main cause of the drop in dust concentration, with changes related to transportation (5.2%) and deposition (14.5%) making only minor contributions. The weakening of surface wind speed in the desert and increasing precipitation in both the desert and glacier were the primary factors driving the decrease in Guliya dust concentration.

Modeling and Control of Strip Transport in Metal Peeling

Modeling and Control of Strip Transport in Metal Peeling

Numerical simulation of deuterium retention in tungsten under ELM-like conditions

Retention of hydrogen isotopes (HI) in plasma-facing components (PFCs) is a crucial process influencing the operation of a fusion device. The dynamics of HI retention is mainly determined by irradiation conditions of PFCs. These conditions can change due to the onset of fast transient events, such as edge-localised modes (ELMs). The development of ELMs results in repetitive short-term plasma bursts on PFCs, affecting the exposure regime. In this work, the effect of ELM-like loads on the fuel retention is numerically analysed using a one-dimensional diffusion model, implemented in the FESTIM code. As a representative simulation case, the deuterium (D) retention in tungsten (W) is considered with the geometry and inter-ELM exposure conditions relevant for large fusion devices. Temporal dependencies of heat and particle fluxes during intra-ELM stages are accounted for by applying the free-streaming model for an ELM filament transport. The simulations were conducted for three inter-ELM exposure regimes with various properties of transient loads, D trapping sites, and D recombination rates. Compared to the ELM-free irradiation, the onset of transients is shown to mainly reduce the D retention rate because of significant material heating induced by the arrival of energetic ELMy particles. This relative difference can decrease if the material is characterised by strong trapping sites or limited desorption from a surface. The findings also demonstrate that additional heating during transients provides conditions for a faster D migration into the bulk. In addition, we discuss the possibility of using ELM-average loads to obtain quick assessments of the D content. The approach is shown to be applicable for the case of small ELMs and is used to derive the analytical expressions for the D distribution in W within the steady-state approximation.

Crude oil-water emulsions formed by nozzle type inflow control device

Conventional water-in-oil emulsion studies where shear energy intensity and duration are applied to achieve ‘equilibrium’ emulsion droplet properties may not represent partial emulsion conditions observed in single pass flow through Inflow Control Device (ICD) orifices. Our experiments provide data expected to improve emulsion transport models for partial emulsification processes expected in well completion and production tubulars. Development of a novel flow adaptor apparatus is described together with a laboratory homogeniser enables independent control of shear intensity and duration together with input water volume fraction. The emulsions formed are characterised with a benchtop Nuclear Magnetic Resonance spectrometer for water content and droplet size, viscosity by a rotating bob viscometer and free water/emulsion proportion by optical photography. Experimental results confirmed partial emulsification above 30% V/V, as observed in both previous ICD flow experiments and field observations. Electrical conductivity measurements of the emulsions produced indicated that partial emulsification comprising an oil external emulsion phase together with free water, rather than emulsion inversion, occurs in this system.

Investigation on aeration efficiency and energy efficiency optimization in recirculating aquaculture coupling CFD with Euler-Euler and species transport model

Aeration plays a crucial role in maintaining the overall health and sustainability of the recirculating aquaculture system. Conventional aeration operates continuous and consumes 35.06 % of equipment energy consumption. How to optimize aeration while ensuring normal fish growth remains a great challenge. To address these challenges, this research paper established a coupled model of computational fluid dynamics and species transport model (CFD-STM). Aeration tests were carried out in a 33.3 m3 culture tank using a 0.31 m diameter aerator to find out the impacts of operating variables (aeration flow rate, bubble diameter, inlet velocity, and aerator position) on oxygen transfer using a microporous aeration system. A CFD-STM model was coupled to obtain the optimum values of aeration parameters. Results revealed that the Standard Aeration Efficiency (SAE) reductions by 32.1 % for 50 % larger aeration flow rate and 86.4 % for 50 % larger bubble diameter, as well as SAE achieved 5.67 kg/kWh with an aerator position scale factor equal to 0.48. In addition, the energy consumption of the Roots blower was compared for a single aeration and dissolved oxygen saturation. Results indicated that smaller bubbles microbubbles can effectively improve oxygen transfer coefficient. The energy consumption was reduced by 15.7 % for 50 % larger aeration flow rate and 38.6 % for 50 % lower bubble diameter. Based on the results, this study demonstrated that CFD is a useful tool to optimize the aeration system design and operation in order to enhance aeration efficiency. It provides a theoretical foundation the future optimization of energy conservation through intermittent aeration practices.

Surrogate based design space exploration and exploitation for an efficient airfoil optimization under uncertainties using transition models

The pursuit of sustainable, zero-emission air travel is heavily dependent on the creation of energy-efficient aircraft. Key strategies for achieving this sustainability in aviation include reducing fuel consumption through low-drag designs harnessing laminar flow. However, designing aircraft with laminar flow characteristics is complex due to their sensitivity to environmental and operational factors. This study tackles the challenge of developing energy-efficient aircraft by using computational fluid dynamics models and sophisticated optimization techniques that account for uncertainty. Our approach demonstrates the effectiveness of surrogate-based optimization and uncertainty quantification in optimizing airfoil drag for a natural laminar airfoil (NLF) design. We use surrogate models, trained with data from detailed airfoil simulations, which include a boundary layer code coupled with a linear stability method and a newly developed transition transport model. Transition location predicted using transition models facilitate an accurate drag prediction used in the optimization process. The accuracy of these surrogate models is enhanced through active sampling strategies. Our robust optimization method considers uncertainties in environmental and operational conditions, offering a deeper insight into their effects on crucial design parameters. Unlike traditional deterministic aerodynamic design optimization, our findings highlight the efficacy and precision of uncertainty-based optimization in achieving robust NLF airfoil designs over large (exploration mode) and small (exploitation mode) design spaces. Investigating design space parameterization based on the size of design variables reveals significant differences in optimal airfoil configurations. The optimized designs we propose favor delayed transition, in contrast to deterministic designs which often result in significant loss of laminarity when facing uncertainties. This study represents a significant advancement in aerospace engineering, providing a practical and effective methodology for creating energy-efficient airfoil designs. The application of these advanced optimization and uncertainty quantification techniques shows great potential for the wider field of aerospace engineering, paving the way for more resilient and robust aircraft designs.

Flat-top plasma operational space of the STEP power plant

Abstract STEP is a spherical tokamak prototype power plant that is being designed to demonstrate net electric power. The design phase involves the exploitation of plasma models to optimise fusion performance subject to satisfying various physics and engineering constraints. A modelling workflow, including integrated core plasma modelling, MHD stability analysis, SOL and pedestal modelling, coil set and free boundary equilibrium solvers, and whole plant design, has been developed to specify the design parameters and to develop viable scenarios. The integrated core plasma model JETTO is used to develop individual flat-top operating points that satisfy imposed criteria for fusion power performance within operational constraints. Key plasma parameters such as normalised beta, Greenwald density fraction, auxiliary power and radiated power have been scanned to scope the operational space and to derive a collection of candidate non-inductive flat-top points. The assumed auxiliary heating and current drive is either from electron cyclotron systems only or a combination of electron cyclotron and electron Bernstein waves. At present stages of transport modelling, there is a large uncertainty in overall confinement for relevant parameter regimes. For each of the two auxiliary heating and current drive systems scenarios, two candidate flat-top points have been developed based on different confinement assumptions, totalling to four operating points. A lower confinement assumption generally suggests operating points in high-density, high auxiliary power regimes, whereas higher confinement would allow access to a broader parameter regime in density and power while maintaining target fusion power performance.

Modeling and simulation of simultaneous transport and incompressible flows on all level sets in a bulk domain

Modeling and simulation of simultaneous transport and incompressible flows on all level sets in a bulk domain

Development of novel osteochondral scaffolds and related in vitro environment with the aid of chemical engineering principles

In tissue engineering, collaboration among experts from different fields is needed to design appropriate cell scaffolds and the required three-dimensional environment. Osteochondral tissue engineering is particularly challenging due to the need to provide scaffolds that imitate structural and compositional differences between two neighboring tissues, articular cartilage and bone, and the required complex biophysical environments for cultivating such scaffolds. This work focuses on two key objectives: first, to develop bilayered osteochondral scaffolds based on gellan gum and bioactive glass and, second, to create a biomimetic environment for scaffold characterization by designing and utilizing novel dual-medium cultivation bioreactor chambers. Basic chemical engineering principles were utilized to help achieve both aims. First, a simple heat transport model based on one-dimensional conduction was applied as a guideline for bilayer scaffold preparation, leading to the formation of a gelatinous upper part and a macroporous lower part with a thin, well-integrated interfacial zone. Second, a novel cultivation chamber was developed to be used in a dynamic compression bioreactor to provide possibilities for flow of two different media, such as chondrogenic and osteogenic. These chambers were utilized for characterization of the novel scaffolds with regard to bioactivity and stability under dynamic compression and fluid perfusion over 14 d, while flow distribution under different conditions was analyzed by a tracer method and residence time distribution analysis.

MFS Transporter as the Molecular Switch Unlocking the Production of Cage-Like Acresorbicillinol C.

Sorbicillinoids are a class of fungal polyketides with diverse structures and distinguished bioactivities. Although remarkable progress has been achieved in their chemistry and biosynthesis, the efflux of sorbicillinoids is poorly understood. Here, we found MFS transporter AcsorT was responsible for the biosynthesis of sorbicillinoids in Acremonium chrysogenum. Combinatorial knockout and subcellular location demonstrated that the plasma membrane-associated AcsorT was responsible for the transportation of sorbicillinol and subsequent formation of oxosorbicillinol and acresorbicillinol C via the berberine bridge enzyme-like oxidase AcsorD in the periplasm. Homology modeling and site-directed mutation revealed that Tyr303 and Arg436 were the key residues of AcsorT, which was further explained by molecular dynamics simulation. Based on our study, it was suggested that AcsorT modulates sorbicillinoid production by coordinating its biosynthesis and export, and a transport model of sorbicillinoids was proposed in A. chrysogenum.

Observed APO Seasonal Cycle in the Pacific: Estimation of Autumn O2 Oceanic Emissions

AbstractIn this work, we investigated the seasonal cycle of atmospheric potential oxygen (APO), a unique tracer of air‐sea gas exchanges of molecular oxygen (O2) and carbon dioxide (CO2), expressed as APO = O2 + 1.1 × CO2. APO data were obtained from flask air samples collected since the late 1990s at three Japanese ground stations and on commercial cargo ships sailing between Japan and Australia/New Zealand, North America, and Southeast Asia. We also analyzed the APO spatial distribution and seasonal cycles with simulations from an atmospheric transport model using climatological oceanic O2 fluxes from an empirical product that relate O2 flux to ocean heat as input. Model simulations reproduced the observed APO seasonal cycles generally well, but with larger amplitudes and earlier occurrence of seasonal minima and maxima than in the observations. Moreover, the observed seasonal cycles exhibited larger APO enhancements than the simulations in autumn and early winter, especially in the North Pacific at 20°N–60°N. These enhancements remained when refining the comparison by adjusting the simulated APO peak‐to‐peak amplitudes and seasonal phases to the observations. This suggests additional O2 emissions in the North Pacific, not well expressed in the air‐sea O2 fluxes used as input for our model simulations. The average autumn enhancement at 40°N–60°N was approximately twice that measured at 20°N–40°N. Confirming previous studies, our results indicate two distinct mechanisms possibly contributing to the additional oceanic O2 emissions: outgassing from a subsurface shallow oxygen maximum at 20°N–40°N and autumn phytoplankton bloom at 40°N–60°N.

A pore-scale lattice Boltzmann model for solute transport coupled with heterogeneous surface reactions and mineral dissolution

In this paper, we propose a pore-scale lattice Boltzmann model to simulate fluid flow coupled with heterogeneous surface reactions and mineral dissolution. The primary innovation lies in the transformation of surface reactions, originally treated as a boundary condition, into a volume-source term through dimensional augmentation within the framework of sharp liquid–solid interfaces. This significantly simplifies the implementation, particularly for reactions occurring in porous media with intricate geometric structures. Several benchmark tests were performed to validate the accuracy of this model, including a reaction–diffusion problem in a rectangular domain, a two-dimensional reaction and dissolution of a circular grain, as well as a three-dimensional calcite crystal dissolution in a micro-channel. All the obtained simulation results agree well with the reference solutions. In addition, a dissolution problem in a three-dimensional porous medium built with the sand pack is then investigated. Cases with different Péclet numbers (Pe) and Damköhler numbers (Da) were simulated, and five dissolution modes were obtained, which were finally summarized in a diagram of Pe and Da.

Assimilating Size Diversity: Population Balance Equations Applied to the Modeling of Microplastic Transport.

Modeling of microplastic (MP) transport in the aquatic environment is complicated by the diverse properties of the plastic particles. Traditional modeling methods such as Lagrangian particle tracking and Eulerian discrete class (DC) methods have limitations as they are not best placed to account for the diverse characteristics of individual particles, namely, size, density, and shape, which are crucial for determining the transport of MPs. In this work, we address the issue of particle size diversity by using the population balance equations (PBE) method. In addition to the advection-diffusion terms, the PBE transport equation involves a deposition sink term. Seven size classes of MPs are modeled in the DC method, which is compared to the PBE method. The evolution of particle size distribution is compared between the two methods using a simplified test case of a schematized estuary with tidal forcing and river discharge. This work successfully demonstrates the applicability and appropriateness of the PBE model in modeling the transport of MPs to track the dynamic and complete size distribution at a reduced computational cost in comparison to the DC model. With the PBE method, it is possible to address other diversities of the MPs such as the shape and density.

Shifting sands: The influence of coral reefs on shoreline erosion from short-term storm protection to long-term disequilibrium

Climate change is exacerbating shoreline erosion and flooding, posing significant risks to coastal communities. Although traditional coastal defenses such as seawalls, dykes, and breakwaters offer protection from these hazards, their high environmental and economic costs are driving interest in cost-competitive nature-based solutions. Coral reef restoration is a nature-based solution that may be particularly apt to mitigate tropical coastal flooding and shoreline erosion while providing benefits to local tourism, fisheries, and nature. However, the novelty of this field requires studies demonstrating the benefits of reefs for coastal protection. While the flood protection benefits of reefs have been well-documented, their effects on shoreline erosion are comparatively less understood. Here, we investigate the effects of coral reefs on shoreline erosion by comparing tropical beach responses at short and long timescales, as well as identifying important reef structural features influencing coastal erosion rates. Our analyses leveraged two key datasets created in this study: the first derived from a literature review on short-term shoreline erosion due to storm events, and another compiling >80 years of long-term erosion rates, bathymetry, habitat, and wave energy for the Hawaiian Islands of Kauaʻi, Oʻahu, and Maui. Our analyses reveal three key findings regarding the effects of reefs on shoreline erosion. Firstly, we find evidence for the role of reefs in mitigating shoreline erosion during storm events, with coral reef-protected beaches experiencing 97 % less beach volume loss than unprotected beaches. Secondly, a linear regression analysis demonstrates that coral reef structure and wave energy are important predictors of long-term shoreline erosion rates, explaining 34 % of the variation across the Hawaiian Islands. Consistent with prior research, we find beaches protected by coral reefs with shallow reef crests, wide reef flats, calmer offshore conditions, and positioned farther from the shore exhibit lower erosion rates than others. Finally, when comparing historical erosion rates of protected and unprotected beaches in Hawai'i, we find a seemingly incongruous pattern where coral reef-protected beaches eroded up to 2x faster than beaches without reefs. While the cause of the enhanced erosion is yet to be fully understood, a combination of coral reef structural degradation and sea-level rise is likely shifting the equilibrium profiles of reef-protected beaches inshore. These results emphasize the role of coral reefs in reducing coastal erosion during storm events while revealing contrasting erosion patterns over long timescales. Future studies would ideally broaden the scope to include various regions, utilize advanced sediment transport models, and undertake field experiments to deepen our understanding of coral reef-coupled shoreline dynamics.

Non-uniform distribution model of rust layer around steel bar circumference and its effect on corrosion-induced concrete cracking

This paper presents a critical review of non-uniform distribution models of rust layer around a reinforcing bar’s circumference and a numerical investigation into the effect of non-uniform distribution of rust on (a) crack patterns and (b) development of concrete surface crack width. Different shapes of rust layer simulated by mathematical models were established based on statistical distribution functions which were compared with the transport-electrochemical model and published testing results. In addition, an advanced 2D finite element model was implemented to simulate crack propagation in concrete caused by rebar corrosion. It is found that non-uniform rust distribution shape around the rebar not only affected concrete surface crack width evolution, but also crack length and crack pattern. Transport-electrochemical model and Gaussian distribution-based model could reasonably simulate rust distribution and associated corrosion-induced cracks. In addition, the authors proposed a simplified semi-empirical model that incorporated an analytical chloride transport model and a Gaussian model to simulate time-dependent advancement of non-uniform rust layer. Furthermore, a parametric study of the effect of elastic stiffness of the rust layer on crack patterns and surface crack width development was also performed.

Insights on influence of polymer molecular weight on gas sorption, transport and plasticization in glassy 6FDA-DAM polyimide membranes

It is well-known that spinning of defect-free asymmetric polyimide hollow fibers is promoted by using high polymer molecular weight. On the other hand, effects of polymer molecular weight on microstructure and separation performance of dense film 6FDA-based polyimide membranes relevant to sheath layers of composite membranes are less well-explored. In this work, gas sorption, diffusion and permeation of CO2, N2 and CH4 for 6FDA-DAM dense films are considered for 45 kDa and 176 kDa weight average molecular weights. Such molecular weights are relevant for practical hollow fiber spinning of membranes. Earlier results for polystyrene samples show similar trends in dual mode sorption model results like those noted here for the 6FDA-DAM polyimide case; however, effects on permeation like those considered here were not addressed for the polystyrene case. Here we also address actual pure and mixed gas CO2 and CH4 permeation and show higher free volumes, higher gas permeability but lower 50:50 CO2/CH4 selectivity associated with the higher molecular weight sample. Analysis using the self-consistent dual-mode sorption and transport models help understand the observed trends. Our results suggest denser packing of polymer segments exists in the Henry's law environment of the low molecular weight sample. Moreover, the higher polymer segmental packing in the Henry's law environment suppresses CO2 plasticization with a dual-mode microstructure controlling separation performance and plasticization behavior for this important member of the 6FDA polyimide family. We show evidence of effects of charge transfer complexes as possible main features in these trends. These dense film results help guide future work on dense sheath layers on composite hollow fiber membranes.

A comparative study of shale oil transport behavior in graphene and kerogen nanopores with various roughness via molecular dynamics simulations

Understanding the transport behavior of molecules within nanopores is crucial in various fields, including nanofluidic, bioengineering and geophysics. To elucidate oil flow within shale organic nanopores, the simplified carbon nanopores (e.g., graphene and CNTs) and realistic kerogen nanopores have been extensively applied. The variation in pore structure models leads to divergent findings in oil flow behavior. This study employs molecular dynamics simulations to compare alkane flows within graphene and kerogen nano-slits with similar surface roughness. We observe that using graphene nanopores with ultra-smooth surfaces may introduce errors of several orders of magnitudes in flux calculations. Conversely, when accounting for similar roughness, the static characteristics and flow behaviors in graphene and kerogen slit pores are largely comparable. The relative error of flux and the flow velocity could be lower than 6% when using graphene to replace the realistic kerogen. Alkane molecules exhibit a parabolic velocity profile, with a no-slip condition adjacent to the rough surfaces. In addition, when applying the Hagen-Poiseullie equation for alkane flow in rough nanopores, Gibbs dividing surface is recommended for defining the pore radius, reducing the error of flux predicted to below 10%. Our work fulfills the gap in accessing the suitability of using graphene as a model for oil transport in shale organic nanopores, and provides valuable insights and guidance for future studies on nanofluidic in shale, including molecular simulations, nanofluidic experiment designs, and multi-scale modeling.

The value of simplified models of radionuclide transport for the safety assessment of nuclear waste repositories: A benchmark study

In order to assess sites for a deep geological repository for storing high-level nuclear waste safely in Germany, various numerical models and tools will be in use. For their interaction within one workflow, their reproducibility, and reliability version-controlled open-source solutions and careful documentation of model setups, results and verifications are of special value. However, spatially fully resolved models including all relevant physical and chemical processes are neither computationally feasible for large domains nor is the data typically available to parameterize such models. Thus, simplified models are crucial for the pre-assessment of possible sites to narrow down the list of suitable candidates for which detailed site investigations and fully resolved models will be done at a later stage. Still, the accuracy of these simplified models is of importance as the pre-assessment of suitable sites will be based on them. In this study, we compare the modelling capabilities of TransPyREnd, a one-dimensional transport code based on finite differences, specifically developed for the fast estimation of radionuclide transport by the German federal company for radioactive waste disposal (BGE), with OpenGeoSys, which is a modelling platform based on finite elements in up to three spatial dimensions. Both codes are used in the site selection procedure for the German nuclear waste repository. The comparison of the model results of TransPyREnd and OpenGeoSys is augmented by comparisons with an analytical solution for a homogeneous material. For the purpose of numerical benchmarking, we consider a geological profile located in southern Germany as an example where the hypothetical repository is located in a clay-stone formation. TransPyREnd and OpenGeoSys yield overall similar results. However, both codes use different discretizations which impact is the highest for strongly sorbing compounds, while the difference gets negligible for less sorbing and more diffusive compounds as higher diffusion tends to blur the initial conditions. Overall, the OpenGeoSys model is more exact whereas the TransPyREnd model has considerable faster run times. We found in our example, that significant substance amounts are only leaving the host rock formation, if apparent diffusion is high, for which case both codes give similar results, while relative differences are considerable for strongly sorbing compounds. However, in the latter case no significant substance amount of radionuclides leaves the host-rock formation, thus deeming the differences in the model results minor for the overall safety assessments of sites.