Asymmetric propagation of airblast from bench blasting

This paper investigates the propagation of airblast from quarry blasting. Peak overpressure is calculated as a function of blasting parameters (explosive mass per delay and velocity at which the detonation sequence proceeds along the bench) and polar coordinates of the point of interest (distance to the blast and azimuth with respect to the free face of the blast). The model is in the form of the product of a classical scaled distance attenuation law times a directional correction factor. The latter considers the influence of the bench face, and attenuates overpressure at the top level and amplifies it at the bottom. Such factor also accounts for the effect of the delay by amplifying the pressure in the direction of the initiation sequence if the velocity of initiation exceeds half the speed of sound and up to an initiation velocity in the range of the speed of sound. The model has been fitted to an empirical data set composed by 134 airblast records monitored in 47 blasts at two quarries. The measurements were made at distances to the blast less than 450 m. The model is statistically significant and has a determination coefficient of 0.869

Coupled model of initiation and propagation of corrosion in reinforced concrete

Corrosion of reinforcing steel in concrete due to chloride ingress is one of the main causes of the deterioration of reinforced concrete structures. Structures most affected by such a corrosion are marine zone buildings and structures exposed to de-icing salts like highways and bridges. Such process is accompanied by an increase in volume of the corrosión products on the rebarsconcrete interface. Depending on the level of oxidation, iron can expand as much as six times its original volume. This increase in volume exerts tensile stresses in the surrounding concrete which result in cracking and spalling of the concrete cover if the concrete tensile strength is exceeded. The mechanism by which steel embedded in concrete corrodes in presence of chloride is the local breakdown of the passive layer formed in the highly alkaline condition of the concrete. It is assumed that corrosion initiates when a critical chloride content reaches the rebar surface. The mathematical formulation idealized the corrosion sequence as a two-stage process: an initiation stage, during which chloride ions penetrate to the reinforcing steel surface and depassivate it, and a propagation stage, in which active corrosion takes place until cracking of the concrete cover has occurred. The aim of this research is to develop computer tools to evaluate the duration of the service life of reinforced concrete structures, considering both the initiation and propagation periods. Such tools must offer a friendly interface to facilitate its use by the researchers even though their background is not in numerical simulation. For the evaluation of the initiation period different tools have been developed: Program TavProbabilidade: provides means to carry out a probability analysis of a chloride ingress model. Such a tool is necessary due to the lack of data and general uncertainties associated with the phenomenon of the chloride diffusion. It differs from the deterministic approach because it computes not just a chloride profile at a certain age, but a range of chloride profiles for each probability or occurrence. Program TavProbabilidade_Fiabilidade: carries out reliability analyses of the initiation period. It takes into account the critical value of the chloride concentration on the steel that causes breakdown of the passive layer and the beginning of the propagation stage. It differs from the deterministic analysis in that it does not predict if the corrosion is going to begin or not, but to quantifies the probability of corrosion initiation. Program TavDif_1D: was created to do a one dimension deterministic analysis of the chloride diffusion process by the finite element method (FEM) which numerically solves Fick’second Law. Despite of the different FEM solver already developed in one dimension, the decision to create a new code (TavDif_1D) was taken because of the need to have a solver with friendly interface for pre- and post-process according to the need of IETCC. An innovative tool was also developed with a systematic method devised to compare the ability of the different 1D models to predict the actual evolution of chloride ingress based on experimental measurements, and also to quantify the degree of agreement of the models with each others. For the evaluation of the entire service life of the structure: a computer program has been developed using finite elements method to do the coupling of both service life periods: initiation and propagation. The program for 2D (TavDif_2D) allows the complementary use of two external programs in a unique friendly interface: • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. This program (TavDif_2D) is responsible to decide in each time step when and where to start applying the boundary conditions of fracture mechanics module in function of the amount of chloride concentration and corrosion parameters (Icorr, etc). This program is also responsible to verify the presence and the degree of fracture in each element to send the Information of diffusion coefficient variation with the crack width. • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. The advantages of the FEM with the interface provided by the tool are: • the flexibility to input the data such as material property and boundary conditions as time dependent function. • the flexibility to predict the chloride concentration profile for different geometries. • the possibility to couple chloride diffusion (initiation stage) with chemical and mechanical behavior (propagation stage). The OOFEM code had to be modified to accept temperature, humidity and the time dependent values for the material properties, which is necessary to adequately describe the environmental variations. A 3-D simulation has been performed to simulate the behavior of the beam on both, action of the external load and the internal load caused by the corrosion products, using elements of imbedded fracture in order to plot the curve of the deflection of the central region of the beam versus the external load to compare with the experimental data.

