Modelling Biogrout: A New Ground Improvement Method Based on Microbial-Induced Carbonate Precipitation

Biogrout is a new soil reinforcement method based on microbial-induced carbonate precipitation. Bacteria are placed and reactants are flushed through the soil, resulting in calcium carbonate precipitation, causing an increase in strength and stiffness of the soil. Due to this precipitation, the porosity of the soil decreases. The decreasing porosity influences the permeability and therefore the flow. To analyse the Biogrout process, a model was created that describes the process. The model contains the concentrations of the dissolved species that are present in the biochemical reaction. These concentrations can be solved from a advection?Cdispersion?Creaction equation with a variable porosity. Other model equations involve the bacteria, the solid calcium carbonate concentration, the (decreasing) porosity, the flow and the density of the fluid. The density of the fluid changes due to the biochemical reactions, which results in density driven flow. The partial differential equations are solved by the Standard Galerkin finite-element method. Simulations are done for some 1D and 2D configurations. A 1D configuration can be used to model a column experiment and a 2D configuration may correspond to a sheet or a cross section of a 3D configuration.     

Macroscale Properties of Porous Media from a Network Model of Biofilm Processes

We demonstrate how a network model can predict porosity and permeability changes in a porous medium as a result of biofilm buildup in the pore spaces. A biofilm consists of bacteria and extracellular polymeric substances (EPS) bonded together and attached to a surface. In this case, the surface consists of the walls of the porous medium, which we model as a random network of pipes.Our model contains five species. Four of these are bacteria and EPS in both fluid and adsorbed phases. The fifth species is nutrient, which we assume to reside in the fluid phase only. Bacteria and EPS transfer between the adsorbed and fluid phases through adsorption and erosion or sloughing. The adsorbed species influence the effective radii of the pipes in the network, which affect the porosity and permeability.We develop a technique for integrating the coupled system of ordinary and partial differential equations that govern transport of these species in the network. We examine ensemble averages of simulations using different arrays of pipe radii having identical statistics. These averages show how different rate parameters in the biofilm transport processes affect the concentration and permeability profiles.     

Suspension-based rheological modeling of crystallizing polymer melts

The applicability of suspension models to polymer crystallization is discussed. Although direct numerical simulations of flowing particle-filled melts are useful for gaining understanding about the rheological phenomena involved, they are computationally expensive. A more coarse-grained suspension model, which can relate the parameters in a constitutive equation for the two-phase material to morphological features, such as the volume fractions of differently shaped crystallites and the rheological properties of both phases, will be more practical in numerical polymer processing simulations. General issues, concerning the modeling of linear and nonlinear viscoelastic phenomena induced by rigid and deformable particles, are discussed. A phenomenological extension of linear viscoelastic suspension models into the nonlinear regime is proposed. A number of linear viscoelastic models for deformable particles are discussed, focusing on their possibilities in the context of polymer crystallization. The predictions of the most suitable model are compared to direct numerical simulation results and experimental data.     

Using Microfluidic Set-Up to Determine the Adsorption Rate of Sporosarcina pasteurii Bacteria on Sandstone

Transport in Porous Media - Microbial-induced carbonate precipitation (MICP) in porous media is a two-step procedure: First, the suspension of bacteria is injected and some of the bacteria get...     

Modeling and control of a nonlinear half-vehicle suspension system: a hybrid fuzzy logic approach

Modeling and control of vehicle suspension system are high noteworthy from safety to comfort. In this paper, an analytical nonlinear half-vehicle model which is included quadratic tire stiffness, cubic suspension stiffness, and coulomb friction is derived based on fundamental physics. A hybrid fuzzy logic approach which combines fuzzy logic and PID controllers is designed for reducing the vibration levels of passenger seat and vehicle body. Performances of designed controllers have been evaluated by numerical simulations. Comparisons with classical PID control, Fuzzy Logic Control (FLC) and Hybrid Fuzzy-PID control (HFPID) have also been provided. Results of numerical simulations are evaluated in terms of time histories of displacement and acceleration responses and ride index comparison. A good performance for the Hybrid Fuzzy-PID controller with coupled rules (HFPIDCR) is achieved in simulation studies despite the nonlinearities.     

Effect of the stiffness and inertia of a rib reinforcement on localized bending waves in semi-infinite strips

