On the Convexity of Nonlinear Elastic Energies in the Right Cauchy-Green Tensor

We present a sufficient condition under which a weak solution of the Euler-Lagrange equations in nonlinear elasticity is already a global minimizer of the corresponding elastic energy functional. This criterion is applicable to energies \(W(F)=\widehat{W}(F^{T}F)=\widehat{W}(C)\) which are convex with respect to the right Cauchy-Green tensor \(C=F^{T}F\), where \(F\) denotes the gradient of deformation. Examples of such energies exhibiting a blow up for \(\det F\to0\) are given.     

Long-time Dynamics of Resonant Weakly Nonlinear CGL Equations

Consider a weakly nonlinear CGL equation on the torus \(\mathbb {T}^d\):

$$\begin{aligned} u_t+i\Delta u=\epsilon [\mu (-1)^{m-1}\Delta ^{m} u+b|u|^{2p}u+ ic|u|^{2q}u]. \end{aligned}$$

(*)

Here \(u=u(t,x)\), \(x\in \mathbb {T}^d\), \(0<\epsilon <<1\), \(\mu \geqslant 0\), \(b,c\in \mathbb {R}\) and \(m,p,q\in \mathbb {N}\). Define \(I(u)=(I_{\mathbf {k}},\mathbf {k}\in \mathbb {Z}^d)\), where \(I_{\mathbf {k}}=v_{\mathbf {k}}\bar{v}_{\mathbf {k}}/2\) and \(v_{\mathbf {k}}\), \(\mathbf {k}\in \mathbb {Z}^d\), are the Fourier coefficients of the function \(u\) we give. Assume that the equation \((*)\) is well posed on time intervals of order \(\epsilon ^{-1}\) and its solutions have there a-priori bounds, independent of the small parameter. Let \(u(t,x)\) solve the equation \((*)\). If \(\epsilon \) is small enough, then for \(t\lesssim {\epsilon ^{-1}}\), the quantity \(I(u(t,x))\) can be well described by solutions of an effective equation:

$$\begin{aligned} u_t=\epsilon [\mu (-1)^{m-1}\Delta ^m u+ F(u)], \end{aligned}$$

where the term \(F(u)\) can be constructed through a kind of resonant averaging of the nonlinearity \(b|u|^{2p}+ ic|u|^{2q}u\).

     

On the Justification of Viscoelastic Elliptic Membrane Shell Equations

We consider a family of linearly viscoelastic shells with thickness \(2\varepsilon\), clamped along their entire lateral face, all having the same middle surface \(S=\boldsymbol{\theta}(\bar{\omega})\subset \mathbb{R}^{3}\), where \(\omega\subset\mathbb{R}^{2}\) is a bounded and connected open set with a Lipschitz-continuous boundary \(\gamma\). We make an essential geometrical assumption on the middle surface \(S\), which is satisfied if \(\gamma\) and \(\boldsymbol{\theta}\) are smooth enough and \(S\) is uniformly elliptic. We show that, if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), the solution of the scaled variational problem in curvilinear coordinates, \(\boldsymbol{u}( \varepsilon)\), defined over the fixed domain \(\varOmega=\omega\times (-1,1)\) for each \(t\in[0,T]\), converges to a limit \(\boldsymbol{u}\) with \(u_{\alpha}(\varepsilon)\rightarrow u_{\alpha}\) in \(W^{1,2}(0,T,H ^{1}(\varOmega))\) and \(u_{3}(\varepsilon)\rightarrow u_{3}\) in \(W^{1,2}(0,T,L^{2}(\varOmega))\) as \(\varepsilon\to0\). Moreover, we prove that this limit is independent of the transverse variable. Furthermore, the average \(\bar{\boldsymbol{u}}= \frac{1}{2}\int_{-1}^{1} \boldsymbol{u}dx_{3}\), which belongs to the space \(W^{1,2}(0,T, V_{M}( \omega))\), where

$$V_{M}(\omega)=H^{1}_{0}(\omega)\times H^{1}_{0}(\omega)\times L ^{2}(\omega), $$

