Modelling of geomagnetic induction in pipelines

International audienceGeomagnetic field variations induce telluric currents in pipelines, which modify the electrochemical conditions at the pipe/soil interface, possibly contributing to corrosion of the pipeline steel. Modelling of geomagnetic induction in pipelines can be accomplished by combining several techniques. Starting with geomagnetic field data, the geoelectric fields in the absence of the pipeline were calculated using the surface impedance derived from a layered-Earth conductivity model. The influence of the pipeline on the electric fields was then examined using an infinitely long cylinder (ILC) model. Pipe-to-soil potentials produced by the electric field induced in the pipeline were calculated using a distributed source transmission line (DSTL) model. The geomagnetic induction process is frequency dependent; therefore, the calculations are best performed in the frequency domain, using a Fourier transform to go from the original time domain magnetic data, and an inverse Fourier transform at the end of the process, to obtain the pipe-to-soil potential variation in the time domain. Examples of the model calculations are presented and compared to observations made on a long pipeline in the auroral zone

Solitary flexural–gravity waves in three dimensions

The focus of this work is on three-dimensional nonlinear flexural–gravity waves, propagating at the interface between a fluid and an ice sheet. The ice sheet is modelled using the special Cosserat theory of hyperelastic shells satisfying Kirchhoff's hypothesis, presented in (Plotnikov & Toland. 2011 Phil. Trans. R. Soc. A 369, 2942–2956 (doi:10.1098/rsta.2011.0104)). The fluid is assumed inviscid and incompressible, and the flow irrotational. A numerical method based on boundary integral equation techniques is used to compute solitary waves and forced waves to Euler's equations. This article is part of the theme issue ‘Modelling of sea-ice phenomena’

Stability of periodic traveling flexural‐gravity waves in two dimensions

In this work, we solve the Euler’s equations for periodic waves travelling under a sheet of ice using a reformulation introduced in [1]. These waves are referred to as flexural-gravity waves. We compare and contrast two models for the effect of the ice: a linear model and a nonlinear model. The benefit of this reformulation is that it facilitates the asymptotic analysis. We use it to derive the nonlinear Schrödinger equation that describes the modulational instability of periodic travelling waves. We compare this asymptotic result with the numerical computation of stability using the Fourier-Floquet-Hill method to show they agree qualitatively. We show that different models have different stability regimes for large values of the flexural rigidity parameter. Numerical computations are also used to analyse high frequency instabilities in addition to the modulational instability. In the regions examined, these are shown to be the same regardless of the model representing ice

Earth conductivity structures and their effects on geomagnetic induction in pipelines

Anomalous, large pipe-to-soil potentials (PSP) have been observed along a natural gas pipeline in eastern Ontario, Canada, where there is a major geological contact between the highly resistive rocks of the Precambrian Shield to the west and the more conductive Paleozoic sediments to the east. This study tested the hypothesis that large variations of PSP are related to lateral changes of Earth conductivity under the pipeline. Concurrent and co-located PSP and magnetotelluric (MT) geophysical data were acquired in the study area. Results from the MT survey were used to model PSP variations based on distributed-source transmission line theory, using a spatially-variant surface geoelectric field. Different models were built to investigate the impact of different subsurface features. Good agreement between modelled and observed PSP was reached when impedance peaks related to major changes of subsurface geological conditions were included. The large PSP could therefore be attributed to the presence of resistive intrusive bodies in the upper crust and/or boundaries between tectonic terranes. This study demonstrated that combined PSP-MT investigations are a useful tool in the identification of potential hazards caused by geomagnetically induced currents in pipelines

Modelling Li+ Ion Battery Electrode Properties

We formulated two detailed models for an electrolytic cell with particulate electrodes based on a lithium atom concentration dependent Butler-Volmer condition at the interface between electrode particles and the electrolyte. The first was based on a dilute-ion assumption for the electrolyte, while the second assumed that Li ions are present in excess. For the first, we used the method of multiple scales to homogenize this model over the microstructure, formed by the small lithium particles in the electrodes. For the second, we gave rigorous bounds for the effective electrochemical conductivity for a linearized case. We expect similar results and bounds for the "full nonlinear problem" because variational results are generally not adversely affected by a sinh term. Finally we used the asymptotic methods, based on parameters estimated from the literature, to attain a greatly simplified one-dimensional version of the original homogenized model. This simplified model accounts for the fact that diffusion of lithium atoms within individual electrode particles is relatively much faster than that of lithium ions across the whole cell so that lithium ion diffusion is what limits the performance of the battery. However, since most of the potential drop occurs across the Debye layers surrounding each electrode particle, lithium ion diffusion only significantly affects cell performance if there is more or less complete depletion of lithium ions in some region of the electrolyte which causes a break in the current flowing across the cell. This causes catastrophic failure. Providing such failure does not occur the potential drop across the cell is determined by the concentration of lithium atoms in the electrode particles. Within each electrode lithium atom concentration is, to leading order, a function of time only and not of position within the electrode. The depletion of electrode lithium atom concentration is directly proportional to the current being drawn off the cell. This leads one to expect that the potential of the cell gradually drops as current is drawn of it. We would like to emphasize that all the homogenization methods employed in this work give a systematic approach for investigating the effect that changes in the microstructure have on the behaviour of the battery. However, due to lack of time, we have not used this method to investigate particular particle geometries