This paper reports an analytical and numerical study of localized bending waves in a thin elastic isotropic semi-infinite strip with a rib reinforment. Such waves can be considered as spatially non-uniform bending perturbations localized near free edges, similar to Rayleigh waves decaying exponentially with the distance. From the analysis of localized bending waves in thin elastic structures it may be possible to infer the presence of imperfections, cracks or inclusions, in the structure. In this study, a semi-infinite strip, simply supported on two edges and reinforced on the free one, is considered. A general solution is presented and the conditions under which localized bending waves exist are derived. The mathematical model, accompanied by numerical simulations, reveals that the presence of a reinforcement rib can suppress localized bending waves. This effect is mainly due to the stiffness coming from the rib, rather than from its added inertia terms.     

端承粘弹性桩纵向振动的轴对称解析解

将桩和土体视为粘弹性介质,在频率域求得了桩轴对称纵向振动时动力响应的解析解,进而得到桩头复刚度的解析表达式.通过数值计算,给出了桩头动刚度因子和等效阻尼随激励频率的响应,考察了土层粘性、桩长径比、桩土模量比、桩泊松比和粘性等参数对桩头刚度因子和等效阻尼的影响.研究表明:由于考虑了桩的径向变形效应以及土层对桩的径向力作用,轴对称解析解的桩头动刚度因子和等效阻尼与经典Euler杆模型桩的桩头动刚度因子和阻尼有较大区别.因此,经典Euler杆模型桩的适用范围具有一定的局限性,在若干情形下应采用基于三维粘弹性的桩纵向振动轴对称模型进行分析.     

Diffusive Effects on Recovery of Light Oil by Medium Temperature Oxidation

Volatile oil recovery by means of air injection is studied as a method to improve recovery from low permeable reservoirs. We consider the case in which the oil is directly combusted into small products, for which we use the term medium temperature oil combustion. The two-phase model considers evaporation, condensation and reaction with oxygen. In the absence of thermal, molecular and capillary diffusion, the relevant transport equations can be solved analytically. The solution consists of three waves, i.e., a thermal wave, a medium temperature oxidation (MTO) wave and a saturation wave separated by constant state regions. A striking feature is that evaporation occurs upstream of the combustion reaction in the MTO wave. The purpose of this paper is to show the effect of diffusion mechanisms on the MTO process. We used a finite element package (COMSOL) to obtain a numerical solution; the package uses fifth-order Lagrangian base functions, combined with a central difference scheme. This makes it possible to model situations at realistic diffusion coefficients. The qualitative behavior of the numerical solution is similar to the analytical solution. Molecular diffusion lowers the temperature of the MTO wave, but creates a small peak near the vaporization region. The effect of thermal diffusion smoothes the thermal wave and widens the MTO region. Capillary diffusion increases the temperature in the upstream part of the MTO region and decreases the efficiency of oil recovery. At increasing capillary diffusion the recovery by gas displacement gradually becomes higher, leaving less oil to be recovered by combustion. Consequently, the analytical solution with no diffusion and numerical solutions at a high capillary diffusion coefficient become different. Therefore high numerical diffusion, significant in numerical simulations especially in coarse gridded simulations, may conceal the importance of combustion in recovering oil.     

The reinforcement and defect interaction of two-dimensional lattice materials with imperfections

The defect interaction and reinforcement of imperfect two-dimensional lattice materials are studied by theoretical investigations and finite element (FE) simulations. An analytical model is proposed to predict the interaction of two defects in lattice materials based on a single defect model. An interaction coefficient is introduced to characterize the degree of interaction. The effects of defect type and defect distance on interaction coefficients are studied. The critical interaction distance of defects, beyond which the interaction of two defects can be neglected, is derived. FE calculations are performed to validate the theoretical model. The simulated results indicate that increasing the number of defects can reduce the stress concentration rather than weakening the strength of the residual parts in certain circumstances. Subsequently, several reinforcement methods are proposed to reduce the stress concentration in the triangular and Kagome lattice for the single-bar-missing defect and single-joint-missing defect. An analytical model is developed for the reinforced lattices, and the predicted stress concentration factors are in good agreement with those of FE simulations. By theoretical studies and FE simulations, optimal reinforcement methods are derived for the triangular and Kagome lattice under planar loading conditions.     

Theoretical analysis and numerical simulation for the spill procedure of liquid fuel of fuel air explosive with shell

Addressing the problem of spill of liquid fuel of fuel air explosive (FAE) with shell, this paper deduces variation relationship of fracture radius and fracture velocity based on continuity and motion equations of shell using the Taylor fracture criteria of cylindrical shells, and analyzes the correctness of the approximating process in detail. Furthermore, analytical solutions and numerical simulations for the liquid fuel spill process were obtained using Adomian's decomposition method (ADM) and the Runge-Kutta numerical method. Results show that both the numerical solution and the analytical solution based on ADM are closer to the result from experiments and computations from commercial dynamics numerical software.     

