Experimental and numerical studies on protection of buried pipelines and underground utilities using geocells

This paper presents the results of the laboratory model tests and the numerical studies conducted on small diameter PVC pipes, buried in geocell reinforced sand beds. The aim of the study was to evaluate the suitability of the geocell reinforcement in protecting the underground utilities and buried pipelines. In addition to geocells, the efficacy of only geogrid and geocell with additional basal geogrid cases were also studied. A PVC (Poly Vinyl Chloride) pipe with external diameter 75 mm and thickness 1.4 mm was used in the experiments. The vehicle tire contact pressure was simulated by applying the pressure on the top of the bed with the help of a steel plate. Results suggest that the use of geocells with additional basal geogrid considerably reduces the deformation of the pipe as compared to other types of reinforcements. Further, the depth of placement of pipe was also varied between 1B to 2B (B is the width of loading plate) below the plate in the presence of geocell with additional basal geogrid. More than 50% reduction in the pressure and more than 40% reduction in the strain values were observed in the presence of reinforcements at different depths as compared to the unreinforced beds. Conversely, the performance of the subgrade soil was also found to be marginally influenced by the position of the pipe, even in the presence of the relatively stiff reinforcement system. Further, experimental results were validated with 3-dimensional numerical studies using FLAC(3D) (Fast Lagrangian Analysis of Continua in 3D). A good agreement in the measured pipe stain values were observed between the experimental and numerical studies. Numerical studies revealed that the geocells distribute the stresses in the lateral direction and thus reduce the pressure on the pipe. In addition, the results of the 1-g model tests were scaled up to the prototype case of the shallow buried pipeline below the pavement using the appropriate scaling laws. (C) 2015 Elsevier Ltd. All rights reserved.

Three-dimensional numerical analysis of geocell-reinforced soft clay beds by considering the actual geometry of geocell pockets

Due to its complex honeycomb structure, the numerical modeling of the geocell has always been a big challenge. Generally, the equivalent composite approach is used to model the geocells. In the equivalent composite approach, the geocellsoil composite is treated as the soil layer with improved strength and stiffness values. Though this approach is very simple, it is unrealistic to model the geocells as the soil layer. This paper presents a more realistic approach of modeling the geocells in three-dimensional (3D) framework by considering the actual curvature of the geocell pocket. A square footing resting on geocell reinforced soft clay bed was modeled using the ``fast Lagrangian analysis of continua in 3D'' (FLAC(3D)) finite difference package. Three different material models, namely modified Cam-clay, Mohr-Coulomb, and linear elastic were used to simulate the behaviour of foundation soil, infill soil and the geocell, respectively. It was found that the geocells distribute the load laterally to the wider area below the footing as compared to the unreinforced case. More than 50% reduction in the stress was observed in the clay bed in the presence of geocells. In addition to geocells, two other cases, namely, only geogrid and geocell with additional basal geogrid cases were also simulated. The numerical model was systematically validated with the results of the physical model tests. Using the validated numerical model, parametric studies were conducted to evaluate the influence of various geocell properties on the performance of reinforced clay beds.

Numerical study of laminar, standing hydraulic jumps in a planar geometry

We solve the two-dimensional, planar Navier-Stokes equations to simulate a laminar, standing hydraulic jump using a Volume-of-Fluid method. The geometry downstream of the jump has been designed to be similar to experimental conditions by including a pit at the edge of the platform over which liquid film flows. We obtain jumps with and without separation. Increasing the inlet Froude number pushes the jump downstream and makes the slope of the jump weaker, consistent with experimental observations of circular jumps, and decreasing the Reynolds number brings the jump upstream while making it steeper. We study the effect of the length of the domain and that of a downstream obstacle on the structure and location of the jump. The transient flow which leads to a final steady jump is described for the first time to our knowledge. In the moderate Reynolds number regime, we obtain steady undular jumps with a separated bubble underneath the first few undulations. Interestingly, surface tension leads to shortening of wavelength of these undulations. We show that the undulations can be explained using the inviscid theory of Benjamin and Lighthill (Proc. R. Soc. London, Ser. A, 1954). We hope this new finding will motivate experimental verification.