satisfies what we have identified as (scaled) two-dimensional equations of a viscoelastic membrane elliptic shell, which includes a long-term memory that takes into account previous deformations. We finally provide convergence results which justify those equations.

     

Boundedness for Some Doubly Nonlinear Parabolic Equations in Measure Spaces

In the context of measure spaces equipped with a doubling non-trivial Borel measure supporting a Poincaré inequality, we derive local and global sup bounds of the nonnegative weak subsolutions of

$$\begin{aligned} (u^{q})_t-\nabla \cdot {(|\nabla u|^{p-2}\nabla u)}=0, \quad \mathrm {in} \ U_\tau = U \times (\tau _1, \tau _2] , \quad p>1,\quad q>1 \end{aligned}$$

and of its associated Dirichlet problem, respectively. For particular ranges of the exponents p and q, we show that any locally nonnegative weak subsolution, taken in \(Q (\subset \bar{Q}\subset U_\tau )\), is controlled from above by the \(L^\alpha (\bar{Q}) \)-norm, for \(\alpha = \max \{p, q+1\}\). As for the global setting, under suitable assumptions on the boundary datum g and on the initial datum, we obtain sup bounds for u, in \(U \times \{ t\}\), which depend on the \(\sup g\) and on the \(L^{q+1}(U \times (\tau _1, \tau _1+t])\)-norm of \((u-\sup g)_+\), for all \(t \in (0, \tau _2-\tau _1]\). On the critical ranges of p and q, a priori local and global \(L^\infty \) estimates require extra qualitative information on u.

     

On the Damped Harmonic Oscillator with Time Dependent Damping Coefficient

Conditions guaranteeing asymptotic stability for the differential equation

$$\begin{aligned} x''+h(t)x'+\omega ^2x=0 \qquad (x\in \mathbb {R}) \end{aligned}$$

are studied, where the damping coefficient \(h:[0,\infty )\rightarrow [0,\infty )\) is a locally integrable function, and the frequency \(\omega >0\) is constant. Our conditions need neither the requirement \(h(t)\le \overline{h}<\infty \) (\(t\in [0,\infty )\); \(\overline{h}\) is constant) (“small damping”), nor \(0< \underline{h}\le h(t)\) (\(t\in [0,\infty )\); \(\underline{h}\) is constant) (“large damping”); in other words, they can be applied to the general case \(0\le h(t)<\infty \) (\(t\in [0,\infty \))). We establish a condition which combines weak integral positivity with Smith’s growth condition

$$\begin{aligned} \int ^\infty _0 \exp [-H(t)]\int _0^t \exp [H(s)]\,\mathrm{{d}}s\,\mathrm{{d}}t=\infty \qquad \left( H(t):=\int _0^t h(\tau )\,\mathrm{{d}}\tau \right) , \end{aligned}$$

so it is able to control both the small and the large values of the damping coefficient simultaneously.

     

One-Dimensional Super-Fast Diffusion: Persistence Versus Extinction Revisited—Extinction at Spatial Infinity

We consider positive classical solutions of

$$\begin{aligned} v_t=(v^{m-1}v_x)_x, \qquad x\in {\mathbb {R}}, \ t>0, \qquad (\star ) \end{aligned}$$

in the super-fast diffusion range \(m<-1\). Our main interest is in smooth positive initial data \(v_0=v(\cdot ,0)\) which decay as \(x\rightarrow +\infty \), but which are possibly unbounded as \(x\rightarrow -\infty \), having in mind monotonically decreasing data as prototypes. It is firstly proved that if \(v_0\) decays sufficiently fast only in one direction by satisfying

$$\begin{aligned} v_0(x) \le cx^{-\beta } \qquad \text{ for } \text{ all } ~x>0 \quad \hbox { with some }\quad \beta >\frac{2}{1-m} \end{aligned}$$

