Non-darcy double-diffusive mixed convection from heated vertical and horizontal plates in saturated porous media

The non-darcy mixed convection flows from heated vertical and horizontal plates in saturated porous media have been considered using boundary layer approximations. The flows are considered to be driven by multiple buoyancy forces. The similarity solutions for both vertical and horizontal plates have been obtained. The governing equations have been solved numerically using a shooting method. The heat transfer, mass transfer and skin friction are reduced due to inertial forces. Also, they increase with the buoyancy parameter for aiding flow and decrease for the opposing flow. For aiding flow, the heat and mass transfer coefficients are found to approach asymptotically the forced or free convection values as the buoyancy parameter approaches zero or infinity.

Comparative study of different formulations and solution methods for two phase flow in saturated porous media

Six models (Simulators) are formulated and developed with all possible combinations of pressure and saturation of the phases as primary variables. A comparative study between six simulators with two numerical methods, conventional simultaneous and modified sequential methods are carried out. The results of the numerical models are compared with the laboratory experimental results to study the accuracy of the model especially in heterogeneous porous media. From the study it is observed that the simulator using pressure and saturation of the wetting fluid (PW, SW formulation) is the best among the models tested. Many simulators with nonwetting phase as one of the primary variables did not converge when used along with simultaneous method. Based on simulator 1 (PW, SW formulation), a comparison of different solution methods such as simultaneous method, modified sequential and adaptive solution modified sequential method are carried out on 4 test problems including heterogeneous and randomly heterogeneous problems. It is found that the modified sequential and adaptive solution modified sequential methods could save the memory by half and as also the CPU time required by these methods is very less when compared with that using simultaneous method. It is also found that the simulator with PNW and PW as the primary variable which had problem of convergence using the simultaneous method, converged using both the modified sequential method and also using adaptive solution modified sequential method. The present study indicates that pressure and saturation formulation along with adaptive solution modified sequential method is the best among the different simulators and methods tested.

Parameters Inversion of Fluid-Saturated Porous Media Based on Neural Networks

The multi-layers feedforward neural network is used for inversion of material constants of fluid-saturated porous media. The direct analysis of fluid-saturated porous media is carried out with the boundary element method. The dynamic displacement responses obtained from direct analysis for prescribed material parameters constitute the sample sets training neural network. By virtue of the effective L-M training algorithm and the Tikhonov regularization method as well as the GCV method for an appropriate selection of regularization parameter, the inverse mapping from dynamic displacement responses to material constants is performed. Numerical examples demonstrate the validity of the neural network method.

A coupled hydro-mechanical analysis for prediction of hydraulic fracture propagation in saturated porous media using EFG mesh-less method

The details of the Element Free Galerkin (EFG) method are presented with the method being applied to a study on hydraulic fracturing initiation and propagation process in a saturated porous medium using coupled hydro-mechanical numerical modelling. In this EFG method, interpolation (approximation) is based on nodes without using elements and hence an arbitrary discrete fracture path can be modelled.The numerical approach is based upon solving two governing partial differential equations of equilibrium and continuity of pore water simultaneously. Displacement increment and pore water pressure increment are discretized using the same EFG shape functions. An incremental constrained Galerkin weak form is used to create the discrete system of equations and a fully implicit scheme is used for discretization in the time domain. Implementation of essential boundary conditions is based on the penalty method. In order to model discrete fractures, the so-called diffraction method is used.Examples are presented and the results are compared to some closed-form solutions and FEM approximations in order to demonstrate the validity of the developed model and its capabilities. The model is able to take the anisotropy and inhomogeneity of the material into account. The applicability of the model is examined by simulating hydraulic fracture initiation and propagation process from a borehole by injection of fluid. The maximum tensile strength criterion and Mohr-Coulomb shear criterion are used for modelling tensile and shear fracture, respectively. The model successfully simulates the leak-off of fluid from the fracture into the surrounding material. The results indicate the importance of pore fluid pressure in the initiation and propagation pattern of fracture in saturated soils. © 2013 Elsevier Ltd.

Quantifying mircoorganism/organic matter interactions in saturated porous media.

Proteins and humic acids are common constituents of waste water. Latex colloids (colloids) acted as surrogates for microorganisms in multiple pulse dynamic column experiments (MPEs) that permitted colloid mobility to be quantified before and after the injection of either BSA (a protein), or Suwannee River humic acid (SRHA).
At low OM coverage colloid breakthrough curves demonstrated both BSA and SRHA reduced colloid deposition rates, but did not affect colloid irreversible deposition mechanisms. By contrast, high levels of SRHA surface coverage not only further reduced the matrix’s ability to attenuate colloids, but also resulted in reversible adsorption of a significant fraction of colloids deposited. Modelling of colloid responses using random sequential adsorption modelling suggested that 1 microgram of SRHA had the same effect as the deposition of 5.90±0.14 x109 colloids; the model suggested that adsorption of the same mass of BSA was equivalent to the deposition of between 7.1x108 and 2.3x109 colloids.
Colloid responses in MPEs where BSA coverage of colloid deposition sites approached saturation demonstrated the sand matrix remained capable of adsorbing colloids. However, in contrast to responses observed in MPEs at low surface coverage, continued colloid injection showed that the sand’s attenuation capacity increased with time, i.e. colloid concentrations declined as more were deposited (filter ripening).
Importance: Study results highlight the contrasting responses that may arise due to the interactions between colloids and OM in porous media. Results not only underscore that colloids can interact differently with various forms of deposited OM, but also that a single type of OM may generate dramatically different responses depending on the degree of surface coverage. The MPE method provides a means of quantifying the influence of OM on microorganism mobility in porous media such as filter beds, which may be used for either drinking water treatment or waste water treatment. In the wider environment study findings have potential to allow more confident predictions of the mobility of sewage derived pathogens discharging to groundwater.