TURBULENCE-TRANSPORT-CHEMISTRY INTERACTION IN STATISTICALLY PLANAR PREMIXED FLAMES AND IGNITION KERNELS IN NEAR ISOTROPIC TURBULENCE

Turbulence-transport-chemistry interaction plays a crucial role on the flame surface geometry, local and global reactionrates, and therefore, on the propagation and extinction characteristics of intensely turbulent, premixed flames encountered in LPP gas-turbine combustors. The aim of the present work is to understand these interaction effects on the flame surface annihilation and extinction of lean premixed flames, interacting with near isotropic turbulence. As an example case, lean premixed H-2-air mixture is considered so as to enable inclusion of detailed chemistry effects in Direct Numerical Simulations (DNS). The work is carried out in two phases namely, statistically planar flames and ignition kernel, both interacting with near isotropic turbulence, using the recently proposed Flame Particle Tracking (FPT) technique. Flame particles are surface points residing and commoving with an iso-scalar surface within a premixed flame. Tracking flame particles allows us to study the evolution of propagating surface locations uniquely identified with time. In this work, using DNS and FPT we study the flame speed, reaction rate and transport histories of such flame particles residing on iso-scalar surfaces. An analytical expression for the local displacement flame speed (SO is derived, and the contribution of transport and chemistry on the displacement flame speed is identified. An examination of the results of the planar case leads to a conclusion that the cause of variation in S-d may be attributed to the effects of turbulent transport and heat release rate. In the second phase of this work, the sustenance of an ignition kernel is examined in light of the S-curve. A newly proposed Damkohler number accounting for local turbulent transport and reaction rates is found to explain either the sustenance or otherwise propagation of flame kernels in near isotropic turbulence.

Systematic Analysis of Mycobacterial Acylation Reveals First Example of Acylation-mediated Regulation of Enzyme Activity of a Bacterial Phosphatase

Protein lysine acetylation is known to regulate multiple aspects of bacterial metabolism. However, its presence in mycobacterial signal transduction and virulence-associated proteins has not been studied. In this study, analysis of mycobacterial proteins from different cellular fractions indicated dynamic and widespread occurrence of lysine acetylation. Mycobacterium tuberculosis proteins regulating diverse physiological processes were then selected and expressed in the surrogate host Mycobacterium smegmatis. The purified proteins were analyzed for the presence of lysine acetylation, leading to the identification of 24 acetylated proteins. In addition, novel lysine succinylation and propionylation events were found to co-occur with acetylation on several proteins. Protein-tyrosine phosphatase B (PtpB), a secretory phosphatase that regulates phosphorylation of host proteins and plays a critical role in Mycobacterium infection, is modified by acetylation and succinylation at Lys-224. This residue is situated in a lid region that covers the enzyme's active site. Consequently, acetylation and succinylation negatively regulate the activity of PtpB.

As-You-Go Deployment of a 2-Connected Wireless Relay Network for Sensor-Sink Interconnection

A person walks along a line (which could be an idealisation of a forest trail, for example), placing relays as he walks, in order to create a multihop network for connecting a sensor at a point along the line to a sink at the start of the line. The potential placement points are equally spaced along the line, and at each such location the decision to place or not to place a relay is based on link quality measurements to the previously placed relays. The location of the sensor is unknown apriori, and is discovered as the deployment agent walks. In this paper, we extend our earlier work on this class of problems to include the objective of achieving a 2-connected multihop network. We propose a network cost objective that is additive over the deployed relays, and accounts for possible alternate routing over the multiple available paths. As in our earlier work, the problem is formulated as a Markov decision process. Placement algorithms are obtained for two source location models, which yield a discounted cost MDP and an average cost MDP. In each case we obtain structural results for an optimal policy, and perform a numerical study that provides insights into the advantages and disadvantages of multi-connectivity. We validate the results obtained from numerical study experimentally in a forest-like environment.

Weakly corrected numerical solutions to stochastically driven nonlinear dynamical systems

A method to weakly correct the solutions of stochastically driven nonlinear dynamical systems, herein numerically approximated through the Eule-Maruyama (EM) time-marching map, is proposed. An essential feature of the method is a change of measures that aims at rendering the EM-approximated solution measurable with respect to the filtration generated by an appropriately defined error process. Using Ito's formula and adopting a Monte Carlo (MC) setup, it is shown that the correction term may be additively applied to the realizations of the numerically integrated trajectories. Numerical evidence, presently gathered via applications of the proposed method to a few nonlinear mechanical oscillators and a semi-discrete form of a 1-D Burger's equation, lends credence to the remarkably improved numerical accuracy of the corrected solutions even with relatively large time step sizes. (C) 2015 Elsevier Inc. All rights reserved.

Numerical Study of Counterflow Diffusion Flame and Water Spray Interaction

This paper reports numerical investigation concerning the interaction of a laminar methane-air counterflow diffusion flame with monodisperse and polydisperse water spray. Commercial code ANSYS FLUENT with reduced chemistry has been used for investigation. Effects of strain rate, Sauter mean diameter (SMD), and droplet size distribution on the temperature along stagnation streamline have been studied. Flame extinction using polydisperse water spray has also been explored. Comparison of monodisperse and polydisperse droplet distribution on flame properties reveals suitability of polydisperse spray in flame temperature reduction beyond a particular SMD. This study also provides a numerical framework to study flame-spray interaction and extinction.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

Scaling analysis and numerical studies on water vapour adsorption in a columnar porous silica gel bed

This article presents a theoretical analysis of heat and mass transfer in a silica gel + water adsorption process using scaling principles. A two-dimensional columnar packed adsorber domain is chosen for the study, with side and bottom walls cooled and vapour inlet from the top. The adsorption process is initiated from the cold walls with a temperature jump of 15 K, whereas the water vapour supply is maintained at a constant inlet pressure of 1 kPa. The first part of the study is dedicated to deriving relevant scales for the adsorption process by an order of magnitude analysis of energy, continuity and momentum equations. In the latter part, the derived scales are compared with the outcome of numerical studies performed for various domain widths and aspect ratio of bed. A good correlation between scaling and simulation results is observed, thereby validating the scaling approach. (C) 2015 Elsevier Ltd. All rights reserved.

Numerical analysis of segmented-electrode Orbitraps

Simple geometries which are possible alternatives for the Orbitrap are studied in this paper. We have taken up for numerical investigation two segmented-electrode structures, ORB1 and ORB2, to mimic the electric field of the Orbitrap. In the ORB1, the inner spindle-like electrode and the outer barrel-like electrode of the Orbitrap have been replaced by 35 rings and 35 discs of fixed radii, respectively. In this structure two segmented end cap electrodes have been added. In this geometry, different potentials are applied to the different electrodes keeping top-bottom symmetry intact. In the second geometry, ORB2, the inner and outer electrodes of the Orbitrap were replaced by an approximate step structure which follows the profile of the Orbitrap electrodes. In the present study 45 steps have been used. In the ORB2, like the Orbitrap, the inner electrode is held at a negative potential and the outer electrode is at ground potential. For the purpose of comparing the performance of ORB1 and ORB2 with that of the Orbitrap, the following studies have been undertaken: (1) variation of electric potential, (2) computation of ion trajectories, (3) simulation of image currents. These studies have been carried out using both 2D and 3D Boundary Element Method (BEM), the 3D BEM was developed specifically for this study. It has been seen in these investigations that ORB1 and ORB2 have performance similar to that of the Orbitrap, with the performance of the ORB1 being seen to be marginally superior to that of the ORB2. It has been shown that with proper optimization, geometries containing far fewer electrodes can be used as mass analyzers. A novel technique of optimization of the electric field has been proposed with the objective of minimizing the dependence of axial frequency of ion motion on the initial position of an ion. The results on the optimization of 9 and 15 segmented-electrode traps having the same design as ORB1 show that it can provide accurate mass analysis. (C) 2015 Elsevier B.V. All rights reserved.

Numerical Investigations for Leakage and Windage Heating in Straight-Through Labyrinth Seals

The ability to quantify leakage flow and windage heating for labyrinth seals with honeycomb lands is critical in understanding gas turbine engine system performance and predicting its component life. Variety of labyrinth seal configurations (number of teeth, stepped or straight, honeycomb cell size) are in use in gas turbines, and for each configuration, there are many geometric factors that can impact a seal's leakage and windage characteristics. This paper describes the development of a numerical methodology aimed at studying the effect of honeycomb lands on leakage and windage heating. Specifically, a three-dimensional computational fluid dynamics (CFD) model is developed utilizing commercial finite volume-based software incorporating the renormalization group (RNG) k-epsilon turbulence model with modified Schmidt number. The modified turbulence model is benchmarked and fine-tuned based on several experiments. Using this model, a broad parametric study is conducted by varying honeycomb cell size, pressure ratio (PR), and radial clearance for a four-tooth straight-through labyrinth seal. The results show good agreement with available experimental data. They further indicate that larger honeycomb cells predict higher seal leakage and windage heating at tighter clearances compared to smaller honeycomb cells and smooth lands. However, at open seal clearances larger honeycomb cells have lower leakage compared to smaller honeycomb cells.

Numerical prediction of nozzle flow separation: Issue of turbulence modeling

Numerical simulation of separated flows in rocket nozzles is challenging because existing turbulence models are unable to predict it correctly. This paper addresses this issue with the Spalart-Allmaras and Shear Stress Transport (SST) eddy-viscosity models, which predict flow separation with moderate success. Their performances have been compared against experimental data for a conical and two contoured subscale nozzles. It is found that they fail to predict the separation location correctly, exhibiting sensitivity to the nozzle pressure ratio (NPR) and nozzle type. A careful assessment indicated how the model had to be tuned for better, consistent prediction. It is learnt that SST model's failure is caused by limiting of the shear stress inside boundary layer according to Bradshaw's assumption, and by over prediction of jet spreading rate. Accordingly, SST's coefficients were empirically modified to match the experimental wall pressure data. Results confirm that accurate RANS prediction of separation depends on the correct capture of the jet spreading rate, and that it is feasible over a wide range of NPRs by modified values of the diffusion coefficients in the turbulence model. (C) 2015 Elsevier Masson SAS. All rights reserved.

Benchmark analytical solutions from beams with shared eigenpair

Structures with governing equations having identical inertial terms but somewhat differing stiffness terms can be termed flexurally analogous. An example of such a structure includes an axially loaded non-uniform beam and an unloaded uniform beam, for which an exact solution exists. We find that there exist shared eigenpairs (frequency and mode shapes) for a particular mode between such structures. Non-uniform beams with uniform axial loads, gravity loaded beams and rotating beams are considered and shared eigenpairs with uniform beams are found. In general, the derived flexural stiffness functions (FSF's) for the non-uniform beams required for the existence of shared eigenpair have internal singularities, but some of the singularities can be removed by an appropriate selection of integration constants using the theory of limits. The derived functions yield an insight into the relationship between the axial load and flexural stiffness of axially loaded beam structures. The derived functions can serve as benchmark solutions for numerical methods. (C) 2016 Elsevier Ltd. All rights reserved.

Numerical study of driven flows of shear thinning viscoelastic fluids in rectangular cavities

In this work, we present a numerical study of flow of shear thinning viscoelastic fluids in rectangular lid driven cavities for a wide range of aspect ratios (depth to width ratio) varying from 1/16 to 4. In particular, the effect of elasticity, inertia, model parameters and polymer concentration on flow features in rectangular driven cavity has been studied for two shear thinning viscoelastic fluids, namely, Giesekus and linear PTT. We perform numerical simulations using the symmetric square root representation of the conformation tensor to stabilize the numerical scheme against the high Weissenberg number problem. The variation in flow structures associated with merging and splitting of elongated vortices in shallow cavities and coalescence of corner eddies to yield a second primary vortex in deep cavities with respect to the variation in flow parameters is discussed. We discuss the effect of the dominant eigenvalues and the corresponding eigenvectors on the location of the primary eddy in the cavity. We also demonstrate, by performing numerical simulations for shallow and deep cavities, that where the Deborah number (based on convective time scale) characterizes the elastic behaviour of the fluid in deep cavities, Weissenberg number (based on shear rate) should be used for shallow cavities. (C) 2016 Elsevier B.V. All rights reserved.