An efficient model for the dynamic behaviour of a single pile in viscoelastic soil

This paper presents a novel, three-dimensional, single-pile model, formulated in the wavenumber domain and adapted to account for boundary conditions using the superposition of loading cases. The pile is modelled as a column in axial vibration, and a Euler-Bernoulli beam in lateral vibration. The surrounding soil is treated as a viscoelastic continuum. The response of the pile is presented in terms of the stiffness and damping coefficients, and also the magnitude and phase of the pile-head frequency-response function. Comparison with existing models shows that excellent agreement is observed between this model, a boundary-element formulation, and an elastic-continuum-type formulation. This three-dimensional model has an accuracy equivalent to a 3D boundary-element model, and a runtime similar to a 2D plane-strain analytical model. Analysis of the response of the single pile illustrates a difference in axial and lateral vibration behaviour; the displacement along the pile is relatively invariant under axial loads, but in lateral vibration the pile exhibits localised deformations. This implies that a plane-strain assumption is valid for axial loadings and only at higher frequencies for lateral loadings. © 2013 Elsevier Ltd.

An Arbitrary Lagrangian Eulerian Method for Three-Phase Flows with Triple Junction Points

We present a moving mesh method suitable for solving two-dimensional and axisymmetric three-liquid flows with triple junction points. This method employs a body-fitted unstructured mesh where the interfaces between liquids are lines of the mesh system, and the triple junction points (if exist) are mesh nodes. To enhance the accuracy and the efficiency of the method, the mesh is constantly adapted to the evolution of the interfaces by refining and coarsening the mesh locally; dynamic boundary conditions on interfaces, in particular the triple points, are therefore incorporated naturally and accurately in a Finite- Element formulation. In order to allow pressure discontinuity across interfaces, double-values of pressure are necessary for interface nodes and triple-values of pressure on triple junction points. The resulting non-linear system of mass and momentum conservation is then solved by an Uzawa method, with the zero resultant condition on triple points reinforced at each time step. The method is used to investigate the rising of a liquid drop with an attached bubble in a lighter liquid.

Numerical simulation of post-flame oxidation of hydrocarbons in spark ignition engines

About 50-90 percent of the hydrocarbons that escape combustion during flame passage in spark-ignition engine operation are oxidized in the cylinder before leaving the system. The process involves the transport of unreacted fuel from cold walls towards the hotter burned gas regions and subsequent reaction. In order to understand controlling factors in the process, a transient one-dimensional reactive-diffusive model has been formulated for simulating the oxidation processes taking place in the reactive layer between hot burned gases and cold unreacted air/fuel mixture, with initial and boundary conditions provided by the emergence of hydrocarbons from the piston top land crevice. Energy and species conservation equations are solved for the entire process, using a detailed chemical kinetic mechanism for propane. Simulation results show that the post-flame oxidation process takes place within a reactive layer where intermediate hydrocarbon products are formed at temperatures above 1100-1200 K, followed by a carbon monoxide conversion region closer to the hot burned gases. Model results show that most of hydrocarbons leaving the crevice are completely oxidized inside the cylinder. The largest contribution of remaining hydrocarbons are those leaving the crevice at temperatures below 1400 K. The largest fraction of non-fuel (intermediate) hydrocarbons results from hydrocarbons leaving the crevice when core temperatures are around 1400 K Copyright © 1997 Society of Automotive Engineers, Inc.

Wave propagation under vertically excited surface foundation

Recently developed equipment allows measurement of the shear modulus of soil in situ as a function of level of strain. In these field experiments, the excitation is applied on the ground surface using large scale shakers, and the response of the soil deposit is recorded through embedded receivers. The focus of this paper is on the simulation of signals which would be recorded at the receiver locations in idealized conditions to provide guidelines on the interpretation of field measurements. Discrete and finite element methods are employed to model one dimensional and three dimensional geometries, respectively, under various lateral boundary conditions. When the first times of arrival are detected by receivers under the vertical impulse, they coincide with the arrival of the P wave, related to the constrained modulus of the material, regardless of lateral boundary conditions. If one considers, on the other hand, phase differences between the motions at two receivers the picture is far more complicated and one would obtain propagation velocities, function of frequency and depth, which do not correspond to either the constrained modulus or Young's modulus. It is thus necessary to apply some care when interpreting the data from field tests based on vertical steady state vibrations. The use of inverse analysis can be considered as a way of extracting the shear modulus of soil from the field test measurements. © 2008 ASCE.