BEM modeling of saturated porous media susceptible to damage

This paper deals with the numerical analysis of saturated porous media, taking into account the damage phenomena on the solid skeleton. The porous media is taken into poro-elastic framework, in full-saturated condition, based on Biot's Theory. A scalar damage model is assumed for this analysis. An implicit boundary element method (BEM) formulation, based on time-independent fundamental solutions, is developed and implemented to couple the fluid flow and two-dimensional elastostatic problems. The integration over boundary elements is evaluated using a numerical Gauss procedure. A semi-analytical scheme for the case of triangular domain cells is followed to carry out the relevant domain integrals. The non-linear problem is solved by a Newton-Raphson procedure. Numerical examples are presented, in order to validate the implemented formulation and to illustrate its efficacy. (C) 2011 Elsevier Ltd. All rights reserved.

Stability, viscous dissipation and local thermal non-equilibrium in fluid saturated porous media

In this thesis, the field of study related to the stability analysis of fluid saturated porous media is investigated. In particular the contribution of the viscous heating to the onset of convective instability in the flow through ducts is analysed. In order to evaluate the contribution of the viscous dissipation, different geometries, different models describing the balance equations and different boundary conditions are used. Moreover, the local thermal non-equilibrium model is used to study the evolution of the temperature differences between the fluid and the solid matrix in a thermal boundary layer problem. On studying the onset of instability, different techniques for eigenvalue problems has been used. Analytical solutions, asymptotic analyses and numerical solutions by means of original and commercial codes are carried out.

Three-dimensional simulation of the rolling of a saturated fibro-porous media

This paper describes recent advances made in computational modelling of the sugar cane liquid extraction process. The saturated fibro-porous material is rolled between circumferentially grooved rolls, which enhance frictional grip and provide a low-resistance path for liquid flow during the extraction process. Previously reported two-dimensional (2D) computational models, account for the large deformation of the porous material by solving the fully coupled governing fibre stress and fluid-flow equations using finite element techniques. While the 2D simulations provide much insight into the overarching cause-effect relationships, predictions of mechanical quantities such as roll separating force and particularly torque as a function of roll speed and degree of compression are not satisfactory for industrial use. It is considered that the unsatisfactory response in roll torque prediction may be due to the stress levels that exist between the groove tips and roots which have been largely neglected in the geometrically simplified 2D model. This paper gives results for both two- and three-dimensional finite element models and highlights their strengths and weaknesses in predicting key milling parameters. (c) 2005 Elsevier B.V. All rights reserved.

Free convection in a triangular enclosure with fluid-saturated porous medium and internal heat generation

Unsteady natural convection inside a triangular cavity has been studied in this study. The cavity is filled with a saturated porous medium with non-isothermal left inclined wall while the bottom surface is isothermally heated and the right inclined surface is isothermally cold. An internal heat generation is also considered which is dependent of the fluid temperature. The governing equations are solved numerically by finite element method. The Prandtl number of the fluid is considered as 0.7 (air) while the aspect ratio and the Rayleigh number are considered as 0.5 and 105 respectively. The effect of the porosity of the medium and heat generation on the fluid flow and heat transfer have been presented as a form of streamlines and isotherms. The rate of heat transfer through three surfaces of the enclosure is also presented.

Unsteady Natural Convection Flow over a Heated Cylinder Buried in a Fluid Saturated Porous Medium

Transient natural convection flow on a heated cylinder buried in a semi-infinite liquid-saturated porous medium has been studied. The unsteadiness in the problem arises due to the cylinder which is heated (cooled) suddenly and then maintained at that temperature. The coupled partial differential equations governing the flow and heat transfer are cast into stream function-temperature formulation, and the solutions are obtained from the initial time to the time when steady state is reached. The heat transfer is found to change significantly with increasing time in a small time interval immediately after the start of the impulsive change, and steady state is reached after some time. The average Nusselt number is found to increase with Rayleigh number When the surface of the cylinder is suddenly cooled, there is a change in the direction of the heat transfer in a small time interval immediately after the start of the impulsive change in the surface temperature;however when the surface temperature is suddenly increased, no such phenomenon is observed.

Modeling the co-transport of viruses and colloids in unsaturated porous media

