Numerical Unsaturated Flow Model of Railway Drainage Systems

Substandard drainage assets are considered to be a major cause of flooding, earthwork failures, and deficient track geometry. Considering the deterioration of track materials due to cyclic loads and tamping forces, the impact of more frequent extreme rainfall events is likely to lead towards higher rates of hydraulic overloads in the drainage system, earthwork failures, and service disruptions. Therefore, the development of a numerical model could be able to describe the ageing track bed materials and provide an alternative tool for the simulation of the flow through the porous media used in the construction of railway tracks. In this paper the model HYDRUS is tested to simulate the drainage of trackbed materials under laboratory controlled conditions prior its application on actual railway drainage case studies

Lime Cake as an Alternative Stabiliser for Loose Clayey Loams

Lime Cake (precipitated calcium carbonate PCC), a by-product of sugar production, is proposed as a stabiliser for improvement of loose silty clayey loams. Two inorganic pedogenic and organic precipitated calcium carbonate polymorphs are artificially synthesized into a base loosely compacted loamy soil. Formation, micromorphology, quality of cementing bonds, and physiochemical interactions in the interlayer are modelled at molecular level and verified by a suite of micro-analytical spectrometry techniques. Emphasis is put into determining the impacts of polysaccharides on soil strength and implications on soil pore anatomy. Erodibility, compressibility, volumetric change, and hydro-mechanical behaviour of base, and modified soils at yield and post-yield states are studied. Anomalies in suction-controlled post-yield stress–strain behaviour of modified soils are discussed and explained within the tenets of mechanics of composite soils with double porosity. PCC-reinforcement offers the closest possible packing at optimum water content. Desiccation cracking remains likely, but at relatively higher lower-bound water contents. Under low confinement levels and unsaturated state, strain-hardening prevails. Loss of shear strength on saturation is minimal. When saturated, PCC-reinforced soil develops substantially high levels of shear strength at all strain levels. Higher levels of confinement are needed for organic fibrous and onion-skin coating matters to effectively encrust the soil pore network; such high levels, however, leads to formation of an unwelcomed brittle, strain–softening stress–stress behaviour

A two-phase flow model to simulate mold filling and saturation in Resin Transfer Molding