Simulations of thermal upgrading methods for oil shale reservoirs

This paper analyzes reaction and thermal front development in porous reservoirs with reacting flows, such as those encountered in shale oil extraction. A set of dimensionless parameters and a 3D code are developed in order to investigate the important physical and chemical variables of such reservoirs when heated by in situ methods. This contribution builds on a 1D model developed for the precursor study to this work. Theory necessary for this study is presented, namely shale decomposition chemical mechanisms, governing equations for multiphase flow in porous media and necessary closure models. Plotting the ratio of the thermal wave speed to the fluid speed allows one to infer that the reaction wave front ends where this ratio is at a minimum. The reaction front follows the thermal front closely, thus allowing assumptions to be made about the extent of decomposition solely by looking at thermal wave progression. Furthermore, this sensitivity analysis showed that a certain minimum permeability is required in order to ensure the formation of a traveling thermal wave. It was found that by studying the non-dimensional governing parameters of the system one can ascribe characteristic values for these parameters for given initial and boundary conditions. This allows one to roughly predict the performance of a particular method on a particular reservoir given approximate values for initial and boundary conditions. Channelling and flow blockage due to carbon residue buildup impeded each method's performance. Blockage was found to be a result of imbalanced heating. Copyright 2012, Society of Petroleum Engineers.

Guidance of collective cell migration by substrate geometry.

Collective behavior refers to the emergence of complex migration patterns over scales larger than those of the individual elements constituting a system. It plays a pivotal role in biological systems in regulating various processes such as gastrulation, morphogenesis and tissue organization. Here, by combining experimental approaches and numerical modeling, we explore the role of cell density ('crowding'), strength of intercellular adhesion ('cohesion') and boundary conditions imposed by extracellular matrix (ECM) proteins ('constraints') in regulating the emergence of collective behavior within epithelial cell sheets. Our results show that the geometrical confinement of cells into well-defined circles induces a persistent, coordinated and synchronized rotation of cells that depends on cell density. The speed of such rotating large-scale movements slows down as the density increases. Furthermore, such collective rotation behavior depends on the size of the micropatterned circles: we observe a rotating motion of the overall cell population in the same direction for sizes of up to 200 μm. The rotating cells move as a solid body, with a uniform angular velocity. Interestingly, this upper limit leads to length scales that are similar to the natural correlation length observed for unconfined epithelial cell sheets. This behavior is strongly altered in cells that present a downregulation of adherens junctions and in cancerous cell types. We anticipate that our system provides a simple and easy approach to investigate collective cell behavior in a well-controlled and systematic manner.

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.

Climb-enabled discrete dislocation plasticity

A small strain two-dimensional discrete dislocation plasticity framework coupled to vacancy diffusion is developed wherein the motion of edge dislocations is by a combination of glide and climb. The dislocations are modelled as line defects in a linear elastic medium and the mechanical boundary value problem is solved by the superposition of the infinite medium elastic fields of the dislocations and a complimentary non-singular solution that enforces the boundary conditions. Similarly, the climbing dislocations are modelled as line sources/sinks of vacancies and the vacancy diffusion boundary value problem is also solved by a superposition of the fields of the line sources/sinks in an infinite medium and a complementary non-singular solution that enforces the boundary conditions. The vacancy concentration field along with the stress field provides the climb rate of the dislocations. Other short-range interactions of the dislocations are incorporated via a set of constitutive rules. We first employ this formulation to investigate the climb of a single edge dislocation in an infinite medium and illustrate the existence of diffusion-limited and sink-limited climb regimes. Next, results are presented for the pure bending and uniaxial tension of single crystals oriented for single slip. These calculations show that plasticity size effects are reduced when dislocation climb is permitted. Finally, we contrast predictions of this coupled framework with an ad hoc model in which dislocation climb is modelled by a drag-type relation based on a quasi steady-state solution. © 2013 Elsevier Ltd. All rights reserved.