A mathematical model is developed to simulate the co-transport of viruses and colloids in unsaturated porous media under steady-state flow conditions. The virus attachment to the mobile and immobile colloids is described using a linear reversible kinetic model. Colloid transport is assumed to be decoupled from virus transport; that is, we assume that colloids are not affected by the presence of attached viruses on their surface. The governing equations,are solved numerically using an alternating three-step operator splitting approach. The model is verified by fitting three sets of experimental data published in the literature: (1) Syngouna and Chrysikopoulos (2013) and (2) Walshe et al. (2010), both on the co-transport of viruses and clay colloids under saturated conditions, and (3) Syngouna and Cluysikopoulos (2015) for the co-transport of viruses and clay colloids under unsaturated conditions. We found a good agreement between observed and fitted breakthrough curves (BTCs) under both saturated and unsaturated conditions. Then, the developed model was used to simulate the co-transport of viruses and colloids in porous media under unsaturated conditions, with the aim of understanding the relative importance of various processes on the co-transport of viruses and colloids in unsaturated porous media. The virus retention in porous media in the presence of colloids is greater during unsaturated conditions as compared to the saturated conditions due to: (1) virus attachment to the air-water interface (AWI), and (2) co-deposition of colloids with attached viruses on its surface to the AWL A sensitivity analysis of the model to various parameters showed that the virus attachment to AWI is the most sensitive parameter affecting the BTCs of both free viruses and total mobile viruses and has a significant effect on all parts of the BTC. The free and the total mobile viruses BTCs are mainly influenced by parameters describing virus attachment to the AIM, virus interaction with mobile and immobile colloids, virus attachment to solid-water interface (SWI), and colloid interaction with SWI and AWL The virus BTC is relatively insensitive to parameters describing the maximum adsorption capacity of the AWI for colloids, inlet colloid concentration, virus detachment rate coefficient from the SW!, maximum adsorption capacity of the AWI for viruses and inlet virus concentration. (C) 2015 Elsevier B.V. All rights reserved.

Permeability of Fractal Porous Media by Monte Carlo Simulations

The permeability of the fractal porous media is simulated by Monte Carlo technique in this work. Based oil the fractal character of pore size distribution in porous media, the probability models for pore diameter and for permeability are derived. Taking the bi-dispersed fractal porous media as examples, the permeability calculations are performed by the present Monte Carlo method. The results show that the present simulations present a good agreement compared with the existing fractal analytical solution in the general interested porosity range. The proposed simulation method may have the potential in prediction of other transport properties (such as thermal conductivity, dispersion conductivity and electrical conductivity) in fractal porous media, both saturated and unsaturated.

Can Pickering emulsion formation aid the removal of creosote DNAPL from porous media?

The purpose of this investigation was to examine the proposition that creosote, emplaced in an initially water saturated porous system, can be removed from the system through Pickering emulsion formation. Pickering emulsions are dispersions of two immiscible fluids in which coalescence of the dispersed phase droplets is hindered by the presence of colloidal particles adsorbed at the interface between the two immiscible fluid phases. Particle trapping is strongly favoured when the wetting properties of the particles are intermediate between strong water wetting and strong oil wetting. In this investigation the necessary chemical conditions for the formation of physically stable creosote-in-water emulsions protected against coalescence by bentonite particles were examined. It was established that physically stable emulsions could be formed through the judicious addition of small amounts of sodium chloride and the surfactant cetyl-trimethylammonium bromide. The stability of the emulsions was initially established by visual inspection. However, experimental determinations of emulsion stability were also undertaken by use of oscillatory rheology. Measurements of the elastic and viscous responses to shear indicated that physically stable emulsions were obtained when the viscoelastic systems showed a predominantly elastic response to shearing. Once the conditions were established for the formation of physically stable emulsions a "proof-of-concept" chromatographic experiment was carried out which showed that creosote could be successfully removed from a saturated model porous system. (C) 2007 Elsevier Ltd. All rights reserved.

Monitoring microbial sulfate reduction in porous media using multi-purpose electrodes

There is growing interest in the application of electrode-based measurements for monitoring microbial processes in the Earth using biogeophysical methods. In this study, reactive electrode measurements were combined to electrical geophysical measurements during microbial sulfate reduction occurring in a column of silica beads saturated with natural river water. Electrodic potential (EP), self potential (SP) and complex conductivity signals were recorded using a dual electrode design (Ag/AgCl metal as sensing/EP electrode, Ag/AgCl metal in KCl gel as reference/SP electrode). Open-circuit potentials, representing the tendency for electrochemical reactions to occur on the electrode surfaces, were recorded between sensing/EP electrode and reference/SP electrode and showed significant spatiotemporal variability associated with microbial activity. The dual electrode design isolates the microbial driven sulfide reactions to the sensing electrode and permits removal of any SP signal from the EP measurement. Based on the known sensitivity of a Ag electrode to dissolved sulfide, we interpret EP signals exceeding 550 mV recorded in this experiment in terms of bisulfide (HS-) concentration near multiple sensing electrodes. Complex conductivity measurements capture an imaginary conductivity (s?) signal interpreted as the response of microbial growth and biomass formation in the column. Our results suggest that the implementation of multipurpose electrodes, combining reactive measurements with electrical geophysical measurements, could improve efforts to monitor microbial processes in the Earth using electrodes.