and some \(c>0\), then the so-called proper solution of (\(\star \)) vanishes identically in \({\mathbb {R}}\times (0,\infty )\), and accordingly no positive classical solution exists in any time interval in this case. Complemented by some sufficient criteria for solutions to remain positive either locally or globally in time, this condition for instantaneous extinction is shown to be optimal at least with respect to algebraic decay of the initial data. This partially extends some known nonexistence results for (\(\star \)) (Daskalopoulos and Del Pino in Arch Rat Mech Anal 137(4):363–380, 1997) in that it does not require any knowledge on the behavior of \(v_0(x)\) for \(x<0\). Next focusing on the phenomenon of extinction in finite time, we show that in this respect a mass influx from \(x=-\infty \) can interact with mass loss at \(x=+\infty \) in a nontrivial manner. Namely, we shall detect examples of monotone initial data, with critical decay as \(x\rightarrow +\infty \) and exponential growth as \(x\rightarrow -\infty \), that lead to solutions of (\(\star \)) which become extinct at a finite positive time, but which have empty extinction sets. This is in sharp contrast to known extinction mechanisms which are such that the corresponding extinction sets coincide with all of \({\mathbb {R}}\).

     

Models of Elastic Shells in Contact with a Rigid Foundation: An Asymptotic Approach

We consider a family of linearly elastic shells with thickness \(2\varepsilon\) (where \(\varepsilon\) is a small parameter). The shells are clamped along a portion of their lateral face, all having the same middle surface \(S\), and may enter in contact with a rigid foundation along the bottom face.We are interested in studying the limit behavior of both the three-dimensional problems, given in curvilinear coordinates, and their solutions (displacements \(\boldsymbol{u}^{\varepsilon}\) of covariant components \(u_{i}^{\varepsilon}\)) when \(\varepsilon\) tends to zero. To do that, we use asymptotic analysis methods. On one hand, we find that if the applied body force density is \(O(1)\) with respect to \(\varepsilon\) and surface tractions density is \(O(\varepsilon)\), a suitable approximation of the variational formulation of the contact problem is a two-dimensional variational inequality which can be identified as the variational formulation of the obstacle problem for an elastic membrane. On the other hand, if the applied body force density is \(O(\varepsilon^{2})\) and surface tractions density is \(O(\varepsilon^{3})\), the corresponding approximation is a different two-dimensional inequality which can be identified as the variational formulation of the obstacle problem for an elastic flexural shell. We finally discuss the existence and uniqueness of solution for the limit two-dimensional variational problems found.     

Multiple Solutions to Neumann Problems with Indefinite Weight and Bounded Nonlinearities

We study the Neumann boundary value problem for the second order ODE

$$\begin{aligned} u^{\prime \prime } + (a^+(t)-\mu a^-(t))g(u) = 0, \qquad t \in [0,T], \end{aligned}$$

(1)

where \(g \in {\mathcal {C}}^1({\mathbb {R}})\) is a bounded function of constant sign, \(a^+,a^-: [0,T] \rightarrow {\mathbb {R}}^+\) are the positive/negative part of a sign-changing weight \(a(t)\) and \(\mu > 0\) is a real parameter. Depending on the sign of \(g^{\prime }(u)\) at infinity, we find existence/multiplicity of solutions for \(\mu \) in a “small” interval near the value

$$\begin{aligned} \mu _c = \frac{\int _0^T a^+(t) \, dt}{\int _0^T a^-(t) \, dt}\,. \end{aligned}$$

The proof exploits a change of variables, transforming the sign-indefinite Eq. (1) into a forced perturbation of an autonomous planar system, and a shooting argument. Nonexistence results for \(\mu \rightarrow 0^+\) and \(\mu \rightarrow +\infty \) are given, as well.

     

Synchronization transitions induced by partial time delay in a excitatory–inhibitory coupled neuronal network

