THE TRANSPORT OF HEAT BY FLOWING GROUNDWATER IN A DISCRETE FRACTURE

In typical theoretical or experimental studies of heat migration in discrete fractures, conduction and thermal dispersion are commonly neglected from the fracture heat transport equation, assuming heat conduction into the matrix is predominant. In this study analytical and numerical models are used to investigate the significance of conduction and thermal dispersion in the plane of the fracture for a point and line sources geometries. The analytical models account for advective, conductive and dispersive heat transport in both the longitudinal and transverse directions in the fracture. The heat transport in the fracture is coupled with a matrix equation in which heat is conducted in the direction perpendicular to the fracture. In the numerical model, the governing heat transport processes are the same as the analytical models; however, the matrix conduction is considered in both longitudinal and transverse directions. Firstly, we demonstrate that longitudinal conduction and dispersion are critical processes that affect heat transport in fractured rock environments, especially for small apertures (eg. 100 μm or less), high flow rate conditions (eg. velocity greater than 50 m/day) and early time (eg. less than 10 days). Secondly, transverse thermal dispersion in the fracture plane is also observed to be an important transport process leading to retardation of the migrating heat front particularly at late time (eg. after 40 days of hot water injection). Solutions which neglect dispersion in the transverse direction underestimate the locations of heat fronts at late time. Finally, this study also suggests that the geometry of the heat sources has significant effects on the heat transport in the system. For example, the effects of dispersion in the fracture are observed to decrease when the width of the heat source expands.

Settling velocity of sediments at high concentrations

New data on the settling velocity of artificial sediments and natural sands at high concentrations are presented. The data are compared with a widely used semiempirical Richardson and Zaki equation (Trans. Inst. Chem. Eng. 32 (1954) 35), which gives an accurate measure of the reduction in velocity as a function of concentration and an experimentally determined empirical power n. Here, a simple method of determining n is presented using standard equations for the clear water settling velocity and the seepage flow within fixed sediment beds. The resulting values for n are compared against values derived from new and existing laboratory data for beach and filter sands. For sands, the appropriate values of n are found to differ significantly from those suggested by Richardson and Zaki for spheres, and are typically larger, corresponding to a greater reduction in settling velocity at high concentrations. For fine and medium sands at concentrations of order 0.4, the hindered settling velocity reduces to about 70% of that expected using values of n derived for spheres. At concentrations of order 0.15, the hindered settling velocity reduces to less than half of the settling velocity in clear water. These reduced settling velocities have important implications for sediment transport modelling close to, and within, sheet flow layers and in the swash zone.

Prediction of slurry transport in SAG mills using SPH fluid flow in a dynamic DEM based porous media

DEM modelling of the motion of coarse fractions of the charge inside SAG mills has now been well established for more than a decade. In these models the effect of slurry has broadly been ignored due to its complexity. Smoothed particle hydrodynamics (SPH) provides a particle based method for modelling complex free surface fluid flows and is well suited to modelling fluid flow in mills. Previous modelling has demonstrated the powerful ability of SPH to capture dynamic fluid flow effects such as lifters crashing into slurry pools, fluid draining from lifters, flow through grates and pulp lifter discharge. However, all these examples were limited by the ability to model only the slurry in the mill without the charge. In this paper, we represent the charge as a dynamic porous media through which the SPH fluid is then able to flow. The porous media properties (specifically the spatial distribution of porosity and velocity) are predicted by time averaging the mill charge predicted using a large scale DEM model. This allows prediction of transient and steady state slurry distributions in the mill and allows its variation with operating parameters, slurry viscosity and slurry volume, to be explored. (C) 2006 Published by Elsevier Ltd.

Stochastic modelling and simulations for solute transport in porous media

A stochastic model for solute transport in aquifers is studied based on the concepts of stochastic velocity and stochastic diffusivity. By applying finite difference techniques to the spatial variables of the stochastic governing equation, a system of stiff stochastic ordinary differential equations is obtained. Both the semi-implicit Euler method and the balanced implicit method are used for solving this stochastic system. Based on the Karhunen-Loeve expansion, stochastic processes in time and space are calculated by means of a spatial correlation matrix. Four types of spatial correlation matrices are presented based on the hydraulic properties of physical parameters. Simulations with two types of correlation matrices are presented.

Fibre agglomerate transport in a horizontal flow

