Non-classical continuum theory for fluids incorporating internal and Cosserat rotation rates

This paper presents a non-classical continuum theory for fluent continua in which the conservation and balance laws are derived by incorporating both internal rotation rates arising from the velocity gradient tensor and the rotation rates of the Cosserats. Specifically, in this non-classical continuum theory we have (1) the usual velocities (\(\bar{ \pmb {\varvec{v }}}\)), (2) the three internal rotation rates (\({}_i^t\bar{ \pmb {\varvec{\Theta }}}\)) about the axes of a fixed triad whose axes are parallel to the x-frame arising from the velocity gradient tensor \((\bar{ \pmb {\varvec{L }}})\) that are completely defined by the antisymmetric part of the velocity gradient tensor, and (3) three additional rotation rates (\({}_e^t\bar{ \pmb {\varvec{\Theta }}}\)) about the axes of the same triad located at each material point as additional three unknown degrees of freedom, referred to as Cosserat rotation rates. This gives rise to \(\bar{ \pmb {\varvec{v }}}\) and \({}_e^t\bar{ \pmb {\varvec{\Theta }}}\) as six degrees of freedom at a material point. The internal rotation rates \({}_i^t\bar{ \pmb {\varvec{\Theta }}}\), often neglected in classical fluid mechanics, exist in all deforming fluent continua as these are due to velocity gradient tensor. When the internal rotation rates \({}_i^t\bar{ \pmb {\varvec{\Theta }}}\) are resisted by deforming fluent continua, conjugate moment tensor arises that together with \({}_i^t\bar{ \pmb {\varvec{\Theta }}}\) may result in energy storage and/or dissipation, which must be considered in the conservation and balance laws. The Cosserat rotation rations \({}_e^t\bar{ \pmb {\varvec{\Theta }}}\) also result in conjugate moment tensor that together with \({}_e^t\bar{ \pmb {\varvec{\Theta }}}\) may also result in energy storage and/or dissipation. The main focus of this paper is a consistent derivation of conservation and balance laws for fluent continua that incorporate the aforementioned physics and associated constitutive theories for thermofluids using the conditions resulting from the entropy inequality. The material coefficients derived in the constitutive theories are clearly defined and discussed.     

A framework for synthetic validation of 3D echocardiographic particle image velocimetry

Particle image velocimetry (PIV) has been significantly advanced since its conception in early 1990s. With the advancement of imaging modalities, applications of 2D PIV have far expanded into biology and medicine. One example is echocardiographic particle image velocimetry that is used for in vivo mapping of the flow inside the heart chambers with opaque boundaries. Velocimetry methods can help better understanding the biomechanical problems. The current trend is to develop three-dimensional velocimetry techniques that take advantage of modern medical imaging tools. This study provides a novel framework for validation of velocimetry methods that are inherently three dimensional such as but not limited to those acquired by 3D echocardiography machines. This framework creates 3D synthetic fields based on a known 3D velocity field \({\mathbf{V}}\) and a given 3D brightness field \({\mathbf{B}}\). The method begins with computing the inverse flow \({\mathbf{V}}^{\varvec{*}} \) based on the velocity field \({\mathbf{V}}\). Then the transformation of \({\mathbf{B}}\), imposed by \({\mathbf{V}}\), is calculated using the computed inverse flow according to \({\mathbf{B}}^{\varvec{*}} \left( {\mathbf{x}} \right) = {\mathbf{B}}\left( {{\mathbf{x}} + {\mathbf{V}}^{\varvec{*}} \left( {\mathbf{x}} \right)} \right)\), where x is the coordinates of voxels in \({\mathbf{B}}^{\varvec{*}} \), with a 3D weighted average interpolation, which provides high accuracy, low memory requirement, and low computational time. To check the validity of the framework, we generate pairs of 3D brightness fields by employing Hill’s spherical vortex velocity field. \({\mathbf{B}}\) and the generated \({\mathbf{B}}^{\varvec{*}} \) are then processed by our in-house 3D particle image velocimetry software to obtain the interrelated velocity field. The results indicates that the computed and imposed velocity fields are in agreement.     

On the Global Regularity of the Two-Dimensional Density Patch for Inhomogeneous Incompressible Viscous Flow

Regarding P.-L. Lions’ open question in Oxford Lecture Series in Mathematics and its Applications, Vol. 3 (1996) concerning the propagation of regularity for the density patch, we establish the global existence of solutions to the two-dimensional inhomogeneous incompressible Navier–Stokes system with initial density given by \({(1 - \eta){\bf 1}_{{\Omega}_{0}} + {\bf 1}_{{\Omega}_{0}^{c}}}\) for some small enough constant \({\eta}\) and some \({W^{k+2,p}}\) domain \({\Omega_{0}}\), with initial vorticity belonging to \({L^{1} \cap L^{p}}\) and with appropriate tangential regularities. Furthermore, we prove that the regularity of the domain \({\Omega_0}\) is preserved by time evolution.     