Information transmission delays are an inherent factor of neuronal systems as a consequence of the finite propagation speeds and time lapses occurring by both dendritic and synaptic processes. In real neuronal systems, some delay between two neurons is too small and can be ignored, which results in partial time delay. In this paper, we focus on investigating influences of partial time delay on synchronization transitions in a excitatory–inhibitory (E–I) coupled neuronal networks. Here, we suppose time delay between two neurons equals to \(\tau \) with probability \(p_{\mathrm{delay}}\) and investigate effect of partial time delay on synchronization transitions of the neuronal networks by controlling \(\tau \) and \(p_{\mathrm{delay}}\) under three cases. In these three cases, excitatory synapses are always considered to delayed with probability \(p_{\mathrm{delay}}\), while inhibitory synapses are considered to be without delays (case I), delayed with probability \(p_{\mathrm{delay}}\) (case II), and always delayed (case III), respectively. It is revealed that, in the first two cases, partial time delay has little influences on synchronization of the neuronal network for small \(p_{\mathrm{delay}}\), while it could induce synchronization transitions at \(\tau \) around integer multiples of the period of individual neuron T when \(p_{\mathrm{delay}}\) is large enough, while in the case III, partial time delay could induce synchronization transitions at \(\tau \) being around odd integer multiples of T / 2 for small \(p_{\mathrm{delay}}\) and at \(\tau \) being around integer multiples of T for large \(p_{\mathrm{delay}}\). Most interesting observation is that partial time delay could induce frequent synchronization transitions at \(\tau \) being around integer multiples of T / 2 for intermediate \(p_{\mathrm{delay}}\). Moreover, effect of rewiring probability on synchronization transitions induced by partial time delay has been discussed. It is found that synchronization transitions induced by partial time delay are robust to rewiring probability for large \(p_{\mathrm{delay}}\) under the three cases.     

A novel approach for generating multi-direction multi-double-scroll attractors

Attention is focused in this work on quasiperiodic motion of nonlinear systems whose spectrum contains uniformly spaced sideband frequencies with a distance \(\omega _{d}\) apart, around a frequency \(\omega \) with \(\omega \gg \omega _{d}\) and its integer multiples, which are referred to as carrier frequencies. The ratio of the two frequencies \(\omega \) and \(\omega _{d}\) is an irrational number. A new method based on the traditional incremental harmonic balance (IHB) method with multiple timescales, referred to as Lau method, where two timescales, \(\tau _{1}=\omega t\) (a fast timescale) and \(\tau _{2}=\omega _{d}t\) (a slow timescale), are introduced, is presented to analyze quasiperiodic motion of nonlinear systems. An amplitude increment algorithm is adapted to deal with cases where the two frequencies \(\omega \) and \(\omega _{d}\) are    unknown a priori, in order to automatically trace frequency response of quasiperiodic motion of nonlinear systems and accurately calculate all frequency components and their corresponding amplitudes. Results of application of the present IHB method to quasiperiodic free vibration of a hinged–clamped beam with internal resonance between two transverse modes are shown and compared with previously published results with Lau method and those from numerical integration. While differences are noted between results predicted by the present IHB method and Lau method, excellent agreement is achieved between results from the present IHB method and numerical integration even in cases of strongly nonlinear vibration. The present IHB method is also used to analyze quasiperiodic free vibration of high-dimensional models of the hinged–clamped beam.     

Chaotic characteristics analysis of the sintering process system with unknown dynamic functions based on phase space reconstruction and chaotic invariables

We develop a local discontinuous Galerkin finite element method for the distributed-order time and Riesz space-fractional convection–diffusion and Schrödinger-type equations. The stability of the presented schemes is proved and optimal order of convergence \(\mathcal {O}(h^{N+1}+(\Delta t)^{1+\frac{\theta }{2}}+\theta ^{2})\) for the Riesz space-fractional diffusion and Schrödinger-type equations with distributed order in time, an order of convergence of \(\mathcal {O}(h^{N+\frac{1}{2}}+(\Delta t)^{1+\frac{\theta }{2}}\) \(+\theta ^{2})\) is provided for the Riesz space-fractional convection–diffusion equations with distributed order in time where h, \(\theta \) and \(\Delta t\) are space step size, the distributed-order variables and the step sizes in time, respectively. Finally, the performed numerical examples confirm the optimal convergence order and illustrate the effectiveness of the method.     