Hillslope soil erosion process model for natural rainfall events

Based on the momentum theorem and water-balance principle, this paper has derived the basic equations of slope runoff. Soil erosion by raindrop splash and runoff are discussed, and hillslope soil erosion processes model was established. Finally, the model was solved by numerical solution considering the common solution-determining conditions, from which not only runoff and soil erosion yield but also their processes can be obtained. The model is validated by using observation data of Soil Conservation Ecological Science & Technology Demonstration Park of Jiangxi Province. Results show that the maximum relative error between the computation and observation values of the slop runoff processes is 10％ while the maximum relative error between the computation and observation values of hillslope soil erosion processes is 9.2％. This indicates that the model is conceptually realistic reasonable and offers a feasible approach for research on soil erosion process.Key words overland flow﹒raindrop splash﹒runoff erosion﹒process﹒numerical solution      

Forward and Inverse Bio-Geochemical Modeling of Microbially Induced Calcite Precipitation in Half-Meter Column Experiments

Microbially induced calcite precipitation (MICP) offers an alternative solution to a wide range of civil engineering problems. Laboratory tests have shown that MICP can immobilize trace metals and radionuclides through co-precipitation with calcium carbonate. MICP has also been shown to improve the undrained shear response of soils and offers potential benefits over current ground improvement techniques that may pose environmental risks and suffer from low “certainty of execution.” Our objective is to identify an effective means of achieving uniform distribution of precipitate in a one-dimensional porous medium. Our approach involves column experiments and numerical modeling of MICP in both forward and inverse senses, using a simplified reaction network, with the bacterial strain Sporoscarcina pasteurii. It was found that the stop-flow injection of a urea- and calcium-rich solution produces a more uniform calcite distribution as compared to a continuous injection method, even when both methods involve flow in opposite direction to that used for bacterial cell emplacement. Inverse modeling was conducted by coupling the reactive transport code TOUGHREACT to UCODE for estimating chemical reaction rate parameters with a good match to the experimental data. It was found, however, that the choice of parameters and data was not sufficient to determine a unique solution, and our findings suggest that additional time and space-varying analytical data of aqueous species would improve the accuracy of numerical modeling of MICP.     

Effects of Surface Waves and Marine Soil Parameters on Seabed Stability

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｍｏｓｔｉｍｐｏｒｔａｎｔｆｏｒｃｅｓｅｘｅｒｔｉｎｇｏｎｃｏａｓｔａｌｓｔｒｕｃｔｕｒｅｓ,ｓｕｃｈａｓｏｆｆｓｈｏｒｅｄｒｉｌｌｉｎｇｒｉｇｓ,ｏｆｆｓｈｏｒｅｏｉｌｓｔｏｒａｇｅｔａｎｋｓ,ｂｕｒｉｅｄｐｉｐｅｌｉｎｅｓ,ｂｒｅａｋｗａｔｅｒｓ,ｅｔｃ .ａｒｅｄｕｅｔｏｗａｔｅｒｓｕｒｆａｃｅｗａｖｅｓ.Ｔｈｅｗａｔｅｒｗａｖｅｓｐｒｏｐａｇａｔｉｎｇｏｖｅｒｔｈｅｏｃｅａｎｃａｎｃｒｅａｔｅｔｒｅｍｅｎｄｏｕｓｄｙｎａｍｉｃｐｒｅｓｓｕｒｅ .Ｔｈｅｐｒｅｓｓｕｒｅ…     

Analytical solution for a reinforcement layer bonded to an elliptic hole under a remote uniform load

A general solution to a reinforced elliptic hole embedded in an infinite matrix subjected to a remote uniform load is provided in this paper. Investigations on the present elasticity problem are rather tedious due to the presence of material inhomogeneities and complex geometric configurations. Based on the technique of conformal mapping and the method of analytical continuation in conjunction with the alternating technique, the general expressions of the displacement and stresses in a reinforcement layer and the matrix are derived explicitly in a series form. Some numerical results are provided to investigate the effects of the material combinations and geometric configurations on the interfacial stresses. The results show that there exists an optimum design of a reinforcement layer such that both the magnitude of stress concentration and the interfacial stresses could be fairly reduced.     

A two-phase elastoplastic model for unidirectionally-reinforced materials