A Simple Method for the Determination of Adsorption Kinetic Parameters Using Recycle Fixed-Bed Adsorber

This study focuses on the development of a novel analysis technique for determining the intraparticle diffusivity \((D_{\mathrm{S}})\) and fluid film mass transfer coefficient \((k_\mathrm{F})\) from a concentration history curve of a recycle fixed-bed reactor without using the linear driving force approximation and empirical equations for the estimation of \(k_\mathrm{F}\). The recycle fixed-bed method requires lesser amounts of working fluid for experiment purposes, which is an advantage over the usual fixed-bed method. Based on the characterization results of the concentration history curves, simulated by rigorous numerical calculations, a novel analysis technique was established. The \(D_{\mathrm{S}}\) value can be determined from the experimentally obtained time at which the concentration of the curve is minimal \((t_{C\mathrm{min}})\). The \(k_\mathrm{F}\) value can be determined from the \(D_{\mathrm{S}}\) value and Biot number (Bi), which can be estimated from the experimental ratio of the maximum concentration to the minimum concentration \((c_{\max }/c_{\min })\). The \(D_{\mathrm{S}}\) and \(k_\mathrm{F}\) values of phenol adsorption on an activated carbon material were determined experimentally using the proposed analysis method. This method enables the determination of reliable adsorption kinetic parameters through a simple and economical experiment. However, appropriate experimental data must be acquired under regulated experimental conditions, especially in the case of fluctuation of the concentration history curve.     

Exact solution of 3D Timoshenko beam problem using linked interpolation of arbitrary order

For arbitrary polynomial loading and a sufficient finite number of nodal points N, the solution for the 3D Timoshenko beam differential equations is polynomial and given as \({{\varvec \theta} = \sum_{i=1}^N I_i {\varvec \theta}_i}\) for the rotation field and \({{\bf u} = \sum_{i=1}^{N+1} J_i {\bf u}_i}\) for the displacement field, where I _i and J _i are the Lagrangian polynomials of order N?1 and N, respectively. It has been demonstrated in this work that the exact solution for the displacement field may be also written in a number of alternative ways involving contributions of the nodal rotations including \({{\bf u} = \sum_{i=1}^N I_i \left[ {\bf u}_i + \frac 1 N ( {\varvec \theta} - {\varvec \theta}_i ) \times {\bf R}_i \right]}\), where R _i are the beam nodal positions.     

Uniaxial experimental study of the acoustic emission and deformation behavior of composite rock based on 3D digital image correlation (DIC)

In this paper, uniaxial compression tests were carried out on a series of composite rock specimens with different dip angles, which were made from two types of rock-like material with different strength. The acoustic emission technique was used to monitor the acoustic signal characteristics of composite rock specimens during the entire loading process. At the same time, an optical non-contact 3 D digital image correlation technique was used to study the evolution of axial strain field and the maximal strain field before and after the peak strength at different stress levels during the loading process. The effect of bedding plane inclination on the deformation and strength during uniaxial loading was analyzed. The methods of solving the elastic constants of hard and weak rock were described. The damage evolution process, deformation and failure mechanism, and failure mode during uniaxial loading were fully determined. The experimental results show that the θ = 0?–45?specimens had obvious plastic deformation during loading, and the brittleness of the θ = 60?–90?specimens gradually increased during the loading process. When the anisotropic angle θincreased from 0?to 90?, the peak strength, peak strain,and apparent elastic modulus all decreased initially and then increased. The failure mode of the composite rock specimen during uniaxial loading can be divided into three categories:tensile fracture across the discontinuities(θ = 0?–30?), slid-ing failure along the discontinuities(θ = 45?–75?), and tensile-split along the discontinuities(θ = 90?). The axial strain of the weak and hard rock layers in the composite rock specimen during the loading process was significantly different from that of the θ = 0?–45?specimens and was almost the same as that of the θ = 60?–90?specimens. As for the strain localization highlighted in the maximum principal strain field, the θ = 0?–30?specimens appeared in the rock matrix approximately parallel to the loading direction,while in the θ = 45?–90?specimens it appeared at the hard and weak rock layer interface.     

Global Well-Posedness for Navier–Stokes Equations with Small Initial Value in $${B^{0}_{n,\infty}(\Omega)}$$