An experimental and theoretical study of the transport of mineral wool fibre agglomerates in nuclear power plant containment sumps is being performed. A racetrack channel was devised to provide data for the validation of numerical models, which are intended to model the transport of fibre agglomerates. The racetrack channel provides near uniform and steady conditions that lead to either the sedimentation or suspension of the agglomerates. Various experimental techniques were used to determine the velocity conditions and the distribution of the fibre agglomerates in the channel. The fibre agglomerates are modelled as fluid particles in the Eulerian reference frame. Simulations of pure sedimentation of a known mass and volume of agglomerations show that the transport of the fibre agglomerates can be replicated. The suspension of the fibres is also replicated in the simulations; however, the definition of the fibre agglomerate phase is strongly dependent on the selected density and diameter. Detailed information on the morphology of the fibre agglomerates is lacking for the suspension conditions, as the fibre agglomerates may undergo breakage and erosion. Therefore, ongoing work, which is described here, is being pursued to improve the experimental characterisation of the suspended transport of the fibre agglomerates.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors, Part A: Eulerian computation of momentum transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase, according to the literature. FLUENT 6.2 has been used as the modelling framework of the simulations with a completely revised drag model, in the form of user defined function (UDF), to calculate the forces exerted on the particle as well as its velocity components. 2-D and 3-D simulations are tested and compared. The study is the first part of a complete pyrolysis model in fluidised bed reactors.

Verification and validation of numerical models of the transport of insulation debris

Damage to insulation materials located near to a primary circuit coolant leak may compromise the operation of the emergency core cooling system (ECCS). Insulation material in the form of mineral wool fiber agglomerates (MWFA) maybe transported to the containment sump strainers, where they may block or penetrate the strainers. Though the impact of MWFA on the pressure drop across the strainers is minimal, corrosion products formed over time may also accumulate in the fiber cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz is investigating the phenomena that maybe observed in the containment vessel during a primary circuit coolant leak. The study entails the generation of fiber agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effect that corrosion and erosion of the containment internals by the coolant has on the strainer pressure drop. The focus of this paper is on the verification and validation of numerical models that can predict the transport of MWFA. A number of pseudo-continuous dispersed phases of spherical wetted agglomerates represent the MWFA. The size, density, the relative viscosity of the fluid-fiber agglomerate mixture and the turbulent dispersion all affect how the fiber agglomerates are transported. In the cases described here, the size is kept constant while the density is modified. This definition affects both the terminal velocity and volume fraction of the dispersed phases. Note that the relative viscosity is only significant at high concentrations. Three single effect experiments were used to provide validation data on the transport of the fiber agglomerates under conditions of sedimentation in quiescent fluid, sedimentation in a horizontal flow and suspension in a horizontal flow. The experiments were performed in a rectangular column for the quiescent fluid and a racetrack type channel that provided a near uniform horizontal flow. The numerical models of sedimentation in the column and the racetrack channel found that the sedimentation characteristics are consistent with the experiments. For channel suspension, the heavier fibers tend to accumulate at the channel base even at high velocities, while lighter phases are more likely to be transported around the channel.

Numerical models used for the modelling of the transport of fibrous insulation debris

Mineral wool insulation material applied to the primary cooling circuit of a nuclear reactor maybe damaged in the course of a loss of coolant accident (LOCA). The insulation material released by the leak may compromise the operation of the emergency core cooling system (ECCS), as it maybe transported together with the coolant in the form of mineral wool fiber agglomerates (MWFA) suspensions to the containment sump strainers, which are mounted at the inlet of the ECCS to keep any debris away from the emergency cooling pumps. In the further course of the LOCA, the MWFA may block or penetrate the strainers. In addition to the impact of MWFA on the pressure drop across the strainers, corrosion products formed over time may also accumulate in the fiber cakes on the strainers, which can lead to a significant increase in the strainer pressure drop and result in cavitation in the ECCS. Therefore, it is essential to understand the transport characteristics of the insulation materials in order to determine the long-term operability of nuclear reactors, which undergo LOCA. An experimental and theoretical study performed by the Helmholtz-Zentrum Dresden-Rossendorf and the Hochschule Zittau/Görlitz1 is investigating the phenomena that maybe observed in the containment vessel during a primary circuit coolant leak. The study entails the generation of fiber agglomerates, the determination of their transport properties in single and multi-effect experiments and the long-term effects that particles formed due to corrosion of metallic containment internals by the coolant medium have on the strainer pressure drop. The focus of this presentation is on the numerical models that are used to predict the transport of MWFA by CFD simulations. A number of pseudo-continuous dispersed phases of spherical wetted agglomerates can represent the MWFA. The size, density, the relative viscosity of the fluid-fiber agglomerate mixture and the turbulent dispersion all affect how the fiber agglomerates are transported. In the cases described here, the size is kept constant while the density is modified. This definition affects both the terminal velocity and volume fraction of the dispersed phases. Only one of the single effect experimental scenarios is described here that are used in validation of the numerical models. The scenario examines the suspension and horizontal transport of the fiber agglomerates in a racetrack type channel. The corresponding experiments will be described in an accompanying presentation (see abstract of Seeliger et al.).

Cross-correlation-aided transport in stochastically driven accretion flows