Glass reinforced glass

Laminated glass units are traditionally used to provide a degree of post-fracture strength, but the residual strength is often limited to relatively low levels suffi cient for holding the glass fragments together for a predetermined amount of time. It is possible to achieve a higher level of residual strength, but this requires specific boundary conditions and/or opaque reinforcing materials. This paper describes the experimental investigations on laminated glass units that can provide a signifi cant degree of post-fracture resistance, without the need of boundary restraints or opaque reinforcing materials. The glass units are composed entirely of combinations of conventional transparent interlayers and commercially available glass (annealed, heat treated and chemically strengthened). The paper also describes an empirical energy based interpretation of the mechanical response of the laminated units.

Bayesian parameter identification for vibro-acoustic models with nonparametric uncertainty

Vibration and acoustic analysis at higher frequencies faces two challenges: computing the response without using an excessive number of degrees of freedom, and quantifying its uncertainty due to small spatial variations in geometry, material properties and boundary conditions. Efficient models make use of the observation that when the response of a decoupled vibro-acoustic subsystem is sufficiently sensitive to uncertainty in such spatial variations, the local statistics of its natural frequencies and mode shapes saturate to universal probability distributions. This holds irrespective of the causes that underly these spatial variations and thus leads to a nonparametric description of uncertainty. This work deals with the identification of uncertain parameters in such models by using experimental data. One of the difficulties is that both experimental errors and modeling errors, due to the nonparametric uncertainty that is inherent to the model type, are present. This is tackled by employing a Bayesian inference strategy. The prior probability distribution of the uncertain parameters is constructed using the maximum entropy principle. The likelihood function that is subsequently computed takes the experimental information, the experimental errors and the modeling errors into account. The posterior probability distribution, which is computed with the Markov Chain Monte Carlo method, provides a full uncertainty quantification of the identified parameters, and indicates how well their uncertainty is reduced, with respect to the prior information, by the experimental data. © 2013 Taylor & Francis Group, London.

Static analysis of Euler-Bernoulli beams with multiple unilateral cracks under combined axial and transverse loads

The paper deals with the static analysis of pre-damaged Euler-Bernoulli beams with any number of unilateral cracks and subjected to tensile or compression forces combined with arbitrary transverse loads. The mathematical representation of cracks with a bilateral behaviour (i.e. always open) via Dirac delta functions is extended by introducing a convenient switching variable, which allows each crack to be open or closed depending on the sign of the axial strain at the crack centre. The proposed model leads to analytical solutions, which depend on four integration constants (to be computed by enforcing the boundary conditions) along with the Boolean switching variables associated with the cracks (whose role is to turn on and off the additional flexibility due to the presence of the cracks). An efficient computational procedure is also presented and numerically validated. For this purpose, the proposed approach is applied to two pre-damaged beams, with different damage and loading conditions, and the results so obtained are compared against those given by a standard finite element code (in which the correct opening of the cracks is pre-assigned), always showing a perfect agreement. © 2013 Elsevier Ltd. All rights reserved.

Non linear constitutive models for lattice materials

We use a computational homogenisation approach to derive a non linear constitutive model for lattice materials. A representative volume element (RVE) of the lattice is modelled by means of discrete structural elements, and macroscopic stress-strain relationships are numerically evaluated after applying appropriate periodic boundary conditions to the RVE. The influence of the choice of the RVE on the predictions of the model is discussed. The model has been used for the analysis of the hexagonal and the triangulated lattices subjected to large strains. The fidelity of the model has been demonstrated by analysing a plate with a central hole under prescribed in plane compressive and tensile loads, and then comparing the results from the discrete and the homogenised models. © 2013 Elsevier Ltd.