We prove global well-posedness for instationary Navier–Stokes equations with initial data in Besov space \({B^{0}_{n,\infty}(\Omega)}\) in whole and half space, and bounded domains of \({{\mathbb R}^{n}}\), \({n \geq 3}\). To this end, we prove maximal \({L^{\infty}_{\gamma}}\) -regularity of the sectorial operators in some Banach spaces and, in particular, maximal \({L^{\infty}_{\gamma}}\) -regularity of the Stokes operator in little Nikolskii spaces \({b^{s}_{q,\infty}(\Omega)}\), \({s \in (-1, 2)}\), which are of independent significance. Then, based on the maximal regularity results and \({b^{s_{1}}_{q_{1},\infty}-B^{s_{2}}_{q_{2,1}}}\) estimates of the Stokes semigroups, we prove global well-posedness for Navier–Stokes equations under smallness condition on \({\|u_{0}\|_{B^{0}_{n,\infty}(\Omega)}}\) via a fixed point argument using Banach fixed point theorem.     

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\).

     

In Situ Chemically-Selective Monitoring of Multiphase Displacement Processes in a Carbonate Rock Using 3D Magnetic Resonance Imaging

Accurate monitoring of multiphase displacement processes is essential for the development, validation and benchmarking of numerical models used for reservoir simulation and for asset characterization. Here we demonstrate the first application of a chemically-selective 3D magnetic resonance imaging (MRI) technique which provides high-temporal resolution, quantitative, spatially resolved information of oil and water saturations during a dynamic imbibition core flood experiment in an Estaillades carbonate rock. Firstly, the relative saturations of dodecane (\(S_{\mathrm{o}})\) and water (\(S_{\mathrm{w}})\), as determined from the MRI measurements, have been benchmarked against those obtained from nuclear magnetic resonance (NMR) spectroscopy and volumetric analysis of the core flood effluent. Excellent agreement between both the NMR and MRI determinations of \(S_{\mathrm{o}}\) and \(S_{\mathrm{w}}\) was obtained. These values were in agreement to 4 and 9% of the values determined by volumetric analysis, with absolute errors in the measurement of saturation determined by NMR and MRI being 0.04 or less over the range of relative saturations investigated. The chemically-selective 3D MRI method was subsequently applied to monitor the displacement of dodecane in the core plug sample by water under continuous flow conditions at an interstitial velocity of \(1.27\times 10^{-6}\,\hbox {m}\,\hbox {s}^{-1}\) (\(0.4\,\hbox {ft}\,\hbox {day}^{-1})\). During the core flood, independent images of water and oil distributions within the rock core plug at a spatial resolution of \(0.31\,\hbox {mm}\times 0.39\,\hbox {mm} \times 0.39\,\hbox {mm}\) were acquired on a timescale of 16 min per image. Using this technique the spatial and temporal dynamics of the displacement process have been monitored. This MRI technique will provide insights to structure–transport relationships associated with multiphase displacement processes in complex porous materials, such as those encountered in petrophysics research.     

A KAM Theorem for Higher Dimensional Forced Nonlinear Schrödinger Equations

In this paper, we prove an infinite dimensional KAM theorem. As an application, it is shown that there are many real-analytic small-amplitude linearly-stable quasi-periodic solutions for higher dimensional nonlinear Schrödinger equations with outer force

$$\begin{aligned} iu_t-\triangle u +M_\xi u+f(\bar{\omega }t)|u|^2u=0, \quad t\in {{\mathbb R}}, x\in {{\mathbb T}}^d \end{aligned}$$

where \(M_\xi \) is a real Fourier multiplier,\(f({\bar{\theta }})({\bar{\theta }}={\bar{\omega }} t)\) is real analytic and the forced frequencies \(\bar{\omega }\) are fixed Diophantine vectors.

     

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.

     

A Shilnikov Phenomenon Due to State-Dependent Delay,by Means of the Fixed Point Index

The first part of this paper is a general approach towards chaotic dynamics for a continuous map \(f:X\supset M\rightarrow X\) which employs the fixed point index and continuation. The second part deals with the differential equation

$$\begin{aligned} x'(t)=-\alpha \,x(t-d_{{\varDelta }}(x_t)). \end{aligned}$$

with state-dependent delay. For a suitable parameter \(\alpha \) close to \(5\pi /2\) we construct a delay functional \(d_{{\varDelta }}\), constant near the origin, so that the previous equation has a homoclinic solution, \(h(t)\rightarrow 0\) as \(t\rightarrow \pm \infty \), with certain regularity properties of the linearization of the semiflow along the flowline \(t\mapsto h_t\). The third part applies the method from the beginning to a return map which describes solution behaviour close to the homoclinic loop, and yields the existence of chaotic motion.

     

Approximation of Flat <Emphasis Type="Italic">W</Emphasis><Superscript>2,2</Superscript> Isometric Immersions by Smooth Ones