Comparison of multi-year Ka-band propagation characteristics at four sites in North America and Europe

Tropospheric phenomena such as clouds and mainly rain cause higher attenuation at Ka-band than at lower frequencies. In this collaborative paper, the main results of four long-term Ka-band propagation campaigns are presented. The experiments are carried out in Ottawa, Canada (satellite Anik F2); Aveiro, Portugal; Madrid, Spain; and Toulouse, France (satellite HotBird 6 in the last three cases) and have been running since 2004 in Aveiro, 2006 in Ottawa and Madrid, and 2008 in Toulouse. After a brief introduction of the experiments, rain rate and excess attenuation results are discussed, first for a common two-year measurement period and later for the whole database available. Seasonal attenuation statistics for Madrid, Ottawa and Aveiro are compared. Finally, fade duration and fade slope statistics derived at three locations are presented and discussed.

Modelling of short runout propagation landslides and debris flows

This paper proposes an extension of methods used to predict the propagation of landslides having a long runout to smaller landslides with much shorter propagation distances. The method is based on: (1) a depth-integrated mathematical model including the coupling between the soil skeleton and the pore fluids, (2) suitable rheological models describing the relation between the stress and the rate of deformation tensors for fluidised soils and (3) a meshless numerical method, Smooth Particle Hydrodynamics, which separates the computational mesh (or set of computational nodes) from the mesh describing the terrain topography, which is of structured type – thus accelerating search operations. The proposed model is validated using two examples for which there are analytical solutions, and then it is applied to two short runout landslides which happened in Hong Kong in 1995, for which there is available information.

Neurite, a finite difference large scale parallel program for the simulation of the electrical signal propagation in neurites under mechanical loading

With the growing body of research on traumatic brain injury and spinal cord injury, computational neuroscience has recently focused its modeling efforts on neuronal functional deficits following mechanical loading. However, in most of these efforts, cell damage is generally only characterized by purely mechanistic criteria, function of quantities such as stress, strain or their corresponding rates. The modeling of functional deficits in neurites as a consequence of macroscopic mechanical insults has been rarely explored. In particular, a quantitative mechanically based model of electrophysiological impairment in neuronal cells has only very recently been proposed (Jerusalem et al., 2013). In this paper, we present the implementation details of Neurite: the finite difference parallel program used in this reference. Following the application of a macroscopic strain at a given strain rate produced by a mechanical insult, Neurite is able to simulate the resulting neuronal electrical signal propagation, and thus the corresponding functional deficits. The simulation of the coupled mechanical and electrophysiological behaviors requires computational expensive calculations that increase in complexity as the network of the simulated cells grows. The solvers implemented in Neurite-explicit and implicit-were therefore parallelized using graphics processing units in order to reduce the burden of the simulation costs of large scale scenarios. Cable Theory and Hodgkin-Huxley models were implemented to account for the electrophysiological passive and active regions of a neurite, respectively, whereas a coupled mechanical model accounting for the neurite mechanical behavior within its surrounding medium was adopted as a link between lectrophysiology and mechanics (Jerusalem et al., 2013). This paper provides the details of the parallel implementation of Neurite, along with three different application examples: a long myelinated axon, a segmented dendritic tree, and a damaged axon. The capabilities of the program to deal with large scale scenarios, segmented neuronal structures, and functional deficits under mechanical loading are specifically highlighted.

Elemental propagation of calcium signals in response-specific patterns determined by environmental stimulus strength

Plant cells can respond qualitatively and quantitatively to a wide range of environmental signals. Ca2+ is used as an intracellular signal for volume regulation in response to external osmotic changes. We show here that the spatiotemporal patterns of hypo-osmotically induced Ca2+ signals vary dramatically with stimulus strength in embryonic cells of the marine alga Fucus. Biphasic or multiphasic Ca2+ signals reflect Ca2+ elevations in distinct cellular domains. These propagate via elemental Ca2+ release in nuclear or peripheral regions that are rich in endoplasmic reticulum. Cell volume regulation specifically requires Ca2+ elevation in apical peripheral regions, whereas an altered cell division rate occurs only in response to stimuli that cause Ca2+ elevation in nuclear regions.

Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage

Ribozymes of hepatitis delta virus have been proposed to use an active-site cytosine as an acid-base catalyst in the self-cleavage reaction. In this study, we have examined the role of cytosine in more detail with the antigenomic ribozyme. Evidence that proton transfer in the rate-determining step involved cytosine 76 (C76) was obtained from examining cleavage activity of the wild-type and imidazole buffer-rescued C76-deleted (C76Δ) ribozymes in D2O and H2O. In both reactions, a similar kinetic isotope effect and shift in the apparent pKa indicate that the buffer is functionally substituting for the side chain in proton transfer. Proton inventory of the wild-type reaction supported a mechanism of a single proton transfer at the transition state. This proton transfer step was further characterized by exogenous base rescue of a C76Δ mutant with cytosine and imidazole analogues. For the imidazole analogues that rescued activity, the apparent pKa of the rescue reaction, measured under kcat/KM conditions, correlated with the pKa of the base. From these data a Brønsted coefficient (β) of 0.51 was determined for the base-rescued reaction of C76Δ. This value is consistent with that expected for proton transfer in the transition state. Together, these data provide strong support for a mechanism where an RNA side chain participates directly in general acid or general base catalysis of the wild-type ribozyme to facilitate RNA cleavage.

Procedure for calculating the chloride diffusion coefficient and surface concentration from a profile having a maximum beyond the concrete surface

Chlorides induce local corrosion in the steel reinforcements when reaching the bar surface. The measurement of the rate of ingress of these ions, is made by mathematically fitting the so called “error function equation” into the chloride concentration profile, obtaining so the diffusion coefficient and the chloride concentration at the concrete surface. However, the chloride profiles do not always follow Fick’s law by having the maximum concentration at the concrete surface, but often the profile shows a maximum concentration more in the interior, which indicates a different composition and performance of the most external concrete layer with respect to the internal zones. The paper presents a procedure prepared during the time of the RILEM TC 178-TMC: “Testing and modeling chloride penetration in concrete”, which suggests neglecting the external layer where the chloride concentration increases and using the maximum as an “apparent” surface concentration, called C max and to fit the error function equation into the decreasing concentration profile towards the interior. The prediction of evolution should be made also from the maximum.

(Table 1) Results from XRD analysis as well as sediment type and selected results from ring shear experiments

The location of the seaward tip of a subduction thrust controls material transfer at convergent plate margins, and hence global mass balances. At approximately half of those margins, the material of the subducting plate is completely underthrust so that no accretion or even subduction erosion takes place. Along the remaining margins, material is scraped off the subducting plate and added to the upper plate by frontal accretion. We here examine the physical properties of subducting sediments off Costa Rica and Nankai, type examples for an erosional and an accretionary margin, to investigate which parameters control the level where the frontal thrust cuts into the incoming sediment pile. A series of rotary-shear experiments to measure the frictional strength of the various lithologies entering the two subduction zones were carried out. Results include the following findings: (1) At Costa Rica, clay-rich strata at the top of the incoming succession have the lowest strength (µres = 0.19) while underlying calcareous ooze, chalk and diatomite are strong (up to µres = 0.43; µpeak = 0.56). Hence the entire sediment package is underthrust. (2) Off Japan, clay-rich deposits within the lower Shikoku Basin inventory are weakest (µres = 0.13-0.19) and favour the frontal proto-thrust to migrate into one particular horizon between sandy, competent turbidites below and ash-bearing mud above. (3) Taking in situ data and earlier geotechnical testing into account, it is suggested that mineralogical composition rather than pore-pressure defines the position of the frontal thrust, which locates in the weakest, clay mineral-rich (up to 85 wt.%) materials. (4) Smectite, the dominant clay mineral phase at either margin, shows rate strengthening and stable sliding in the frontal 50 km of the subduction thrust (0.0001-0.1 mm/s, 0.5-25 MPa effective normal stress). (5) Progressive illitization of smectite cannot explain seismogenesis, because illite-rich samples also show velocity strengthening at the conditions tested.