Spreading in a Shifting Environment Modeled by the Diffusive Logistic Equation with a Free Boundary

We investigate the influence of a shifting environment on the spreading of an invasive species through a model given by the diffusive logistic equation with a free boundary. When the environment is homogeneous and favourable, this model was first studied in Du and Lin (SIAM J Math Anal 42:377–405, 2010), where a spreading–vanishing dichotomy was established for the long-time dynamics of the species, and when spreading happens, it was shown that the species invades the new territory at some uniquely determined asymptotic speed \(c_0>0\). Here we consider the situation that part of such an environment becomes unfavourable, and the unfavourable range of the environment moves into the favourable part with speed \(c>0\). We prove that when \(c\ge c_0\), the species always dies out in the long-run, but when \(0<c<c_0\), the long-time behavior of the species is determined by a trichotomy described by (a) vanishing, (b) borderline spreading, or (c) spreading. If the initial population is written in the form \(u_0(x)=\sigma \phi (x)\) with \(\phi \) fixed and \(\sigma >0\) a parameter, then there exists \(\sigma _0>0\) such that vanishing happens when \(\sigma \in (0,\sigma _0)\), borderline spreading happens when \(\sigma =\sigma _0\), and spreading happens when \(\sigma >\sigma _0\).     

SIF-based fracture criterion for interface cracks

The complex stress intensity factor K governing the stress field of an interface crack tip may be split into two parts, i.e.,■ and s~(-iε), so that K = ■ s~(-iε), s is a characteristic length and ε is the oscillatory index. ■ has the same dimension as the classical stress intensity factor and characterizes the interface crack tip field. That means a criterion for interface cracks may be formulated directly with■, as Irwin(ASME J. Appl. Mech. 24:361–364, 1957) did in 1957 for the classical fracture mechanics. Then, for an interface crack,it is demonstrated that the quasi Mode I and Mode II tip fields can be defined and distinguished from the coupled mode tip fields. Built upon SIF-based fracture criteria for quasi Mode I and Mode II, the stress intensity factor(SIF)-based fracture criterion for mixed mode interface cracks is proposed and validated against existing experimental results.     

Threshold Behavior and Non-quasiconvergent Solutions with Localized Initial Data for Bistable Reaction–Diffusion Equations

We consider bounded solutions of the semilinear heat equation \(u_t=u_{xx}+f(u)\) on \(R\), where \(f\) is of the unbalanced bistable type. We examine the \(\omega \)-limit sets of bounded solutions with respect to the locally uniform convergence. Our goal is to show that even for solutions whose initial data vanish at \(x=\pm \infty \), the \(\omega \)-limit sets may contain functions which are not steady states. Previously, such examples were known for balanced bistable nonlinearities. The novelty of the present result is that it applies to a robust class of nonlinearities. Our proof is based on an analysis of threshold solutions for ordered families of initial data whose limits at infinity are not necessarily zeros of \(f\).     

Whirl Mappings on Generalised Annuli and the Incompressible Symmetric Equilibria of the Dirichlet Energy

In this paper we show a striking contrast in the symmetries of equilibria and extremisers of the total elastic energy of a hyperelastic incompressible annulus subject to pure displacement boundary conditions. Indeed upon considering the equilibrium equations, here, the nonlinear second order elliptic system formulated for the deformation \(u=(u_{1}, \ldots, u_{N})\):