Let \({S\subset\mathbb{R}^2}\) be a bounded Lipschitz domain and denote by \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\) the set of mappings \({u\in W^{2,2}(S;\mathbb{R}^3)}\) which satisfy \({(\nabla u)^T(\nabla u) = Id}\) almost everywhere. Under an additional regularity condition on the boundary \({\partial S}\) (which is satisfied if \({\partial S}\) is piecewise continuously differentiable), we prove that the strong W ^2,2 closure of \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)\cap C^{\infty}(\overline{S};\mathbb{R}^3)}\) agrees with \({W^{2,2}_{\text{iso}}(S; \mathbb{R}^3)}\).     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

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.

     

Measurements of Transitional and Turbulent Flow in a Randomly Packed Bed of Spheres with Particle Image Velocimetry

Particle image velocimetry (PIV) has been used to investigate transitional and turbulent flow in a randomly packed bed of mono-sized transparent spheres at particle Reynolds number, \(20<{{ Re}}_{\mathrm{p}}< 3220\). The refractive index of the liquid is matched with the spheres to provide optical access to the flow within the bed without distortions. Integrated pressure drop data yield that Darcy law is valid at \({{ Re}}_{\mathrm{p}} \approx 80\). The PIV measurements show that the velocity fluctuations increase and that the time-averaged velocity distribution start to change at lower \({{ Re}}_{\mathrm{p}}\). The probability for relatively low and high velocities decreases with \({{ Re}}_{\mathrm{p}}\) and recirculation zones that appear in inertia dominated flows are suppressed by the turbulent flow at higher \({{ Re}}_{\mathrm{p}}\). Hence there is a maximum of recirculation at about \({{ Re}}_{\mathrm{p}} \approx 400\). Finally, statistical analysis of the spatial distribution of time-averaged velocities shows that the velocity distribution is clearly and weakly self-similar with respect to \({{ Re}}_{\mathrm{p}}\) for turbulent and laminar flow, respectively.     

Uniqueness of Positive Ground State Solutions of the Logarithmic Schrödinger Equation

We prove the uniqueness of positive ground state solutions of the problem \({ {\frac {d^{2}u}{dr^{2}}} + {\frac {n-1}{r}}{\frac {du}{dr}} + u \ln(|u|) = 0}\), \({u(r) > 0~\forall r \ge 0}\), and \({(u(r),u'(r)) \to (0, 0)}\) as \({r \to \infty}\). This equation is derived from the logarithmic Schrödinger equation \({{\rm i}\psi_{t} = {\Delta} \psi + u \ln \left(|u|^{2}\right)}\), and also from the classical equation \({{\frac {\partial u}{\partial t}} = {\Delta} u +u \left(|u|^{p-1}\right) -u}\). For each \({n \ge 1}\), a positive ground state solution is \({ u_{0}(r) = \exp \left(-{\frac{r^2}{4}} + {\frac{n}{2}}\right),~0 \le r < \infty}\). We combine \({u_{0}(r)}\) with energy estimates and associated Ricatti equation estimates to prove that, for each \({n \in \left[1, 9 \right]}\), \({u_{0}(r)}\) is the only positive ground state. We also investigate the stability of \({u_{0}(r)}\). Several open problems are stated.     

Further Dense Properties of the Space of Circle Diffeomorphisms with a Liouville Rotation Number

In continuation of Matsumoto’s paper (Nonlinearity 25:1495–1511, 2012) we show that various subspaces are \(C^{\infty }\)-dense in the space of orientation-preserving \(C^{\infty }\)-diffeomorphisms of the circle with rotation number \(\alpha \), where \(\alpha \in {\mathbb {S}}^1\) is any prescribed Liouville number. In particular, for every odometer \({\mathcal {O}}\) of product type we prove the denseness of the subspace of diffeomorphisms which are orbit-equivalent to \({\mathcal {O}}\).     

Resolvent Estimates for High-Contrast Elliptic Problems with Periodic Coefficients

We study the asymptotic behaviour of the resolvents \({(\mathcal{A}^\varepsilon+I)^{-1}}\) of elliptic second-order differential operators \({{\mathcal{A}}^\varepsilon}\) in \({\mathbb{R}^d}\) with periodic rapidly oscillating coefficients, as the period \({\varepsilon}\) goes to zero. The class of operators covered by our analysis includes both the “classical” case of uniformly elliptic families (where the ellipticity constant does not depend on \({\varepsilon}\)) and the “double-porosity” case of coefficients that take contrasting values of order one and of order \({\varepsilon^2}\) in different parts of the period cell. We provide a construction for the leading order term of the “operator asymptotics” of \({(\mathcal{A}^\varepsilon+I)^{-1}}\) in the sense of operator-norm convergence and prove order \({O(\varepsilon)}\) remainder estimates.