Interpretation of the temperature dependence of rate constants in biosensor studies

A comparison is made between Arrhenius and transition-state analyses of the temperature dependence of rate constants reported in four published biosensor studies. Although the Eyring transition-state theory seemingly affords a more definitive solution to the problem of characterizing the activation energetics, the analysis is equivocal because of inherent assumptions about reaction mechanism and the magnitude of the transmission coefficient. In view of those uncertainties it is suggested that a preferable course of action entails reversion to the empirical Arrhenius analysis with regard to the energy of activation and a preexponential factor. The former is essentially equivalent to the enthalpy of activation, whereas the magnitude of the latter indicates directly the extent of disparity between the frequency of product formation and the universal frequency factor (temperature multiplied by the ratio of the Boltzmann and Planck constants) and hence the likelihood of a more complicated kinetic mechanism than that encompassed by the Eyring transition-state theory. (C) 2004 Elsevier Inc. All rights reserved.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.

1220km propagation of 40Gbit/s single channel RZ data over dispersion managed standard (non-dispersion shifted) fibre

Error free propagation of a single polarisation optical time division multiplexed 40 Gbit/s dispersion managed pulsed data stream over dispersion (non-shifted) fibre. This distance is twice the previous record at this data rate.

Inscription and characterization of waveguides written into borosilicate glass by a high-repetition-rate femtosecond laser at 800 nm

A series of waveguides was inscribed in a borosilicate glass (BK7) by an 11 MHz repetition rate femtosecond laser operating with pulse energies from 16 to 30 nJ and focused at various depths within the bulk material. The index modification was measured using a quantitative phase microscopy technique that revealed central index changes ranging from 5×10-3 to 10-2, leading to waveguides that exhibited propagation losses of 0.2 dB/cm at a wavelength of 633 nm and 0.6 dB/cm at a wavelength of 1550 nm with efficient mode matching, less than 0.2 dB, to standard optical fibers. Analysis of the experimental data shows that, for a given inscription energy, the index modification has a strong dependence on inscription scanning velocity. At higher energies, the index modification increases with increasing inscription scanning velocity with other fabrication parameters constant.

Impact of longitudinal power budget in coherent transmission systems employing digital back-propagation

We report the impact of longitudinal signal power profile on the transmission performance of coherently-detected 112 Gb/s m-ary polarization multiplexed quadrature amplitude modulation system after compensation of deterministic nonlinear fibre impairments. Performance improvements up to 0.6 dB (Q(eff)) are reported for a non-uniform transmission link power profile. Further investigation reveals that the evolution of the transmission performance with power profile management is fully consistent with the parametric amplification of the amplified spontaneous emission by the signal through four-wave mixing. In particular, for a non-dispersion managed system, a single-step increment of 4 dB in the amplifier gain, with respect to a uniform gain profile, at similar to 2/3(rd) of the total reach considerably improves the transmission performance for all the formats studied. In contrary a negative-step profile, emulating a failure (gain decrease or loss increase), significantly degrades the bit-error rate.

Nonlinearity compensation in multi-rate 28 Gbaud WDM systems employing optical and digital techniques under diverse link configurations

Digital back-propagation (DBP) has recently been proposed for the comprehensive compensation of channel nonlinearities in optical communication systems. While DBP is attractive for its flexibility and performance, it poses significant challenges in terms of computational complexity. Alternatively, phase conjugation or spectral inversion has previously been employed to mitigate nonlinear fibre impairments. Though spectral inversion is relatively straightforward to implement in optical or electrical domain, it requires precise positioning and symmetrised link power profile in order to avail the full benefit. In this paper, we directly compare ideal and low-precision single-channel DBP with single-channel spectral-inversion both with and without symmetry correction via dispersive chirping. We demonstrate that for all the dispersion maps studied, spectral inversion approaches the performance of ideal DBP with 40 steps per span and exceeds the performance of electronic dispersion compensation by ~3.5 dB in Q-factor, enabling up to 96% reduction in complexity in terms of required DBP stages, relative to low precision one step per span based DBP. For maps where quasi-phase matching is a significant issue, spectral inversion significantly outperforms ideal DBP by ~3 dB.