$$ {\mathbb{E}} {\mathbb{L}}[u, {\mathbf {X}}] = \left \{ \textstyle\begin{array}{l@{\quad}l} \Delta u = \operatorname{div}(\mathscr{P} (x) \operatorname{cof} \nabla u) & \textrm{in }{\mathbf {X}},\\ \det\nabla u = 1 & \textrm{in }{\mathbf {X}},\\ u \equiv\varphi& \textrm{on }\partial{\mathbf {X}}, \end{array}\displaystyle \right . $$

where \({\mathbf {X}}\) is a finite, open, symmetric \(N\)-annulus (with \(N \ge2\)), \(\mathscr{P}=\mathscr{P}(x)\) is an unknown hydrostatic pressure field and \(\varphi\) is the identity mapping, we prove that, despite the inherent rotational symmetry in the system, when \(N=3\), the problem possesses no non-trivial symmetric equilibria whereas in sharp contrast, when \(N=2\), the problem possesses an infinite family of symmetric and topologically distinct equilibria. We extend and prove the counterparts of these results in higher dimensions by way of showing that a similar dichotomy persists between all odd vs. even dimensions \(N \ge4\) and discuss a number of closely related issues.

     

Initial and Boundary Blow-U Problem for $$$$-Lalacian Parabolic Equation with General Absortion

In this article, we investigate the initial and boundary blow-up problem for the \(p\)-Laplacian parabolic equation \(u_t-\Delta _p u=-b(x,t)f(u)\) over a smooth bounded domain \(\Omega \) of \(\mathbb {R}^N\) with \(N\ge 2\), where \(\Delta _pu=\mathrm{div}(|\nabla u|^{p-2}\nabla u)\) with \(p>1\), and \(f(u)\) is a function of regular variation at infinity. We study the existence and uniqueness of positive solutions, and their asymptotic behaviors near the parabolic boundary.     

Time Periodic Traveling Curved Fronts in the Periodic Lotka–Volterra Competition–Diffusion System

This paper is concerned with time periodic traveling curved fronts for periodic Lotka–Volterra competition system with diffusion in two dimensional spatial space

$$\begin{aligned} {\left\{ \begin{array}{ll} \dfrac{\partial u_{1}}{\partial t}=\Delta u_{1} +u_{1}(x,y,t)\left( r_{1}(t)-a_{1}(t)u_{1}(x,y,t)-b_{1}(t)u_{2}(x,y,t)\right) ,\\ \dfrac{\partial u_{2}}{\partial t}=d\Delta u_{2} +u_{2}(x,y,t)\left( r_{2}(t)-a_{2}(t)u_{1}(x,y,t)-b_{2}(t)u_{2}(x,y,t)\right) , \end{array}\right. } \end{aligned}$$

where \(\Delta \) denotes \(\frac{\partial ^{2}}{\partial x^{2} }+ \frac{\partial ^{2}}{\partial y^{2} }\), \(x,y\in {\mathbb {R}}\) and \(d>0\) is a constant, the functions \(r_i(t),a_i(t)\) and \(b_i(t)\) are T-periodic and Hölder continuous. Under suitable assumptions that the corresponding kinetic system admits two stable periodic solutions (p(t), 0) and (0, q(t)), the existence, uniqueness and stability of one-dimensional traveling wave solution \(\left( \Phi _{1}(x+ct,t),\Phi _{2}(x+ct,t)\right) \) connecting two periodic solutions (p(t), 0) and (0, q(t)) have been established by Bao and Wang ( J Differ Equ 255:2402–2435, 2013) recently. In this paper we continue to investigate two-dimensional traveling wave solutions of the above system under the same assumptions. First, we establish the asymptotic behaviors of one-dimensional traveling wave solutions of the system at infinity. Using these asymptotic behaviors, we then construct appropriate super- and subsolutions and prove the existence and non-existence of two-dimensional time periodic traveling curved fronts. Finally, we show that the time periodic traveling curved front is asymptotically stable.

     

Zero group velocity longitudinal modes in an isotropic cylinder

Zero group velocity(ZGV) modes are studied in an isotropic cylinder. The L(0, 2) mode behaves anomalously for the mater√ials with a value of the bulk velocity ratio, κ,in the range ~(1/2) κ 2.64 and normally otherwise. All higher modes, except the first few, have no ZGV point for all isotropic materials. This is explained analytically by finding the slope of phase velocity dispersion curves of modes first when the phase velocity equals κ and then at their initial state.     

On Reducibility of 1d Wave Equation with Quasiperiodic in Time Potentials

In this paper we use a KAM theorem of Grébert and Thomann (Commun Math Phys 307:383–427, 2011) to prove the reducibility of the 1d wave equation with Dirichlet boundery conditions on \([0,\pi ]\) with a quasi-periodic in time potential under some symmetry assumptions. From Mathieu–Hill operator’s known results (Eastham in The spectral theory of periodic differential operators, Hafner, New York, 1974; Magnus and Winkler in Hill’s equation, Wiley-Interscience, London, 1969) and Bourgain’s techniques (Commun Math Phys 204:207–247, 1999), we prove that for any \(\epsilon \) small enough, there exist a \(0<m_{\epsilon }\le 1\) and one solution \(u_{\epsilon }(t,x)\) with

$$\begin{aligned} \Vert u_{\epsilon }(t_n,x)\Vert _{H^1({\mathbb {T}})}\rightarrow \infty , \qquad |t_n|\rightarrow \infty , \end{aligned}$$

where \(u_{\epsilon }(t,x)\) satisfies 1d wave equation

$$\begin{aligned} u_{tt}-u_{xx}+m_{\epsilon }u-\epsilon \cos 2t u=0, \end{aligned}$$

with Dirichlet boundery conditions on \([0,\pi ]\).

     

Insight into Heterogeneity Effects in Methane Hydrate Dissociation via Pore-Scale Modeling

A large number (1253) of high-quality streaming potential coefficient (\(C_\mathrm{sp})\) measurements have been carried out on Berea, Boise, Fontainebleau, and Lochaline sandstones (the latter two including both detrital and authigenic overgrowth forms), as a function of pore fluid salinity (\(C_\mathrm{f})\) and rock microstructure. All samples were saturated with fully equilibrated aqueous solutions of NaCl (10\(^{-5}\) and 4.5 mol/dm\(^{3})\) upon which accurate measurements of their electrical conductivity and pH were taken. These \(C_\mathrm{sp}\) measurements represent about a fivefold increase in streaming potential data available in the literature, are consistent with the pre-existing 266 measurements, and have lower experimental uncertainties. The \(C_\mathrm{sp}\) measurements follow a pH-sensitive power law behaviour with respect to \(C_\mathrm{f}\) at medium salinities (\(C_\mathrm{sp} =-\,1.44\times 10^{-9} C_\mathrm{f}^{-\,1.127} \), units: V/Pa and mol/dm\(^{3})\) and show the effect of rock microstructure on the low salinity \(C_\mathrm{sp}\) clearly, producing a smaller decrease in \(C_\mathrm{sp}\) per decade reduction in \(C_\mathrm{f}\) for samples with (i) lower porosity, (ii) larger cementation exponents, (iii) smaller grain sizes (and hence pore and pore throat sizes), and (iv) larger surface conduction. The \(C_\mathrm{sp}\) measurements include 313 made at \(C_\mathrm{f} > 1\) mol/dm\(^{3}\), which confirm the limiting high salinity \(C_\mathrm{sp}\) behaviour noted by Vinogradov et al., which has been ascribed to the attainment of maximum charge density in the electrical double layer occurring when the Debye length approximates to the size of the hydrated metal ion. The zeta potential (\(\zeta \)) was calculated from each \(C_\mathrm{sp}\) measurement. It was found that \(\zeta \) is highly sensitive to pH but not sensitive to rock microstructure. It exhibits a pH-dependent logarithmic behaviour with respect to \(C_\mathrm{f}\) at low to medium salinities (\(\zeta =0.01133 \log _{10} \left( {C_\mathrm{f} } \right) +0.003505\), units: V and mol/dm\(^{3})\) and a limiting zeta potential (zeta potential offset) at high salinities of \({\zeta }_\mathrm{o} = -\,17.36\pm 5.11\) mV in the pH range 6–8, which is also pH dependent. The sensitivity of both \(C_\mathrm{sp}\) and \(\zeta \) to pH and of \(C_\mathrm{sp}\) to rock microstructure indicates that \(C_\mathrm{sp}\) and \(\zeta \) measurements can only be interpreted together with accurate and equilibrated measurements of pore fluid conductivity and pH and supporting microstructural and surface conduction measurements for each sample.