The origin of linear instability resulting in rotating sheared accretion flows has remained a controversial subject for a long time. While some explanations of such non-normal transient growth of disturbances in the Rayleigh stable limit were available for magnetized accretion flows, similar instabilities in the absence of magnetic perturbations remained unexplained. This dichotomy was resolved in two recent publications by Chattopadhyay and co-workers [Mukhopadhyay and Chattopadhyay, J. Phys. A 46, 035501 (2013)1751-811310.1088/1751-8113/46/3/035501; Nath, Phys. Rev. E 88, 013010 (2013)PLEEE81539-375510.1103/PhysRevE.88.013010] where it was shown that such instabilities, especially for nonmagnetized accretion flows, were introduced through interaction of the inherent stochastic noise in the system (even a "cold" accretion flow at 3000 K is too "hot" in the statistical parlance and is capable of inducing strong thermal modes) with the underlying Taylor-Couette flow profiles. Both studies, however, excluded the additional energy influx (or efflux) that could result from nonzero cross correlation of a noise perturbing the velocity flow, say, with the noise that is driving the vorticity flow (or equivalently the magnetic field and magnetic vorticity flow dynamics). Through the introduction of such a time symmetry violating effect, in this article we show that nonzero noise cross correlations essentially renormalize the strength of temporal correlations. Apart from an overall boost in the energy rate (both for spatial and temporal correlations, and hence in the ensemble averaged energy spectra), this results in mutual competition in growth rates of affected variables often resulting in suppression of oscillating Alfven waves at small times while leading to faster saturations at relatively longer time scales. The effects are seen to be more pronounced with magnetic field fluxes where the noise cross correlation magnifies the strength of the field concerned. Another remarkable feature noted specifically for the autocorrelation functions is the removal of energy degeneracy in the temporal profiles of fast growing non-normal modes leading to faster saturation with minimum oscillations. These results, including those presented in the previous two publications, now convincingly explain subcritical transition to turbulence in the linear limit for all possible situations that could now serve as the benchmark for nonlinear stability studies in Keplerian accretion disks.

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Characteristics of surface-water flows in the ridge and slough landscape of Everglades National Park: implications for particulate transport

Over the last one hundred years, compartmentalization and water management activities have reduced water flow to the ridge and slough landscape of the Everglades. As a result, the once corrugated landscape has become topographically and vegetationally uniform. The focus of this study was to quantify variation in surface flow in the ridge and slough landscape and to relate flow conditions to particulate transport and deposition. Over the 2002–2003 and 2003–2004 wet seasons, surface velocities and particulate accumulation were measured in upper Shark River Slough in Everglades National Park. Landscape characteristics such as elevation, plant density and biomass also were examined to determine their impact on flow characteristics and material transport. The results of this study demonstrate that the release of water during the wet season not only increases water levels, but also increased flow speeds and particulate transport and availability. Further, flow speeds were positively and significantly correlated with water level thereby enhancing particulate transport in sloughs relative to ridges especially during peak flow periods. Our results also indicate that the distribution of biomass in the water column, including floating plants and periphyton, affects velocity magnitude and shape of vertical profiles, especially in the sloughs where Utricularia spp. and periphyton mats are more abundant. Plot clearing experiments suggest that the presence of surface periphyton and Utricularia exert greater control over flow characteristics than the identity (i.e., sawgrass or spike rush) or density of emergent macrophytes, two parameters frequently incorporated into models describing flow through vegetated canopies. Based on these results, we suggest that future modeling efforts must take the presence of floating biomass, such as Utricularia, and presence of periphyton into consideration when describing particulate transport.

Influence of hydrobiogeochemistry on transport of chromium through manganese -containing sediments: Experimental and modeling approaches