The final publication is available at Springer via http://dx.doi.org/10.1007/s12289-015-1225-zThis paper addresses the numerical simulation of void formation and transport during mold filling in Resin Transfer Molding (RTM). The saturation equation, based on a two-phase flow model resin/air, is coupled with Darcy s law and mass conservation to simulate the unsaturated filling flow that takes place in a RTM mold when resin is injected through the fiber bed. These equations lead to a system composed of an advection diffusion equation for saturation including capillary effects and an elliptic equation for pressure taking into account the effect of air residual saturation. The model introduces the relative permeability as a function of resin saturation. When capillary effects are omitted, the hyperbolic nature of the saturation equation and its strong coupling with Darcy equation through relative permeability represent a challenging numerical issue. The combination of the constitutive physical laws relating permeability to saturation with the coupled system of the pressure and saturation equations allows predicting the saturation profiles. The model was validated by comparison with experimental data obtained for a fiberglass reinforcement injected in a RTM mold at constant flow rate. The saturation measured as a function of time during the resin impregnation of the fiber bed compared very well with numerical predictions.The authors acknowledge financial support of the Spanish Government (Projects DPI2010-20333 and DPI2013-44903-R-AR), of the National Science and Research Council of Canada (NSERC) and of the Canada Reseach Chair (CRC) program.Gascón Martínez, ML.; García Manrique, JA.; Lebel, F.; Ruiz, E.; Trochu, F. (2016). A two-phase flow model to simulate mold filling and saturation in Resin Transfer Molding. International Journal of Material Forming. 9(2):229-239. doi:10.1007/s12289-015-1225-zS22923992Patel N, Lee LJ (1996) Modeling of void formation and removal in liquid composite molding. Part I: wettability analysis. Polym Compos 17(1):96–103Ruiz E, Achim V, Soukane S, Trochu F, Bréard J (2006) Optimization of injection flow rate to minimize micro/macro-voids formation in resin transfer molded composites. Compos Sci Technol 66(3–4):475–486Trochu F, Ruiz E, Achim V, Soukane S (2006) Advanced numerical simulation of liquid composite molding for process analysis and optimization. Compos A: Appl Sci Manuf 37(6):890–902Park CH, Lee W (2011) Modeling void formation and unsaturated flow in liquid composite molding processes: a survey and review. J Reinf Plast Compos 30(11):957–977Pillai KM (2004) Modeling the unsaturated flow in liquid composite molding processes: a review and some thoughts. J Compos Mater 38(23):2097–2118Breard J, Saouab A, Bouquet G (2003) Numerical simulation of void formation in LCM. Compos A: Appl Sci Manuf 34:517–523Breard J, Henzel Y, Trochu F, Gauvin R (2003) Analysis of dynamic flows through porous media. Part I: comparison between saturated and unsaturated flows in fibrous reinforcements. Polym Compos 24(3):391–408Parnas RS, Phelan FR Jr (1991) The effect of heterogeneous porous media on mold filling in Resin Transfer Molding. SAMPE Q 22(2):53–60Parseval DY, Pillai KM, Advani SG (1997) A simple model for the variation of permeability due to partial saturation in dual scale porous media. Transp Porous Media 27(3):243–264Pillai KM (2002) Governing equations for unsaturated flow through woven fiber mats. Part 1. Isothermal flows. Compos A: Appl Sci Manuf 33(7):1007–1019Simacek P, Advani SG (2003) A numerical model to predict fiber tow saturation during Liquid Composite Molding. Compos Sci Technol 63:1725–1736García JA, Gascón L, Chinesta F (2010) A flux limiter strategy for solving the saturation equation in RTM process simulation. Compos A: Appl Sci Manuf 41:78–82Chui WK, Glimm J, Tangerman FM, Jardine AP, Madsen JS, Donnellan TM, Leek R (1997) Process modeling in Resin Transfer Molding as a method to enhance product quality. SIAM Rev 39(4):714–727Nordlund M, Michaud V (2012) Dynamic saturation curve measurement for resin flow in glass fibre reinforcement. Compos A: Appl Sci Manuf 43:333–343García JA, Ll G, Chinesta F (2003) A fixed mesh numerical method for modelling the flow in liquid composites moulding processes using a volume of fluid technique. Comput Methods Appl Mech Eng 192(7–8):877–893García JA, Ll G, Chinesta F, Trochu F, Ruiz E (2010) An efficient solver of the saturation equation in liquid composite molding processes. Int J Mater Form 3(2):1295–1302Lebel F (2012) Contrôle de la fabrication des composites par injection sur renforts. École Polytechnique de Montréal, CanadaVan Genuchten MT (1980) Closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J 44(5):892–898Buckley SE, Leverett MC (1942) Mechanism of fluid displacement in sands. Pet Trans AWME 146:107–116Lundstrom TS, Gebart BR (1994) Influence from process parameters on void formation in Resin Transfer Molding. Polym Compos 15(1):25–33Lundstrom TS (1997) Measurement of void collapse during Resin Transfer Molding. Compos A: Appl Sci Manuf 28(3):201–214Lundstrom TS, Frishfelds V, Jakovics A (2010) Bubble formation and motion in non-crimp fabrics with perturbed bundle geometry. Compos A: Appl Sci Manuf 41:83–92Lebel F, Fanaei A, Ruiz E, Trochu F (2012) Experimental characterization by fluorescence of capillary flows in the fiber tows of engineering fabrics. Open J Inorg Non-Metallic Mater 2(3):25–45Brooks RH, Corey AT (1964) Hydraulic properties of porous media. Colorado State University. Hydrology Papers 1–37Corey AT (1954) The interrelation between gas and oil relative permeabilities. Prod Monthly 19(1):38–4