A complete analytical formulation for the elastoplastic behaviour of a composite material comprising one single array of reinforcing inclusions perfectly bonded to the matrix is developed in this paper. Fundamental relationships establish the link between the total stress and strain variables, and those pertaining to the individual constituents (matrix and reinforcement) regarded as superposed continuous phases. Assuming that each constituent behaves as an elastic perfectly plastic material, the constitutive equations governing the evolution of the reinforced material as a whole are derived. They reveal a hardening phenomenon arising from the non-compatibility between matrix and reinforcement plastic strains. It is shown in particular that the obtained constitutive law falls within the formalism of generalized standard plasticity: the reinforcement residual stress plays the role of a hardening parameter which controls the evolution of the yield surface, while the associated kinematic variable is the plastic strain discrepancy between matrix and reinforcement phases.Owing to its inherent simplicity, the model is easily amenable to a numerical treatment for structural analysis. It is shown in particular how the classical iterative algorithm can be modified accordingly, and an illustrative application is finally presented in the field of civil engineering.     

Theoretical statistical solution and numerical simulation of heterogeneous brittle materials

The analytical stress-strain relation with heterogeneous parameters is derived for theheterogeneous brittle materials under a uniaxial extensional load, in which the distributions of theelastic modulus and the failure strength are assumed to be statistically independent. This theoreticalsolution gives an approximate estimate of the equivalent stress-strain relations for 3-D heterogeneousmaterials. In one-dimensional cases it may provide comparatively accurate results. The theoreticalsolution can help us to explain how the heterogeneity influences the mechanical behaviors, Further, anumerical approach is developed to model the non-linear behavior of three-dimensional heterogeneousbrittle materials. The lattice approach and statistical techniques are applied to simulate the initialheterogeneity of heterogeneous materials. The load increment in each loading stage is adaptivelydetermined so that the better approximation of the failure process can be realized. When the maximumtensile principal strain exceeds the failure strain, the elements are considered to be broken, which canbe carried out by replacing its Young‘s modulus with a very small value. A 3-D heterogeneous brittlematerial specimen is simulated during a full failure process. The numerical results are in good agreementwith the analytical solutions and experimental data.     

A Similarity Solution for Oil Lens Redistribution Including Capillary Forces and Oil Entrapment

Redistribution of a LNAPL lens (oil) at the phreatic surface is described using a multi-phase flow model, with emphasis on the effect of oil entrapment by water. The flow process is analyzed under the assumption that the vertical capillary and gravitational forces balance. Vertical integration leads to explicit functions which approximate the relations between the free oil volume per unit lateral area and the vertically averaged oil relative permeability on the one hand and the vertical position of the interface between zones with either two or three phases on the other hand. A linear relation between the trapped and free oil volume per unit lateral area approximates the vertically integrated nonlinear expression for the trapped oil saturation. The resulting differential equation admits a similarity solution describing the lateral spreading of free oil and the amount and location of trapped oil. Comparison with illustrative numerical computations, which are based on the nonreduced flow model in a two-dimensional planar or axisymmetric domain, shows that the analytical solution provides a good approximation of the free oil distribution at all later times.     

横观各向同性弹性半空间地基上四边自由各向异性矩形薄板弯曲解析解

选用更具广泛性的横观各向同性弹性半空间地基模型,来分析四边自由各向异性矩形地基板的弯曲解析解．将异性薄板的弯曲控制方程,与基于横观各向同性弹性半空间地基位移解建立的板与地基变形协调方程相结合,先按对称性分解,然后用三角级数法,得出横观各向同性弹性半空间地基上四边自由各向异性矩形薄板的弯曲解析解,包括地基反力、板的挠度及内力的解析表达式．该解析解克服了数值法的弊端,取消了对地基反力的假设,板的内力及地基反力求解更切实际．算例结果与文献结果吻合良好,证明本文方法的可行性．     

Simulation of crack spacing using a reinforced concrete model with an internal length parameter

Summary A gradient-enhanced smeared crack model and bond-slip interface elements are utilized in finite element simulations of reinforced concrete. The crack model is rooted in an enhanced plasticity theory. It uses the Rankine failure surface dependent on an equivalent inelastic strain measure as well as on its Laplacian. As a result, finitely sized fracture process zones and realistic crack spacings are obtained. A reinforced concrete bar in uniaxial tension is analyzed to demonstrate the regularizing influence of the internal length parameter in the model and to evaluate the influence of the model parameters on the energy dissipation in multiple cracks. A comparison of numerical simulations with experimental results for a beam without shear reinforcement in four-point bending concludes the analysis. Received 4 November 1997; accepted for publication 23 April 1998