Chromium (Cr) is a metal of particular environmental concern, owing to its toxicity and widespread occurrence in groundwater, soil, and soil solution. A combination of hydrological, geochemical, and microbiological processes governs the subsurface migration of Cr. Little effort has been devoted to examining how these biogeochemical reactions combine with hydrologic processes influence Cr migration. This study has focused on the complex problem of predicting the Cr transport in laboratory column experiments. A 1-D reactive transport model was developed and evaluated against data obtained from laboratory column experiments. ^ A series of dynamic laboratory column experiments were conducted under abiotic and biotic conditions. Cr(III) was injected into columns packed with β-MnO 2-coated sand at different initial concentrations, variable flow rates, and at two different pore water pH (3.0 and 4.0). In biotic anaerobic column experiments Cr(VI) along with lactate was injected into columns packed with quartz sand or β-MnO2-coated sand and bacteria, Shewanella alga Simidu (BrY-MT). A mathematical model was developed which included advection-dispersion equations for the movement of Cr(III), Cr(VI), dissolved oxygen, lactate, and biomass. The model included first-order rate laws governing the adsorption of each Cr species and lactate. The equations for transport and adsorption were coupled with nonlinear equations for rate-limited oxidation-reduction reactions along with dual-monod kinetic equations. Kinetic batch experiments were conducted to determine the reduction of Cr(VI) by BrY-MT in three different substrates. Results of the column experiments with Cr(III)-containing influent solutions demonstrate that β-MnO2 effectively catalyzes the oxidation of Cr(III) to Cr(VI). For a given influent concentration and pore water velocity, oxidation rates are higher, and hence effluent concentrations of Cr(VI) are greater, at pH 4 relative to pH 3. Reduction of Cr(VI) by BrY-MT was rapid (within one hour) in columns packed with quartz sand, whereas Cr(VI) reduction by BrY-MT was delayed (57 hours) in presence of β-MnO 2-coated sand. BrY-MT grown in BHIB (brain heart infusion broth) reduced maximum amount of Cr(VI) to Cr(III) followed by TSB (tryptic soy broth) and M9 (minimum media). The comparisons of data and model results from the column experiments show that the depths associated with Cr(III) oxidation and transport within sediments of shallow aquatic systems can strongly influence trends in surface water quality. The results of this study suggests that carefully performed, laboratory column experiments is a useful tool in determining the biotransformation of redox-sensitive metals even in the presence of strong oxidant, like β-MnO2. ^

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Improvement Of The Physical Performance Is Associated With Activation Of No/pgc-1α/mttfa Signaling Pathway And Increased Protein Expressions Of Electron Transport Chain In Gastrocnemius Muscle From Rats Supplemented With L-arginine.

To examine the influence of l-arginine supplementation in combination with physical training on mitochondrial biomarkers from gastrocnemius muscle and its relationship with physical performance. Male Wistar rats were divided into four groups: control sedentary (SD), sedentary supplemented with l-arginine (SDLA), trained (TR) and trained supplemented with l-arginine (TRLA). Supplementation of l-arginine was administered by gavage (62.5mg/ml/day/rat). Physical training consisted of 60min/day, 5days/week, 0% grade, speed of 1.2km/h. The study lasted 8weeks. Skeletal muscle mitochondrial enriched fraction as well as cytoplasmic fractions were obtained for Western blotting and biochemical analyses. Protein expressions of transcriptor coactivator (PGC-1α), transcriptor factors (mtTFA), ATP synthase subunit c, cytochrome oxidase (COXIV), constitutive nitric oxide synthases (eNOS and nNOS), Cu/Zn-superoxide dismutase (SOD) and manganese-SOD (Mn-SOD) were evaluated. We also assessed in plasma: lipid profile, glycemia and malondialdehyde (MDA) levels. The nitrite/nitrate (NOx(-)) levels were measured in both plasma and cytosol fraction of the gastrocnemius muscle. 8-week l-arginine supplementation associated with physical training was effective in promoting greater tolerance to exercise that was accompanied by up-regulation of the protein expressions of mtTFA, PGC-1α, ATP synthase subunit c, COXIV, Cu/Zn-SOD and Mn-SOD. The upstream pathway was associated with improvement of NO bioavailability, but not in NO production since no changes in nNOS or eNOS protein expressions were observed. This combination would be an alternative approach for preventing cardiometabolic diseases given that in overt diseases a profound impairment in the physical performance of the patients is observed.

Characterisation Of The Membrane Transport Of Pilocarpine In Cell Suspension Cultures Of Pilocarpus Microphyllus.

Pilocarpine is an alkaloid obtained from the leaves of Pilocarpus genus, with important pharmaceutical applications. Previous reports have investigated the production of pilocarpine by Pilocarpus microphyllus cell cultures and tried to establish the alkaloid biosynthetic route. However, the site of pilocarpine accumulation inside of the cell and its exchange to the medium culture is still unknown. Therefore, the aim of this study was to determine the intracellular accumulation of pilocarpine and characterise its transport across membranes in cell suspension cultures of P. microphyllus. Histochemical analysis and toxicity assays indicated that pilocarpine is most likely stored in the vacuoles probably to avoid cell toxicity. Assays with exogenous pilocarpine supplementation to the culture medium showed that the alkaloid is promptly uptaken but it is rapidly metabolised. Treatment with specific ABC protein transporter inhibitors and substances that disturb the activity of secondary active transporters suppressed pilocarpine uptake and release suggesting that both proteins may participate in the traffic of pilocarpine to inside and outside of the cells. As bafilomicin A1, a specific V-type ATPase inhibitor, had little effect and NH4Cl (induces membrane proton gradient dissipation) had moderate effect, while cyclosporin A and nifedipine (ABC proteins inhibitors) strongly inhibited the transport of pilocarpine, it is believed that ABC proteins play a major role in the alkaloid transport across membranes but it is not the exclusive one. Kinetic studies supported these results.