Removal of hexavalent chromium of contaminated soil by coupling electrokinetic remediation and permeable reactive biobarriers

PURPOSE: In this study, a novel and ecological alternative have been developed to treat soils contaminated with hexavalent chromium coupling two well-known systems: electrokinetic remediation and permeable reactive biobarriers. The electric field promotes the electromigration of the hexavalent chromium oxyanions towards the anode. The biobarriers were placed before the anode electrode, in order to promote the reduction and retention of the chromium migrating in its direction. Thus, this technology provided a global treatment to soil removal without subsequent treatments of the contaminated effluents. METHODS: The electrokinetic system was coupled with two different permeable reactive biobarriers composed by Arthrobacter viscosus bacteria, supported either in activated carbon or zeolite. An electric field of 10 V was applied and two different treatment times of 9 and 18 days were tested. RESULTS: Removal values of 60% and 79% were obtained when electrokinetic treatment was coupled with zeolite and activated carbon biobarriers, respectively, for a test period of 18 day. The reduction of hexavalent chromium to trivalent chromium was around 45% for both systems. CONCLUSIONS: In this work, two types of biobarriers were efficiently coupled to electrokinetic treatment to decontaminate soil with Cr(VI). Furthermore, the viability of the new coupling technology developed (electrokinetic + biobarriers) to treat low-permeability polluted soils was demonstrated.This work was supported by the Spanish Ministry of Science and Innovation (CTQ2008-03059/PPQ), Xunta de Galicia (08MDS034314PR). The authors are grateful to the Spanish Ministry of Science and Innovation for providing financial support for Marta Pazos under the Ramon y Cajal program and the Fundacao para a Ciencia e Tecnologia, Ministerio da Ciencia e Tecnologia, Portugal through the PhD grant of Bruna Fonseca (SFRH/BD/27780/2006)

Mechanistic framework to link root growth models with weather and soil physical properties, including example applications to soybean growth in Brazil

Background and aimsRoot elongation is generally limited by a combination of mechanical impedance and water stress in most arable soils. However, dynamic changes of soil penetration resistance with soil water content are rarely included in models for predicting root growth. Better modelling frameworks are needed to understand root growth interactions between plant genotype, soil management, and climate. Aim of paper is to describe a new model of root elongation in relation to soil physical characteristics like penetration resistance, matric potential, and hypoxia.MethodsA new diagrammatic framework is proposed to illustrate the interaction between root elongation, soil management, and climatic conditions. The new model was written in Matlab®, using the root architecture model RootBox and a model that solves the 1D Richards equations for water flux in soil. Inputs: root architectural parameters for Soybean; soil hydraulic properties; root water uptake function in relation to matric flux potential; root elongation rate as a function of soil physical characteristics. Simulation scenarios: (a) compact soil layer at 16 to 20 cm; (b) test against a field experiment in Brazil during contrasting drought and normal rainfall seasons.Results(a) Soil compaction substantially slowed root growth into and below the compact layer. (b) Simulated root length density was very similar to field measurements, which was influenced greatly by drought. The main factor slowing root elongation in the simulations was evaluated using a stress reduction function.ConclusionThe proposed framework offers a way to explore the interaction between soil physical properties, weather and root growth. It may be applied to most root elongation models, and offers the potential to evaluate likely factors limiting root growth in different soils and tillage regimes

Long-term TNT and DNT contamination: 1-D modeling of natural attenuation in the vadose zone: case study, Portugal

The vadose zone of a trinitrotoluene (TNT) and dinitrotoluene (DNT) contaminated site was investigated to assess the mobility of those explosives under natural conditions. Located in the left margin of the River Tejo Basin, Portugal, the site is located on unconsolidated sediments. Wastewaters associated with the 50-year explosives production were disposed in excavated ponds, from where water would infiltrate and pollute the unsaturated and saturated parts of the local aquifers. Two boreholes were drilled to 9 m depth in such a former waste pond to investigate the contaminant's fate in the vadose zone. Sediment samples were taken every 1-2 m for analysis of the polynitroaromatics (p-NACs) and organic volatile compounds, pH, organic carbon content, cation exchange capacity and grain size analysis. The main contaminant was TNT representing >70 % of the total p-NACs concentration that peaked approximately 7 mg/kg in one borehole, even if the median in both boreholes was of similar to 1 mg/kg. DNT was 4-30 % of the total p-NACs and nitrotoluene (NT), up to 5 %. No other (volatile) organic compound was detected. The predominance of TNT as the main contaminant implies that any natural mass reduction has been inefficient to clean the site. Several 1-D model simulations of p-NACs cleaning of the vadose zone under natural conditions indicated that the most probable scenario of combined advection and partitioning will only remove TNT after 10's of years, whereas DNT and NT will hardly be removed. Such low concentrations and long times for the p-NACs removal, suggest that by now those compounds have been washed-out to a level below standard limits

Improved Interpretation of Mercury Intrusion and Soil Water Retention Percolation Characteristics by Inverse Modelling and Void Cluster Analysis

This work addresses two continuing fallacies in the interpretation of percolation characteristics of porous solids. The first is that the first derivative (slope) of the intrusion characteristic of the non-wetting fluid or drainage characteristic of the wetting fluid corresponds to the void size distribution, and the second is that the sizes of all voids can be measured. The fallacies are illustrated with the aid of the PoreXpert® inversemodelling package.Anewvoid analysis method is then described, which is an add-on to the inverse modelling package and addresses the second fallacy. It is applied to three widely contrasting and challenging porous media. The first comprises two fine-grain graphites for use in the next-generation nuclear reactors. Their larger void sizes were measured by mercury intrusion, and the smallest by using a grand canonical Monte Carlo interpretation of surface area measurement down to nanometre scale. The second application is to the mercury intrusion of a series of mixtures of ground calcium carbonate with powdered microporous calcium carbonate known as functionalised calcium carbonate (FCC). The third is the water retention/drainage characteristic of a soil sample which undergoes naturally occurring hydrophilic/hydrophobic transitions. The first-derivative approximation is shown to be reasonable in the interpretation of the mercury intrusion porosimetry of the two graphites, which differ only at low mercury intrusion pressures, but false for FCC and the transiently hydrophobic soil. The findings are supported by other experimental characterisations, in particular electron and atomic force microscopy

Modelling of Short-Term Interactions Between Concrete Support and the Excavated Damage Zone Around Galleries Drilled in Callovo–Oxfordian Claystone

peer reviewedProduction of energy from nuclear power plants generates high-level radioactive nuclear waste, harmful during dozens of thousand years. Deep geological disposal of nuclear waste represents the most reliable solutions for its safe isolation. Confinement of radioactive wastes relies on the multi-barrier concept in which isolation is provided by a series of engineered (canister, backfill) and natural (host rock) barriers. Few underground research laboratories have been built all over the world to test and validate storage solutions. The underground drilling process of disposal drifts may generate cracks, fractures/strain localisation in shear bands within the rock surrounding the gallery especially in argillaceous rocks. These degradations affect the hydro-mechanical properties of the material, such as permeability, e.g. creating a preferential flow path for radionuclide migration. Hydraulic conductivity increase within this zone must remain limited to preserve the natural barrier. In addition galleries are currently reinforced by different types of concrete supports such as shotcrete and/or prefab elements. Their purpose is twofold: avoiding partial collapse of the tunnel during drilling operations and limiting convergence of the surrounding rock. Properties of both concrete and rock mass are time dependent, due to shotcrete hydration and hydromechanical couplings within the host rock. By the use of a hydro-mechanical coupled Finite Element Code with a Second Gradient regularization, this paper aims at investigating and predicting support and rock interactions (convergence, stress field). The effect of shotcrete hydration evolution, spraying time and use of compressible wedges is studied in order to determine their relative influence