Divergence-free Vector Research Articles

Localization of Beltrami fields: Global smooth solutions and vortex reconnection for the Navier-Stokes equations

We introduce a class of divergence-free vector fields on R3 obtained after a suitable localization of Beltrami fields. First, we use them as initial data to construct unique global smooth solutions of the three dimensional Navier-Stokes equations. The relevant fact here is that these initial data can be chosen to be large in any critical space for the Navier–Stokes problem, however they satisfy the nonlinear smallness assumption introduced in [10]. As a further application of the method, we use these vector fields to provide analytical example of vortex-reconnection for the three-dimensional Navier-Stokes equations on R3. To do so, we exploit the ideas developed in [13] but differently from this latter we cannot rely on the non-trivial homotopy of the three-dimensional torus. To overcome this obstacle we use a different topological invariant, i.e. the number of hyperbolic zeros of the vorticity field.

New spinorial mass-quasilocal angular momentum inequality for initial data with marginally future trapped surface

We prove a new geometric inequality that relates the Arnowitt–Deser–Misner mass of initial data to a quasilocal angular momentum of a marginally outer trapped surface (MOTS) inner boundary. The inequality is expressed in terms of a 1-spinor, which satisfies an intrinsic first-order Dirac-type equation. Furthermore, we show that if the initial data is axisymmetric, then the divergence-free vector used to define the quasilocal angular momentum cannot be a Killing field of the generic boundary.

Density of weak solutions of the fractional Navier-Stokes equations in the smooth divergence-free vector fields

Density of weak solutions of the fractional Navier-Stokes equations in the smooth divergence-free vector fields

A decomposition theorem of surface vector fields and spectral structure of the Neumann-Poincaré operator in elasticity

We prove that the space of vector fields on the boundary of a bounded domain with the Lipschitz boundary in three dimensions is decomposed into three subspaces: elements of the first one extend to inside the domain as divergence-free and rotation-free vector fields, the second one to the outside as divergence-free and rotation-free vector fields, and the third one to both the inside and the outside as divergence-free harmonic vector fields. We then show that each subspace in the decomposition is infinite-dimensional. We also prove under a mild regularity assumption on the boundary that the decomposition is almost direct in the sense that any intersection of two subspaces is finite-dimensional. We actually prove that the dimension of intersection is bounded by the first Betti number of the boundary. In particular, if the boundary is simply connected, then the decomposition is direct. We apply this decomposition theorem to investigate spectral properties of the Neumann-Poincaré operator in elasticity, whose cubic polynomial is known to be compact. We prove that each linear factor of the cubic polynomial is compact on each subspace of decomposition separately and those subspaces characterize eigenspaces of the Neumann-Poincaré operator. We then prove all the results for three dimensions, decomposition of surface vector fields and spectral structure, are extended to higher dimensions. We also prove analogous but different results in two dimensions.

A divergence-free vector finite-element method for efficient 3D magnetotelluric forward modeling

For large-scale magnetotelluric (MT) forward modeling using vector finite-element methods, iterative solvers are commonly applied. However, as frequency decreases close to zero, the iterative solution process struggles to converge. This is mainly caused by the fact that the weak conductivity term in the curl-curl equation governing the electromagnetic (EM) diffusion in the earth fades during the iterative solution for EM fields. This can lead to violation of the divergence-free condition for current and false jumps in the calculated fields. We develop a regularization technique for vector finite-element MT forward modeling, in which a scaled grad-div operator for electrical fields is included in the curl-curl equation to enforce the divergence-free condition explicitly. Because of its use of basis functions, the direct edge element discretization of the scaled term can lead to zero coefficients. To address this, an equivalent substitute of the scaled grad-div operator based on potential representation of electrical fields is used, discretized with node elements. For one specific element, this is finally reduced to a scaled gradient of the surface integral of current across all the surfaces of the element, resulting in better connectivity of the linear system of equations. The correctness of our algorithm is verified with two synthetic models and an inversion model from real data. The numerical performance is compared to the traditional iterative divergence correction technique (and without the application of the divergence correction for one case) for each of the three models. The results indicate that our algorithm is generally more efficient and stable compared to the traditional technique for all models at all periods considered, with significant improvement at long periods.

On shifting the principal eigenvalue of Dirichlet problem to infinity with non-transversal incompressible drift

We prove that it is always possible to add some divergence free drift vector field to some two dimensional spherical Dirichlet problem, such that the resulting principal eigenvalue lies above a prescribed bound. By construction those drift vector fields vanish on the boundary and their flow lines individually stay away from the boundary. The capacity of those drift vector fields to accelerate diffusivity originates from high frequency oscillation of the associated flow lines. The lower bounds for the spectrum are obtained through isoperimetric inequalities for flow invariant functions.

Error estimates of divergence-free generalized moving least squares (Div-Free GMLS) derivatives approximations in Sobolev spaces

The divergence-free generalized moving least squares (Div-Free GMLS) approximation has recently been utilized to solve some incompressible fluid flows problems. In our recent work (Mohammadi and Dehghan (2021) [28]), we have presented its formulation more precisely, and also the error estimates of derivatives have been carried out in L∞(Ω), where Ω⊂Rd is a bounded set satisfying an interior cone condition. However, the error estimates of this vector-valued approximation in Sobolev spaces are not done. So, in this paper we make the error estimates of Div-Free GMLS derivatives approximations in Lq(Ω), where 1≤q≤∞, using a stable local divergence-free polynomial reproduction property. Note that, the method is a direct approximants of exact derivatives of a divergence-free vector field, which possesses the optimal rates of convergence. This vector-valued technique can also be developed to find the numerical solution of the incompressible fluid flows problems easier than the other available mesh-dependent methods. Finally, we have shown how the proposed approximation can recover the velocity field variable of the well-known Darcy's problem in a two-dimensional space.

Existence of global symmetries of divergence-free fields with first integrals

The relationship between symmetry fields and first integrals of divergence-free vector fields is explored in three dimensions in light of its relevance to plasma physics and magnetic confinement fusion. A Noether-type theorem is known: for each such symmetry, there corresponds a first integral. The extent to which the converse is true is investigated. In doing so, a reformulation of this Noether-type theorem is found for which the converse holds on what is called the toroidal region. Some consequences of the methods presented are quick proofs of the existence of flux coordinates for magnetic fields in high generality, without needing to assume a symmetry such as in the cases of magneto-hydrostatics or quasi-symmetry.

Rigidity Results On Generalized m-Quasi Einstein Manifolds with Associated Affine Killing Vector Field.

We study a non-trivial generalized $m$-quasi Einstein manifold $M$ with finite $m$ and associated divergence-free affine Killing vector field, and show that $M$ reduces to an $m$-quasi Einstein manifold. In addition, if $M$ is complete, then it splits as the product of a line and an $(n-1)$-dimensional negatively Einstein manifold. Finally, we show that the same result holds for a complete non-trivial $m$-quasi Einstein manifold $M$ with finite $m$ and associated affine Killing vector field.

Exploring Finite-Time Singularities and Onsager’s Conjecture with Endpoint Regularity in the Periodic Navier Stokes Equations

It has recently been proposed by the author of the present work that the periodic NS equations (PNS) with high energy assumption can breakdown in finite time but with sufficient low energy scaling the equations may not exhibit finite time blowup.This article gives a general model using specific periodic special functions, that is degenerate elliptic Weierstrass P functions whose presence in the governing equations through the forcing terms simplify the PNS equations at the centers of cells of the 3-Torus. Satisfying a divergence free vector field and periodic boundary conditions respectively with a general spatio-temporal forcing term which is smooth and spatially periodic, the existence of solutions which blowup in finite time for PNS can occur starting with the first derivative and higher with respect to time. P. Isett (2016) has shown that the conservation of energy fails for the 3D incompressible Euler flows with Hlder regularity below 1/3.(Onsager’s second conjecture) The endpoint regularity in Onsager’s conjecture is addressed, and it is found that conservation of energy occurs when the Hlder regularity is exactly 1/3. The endpoint regularity problem has important connections with turbulence theory. Finally very recent developed new governing equations of fluid mechanics are proposed to have no finite time singularities.

A sharp interface immersed edge-based smoothed finite element method with extended fictitious domain scheme

This paper proposes a versatile and robust immersed edge-based smoothed finite element method with the mass conservation algorithm (IESFEM/Mass) to solve partitioned fluid–structure interaction (FSI). A gradient smoothing technique was used to solve the system governing equations, which can improve the calculated capability of the linear triangular elements in two phases. Based on the quadratic sharp interface representation of immersed boundary, an extended fictitious domain constructed by a least squares method approximately corrected the residual flux error. The compatibility for boundary conditions on moving interfaces was satisfied, thus eliminating spurious oscillations. The results from all numerical examples were consistent with those from the existing experiments and published numerical solutions. Furthermore, the present divergence-free vector field had a faster-converged rate in the flow velocity, pressure, and FSI force. Even if in distorted meshes, the proposed algorithm maintained a stable accuracy improvement. The aerodynamics of one- and two-winged flapping motions in insect flight has been investigated through the IESFEM/Mass. It can be seen that the wing–wake interaction mechanism is a vital factor affecting the lift. The applicability of the present method in the biological FSI scenario was also well-demonstrated.

Lyapunov exponents and entropy for divergence-free Lipschitz vector fields

Let mathfrak {X}^{0,1}_nu (M) be the subset of divergence-free Lipschitz vector fields defined on a closed Riemannian manifold M endowed with the Lipschitz topology Vert ,{cdot },Vert _{0,1} where nu is the volume measure. Let mathfrak {X}^{0,1}_{nu ,ell }(M)subset mathfrak {X}^{0,1}_nu (M) be the subset of vector fields satisfying the ell -property, a property that implies C^1-regularity nu -almost everywhere. We prove that there exists a residual subset with respect to Vert ,{cdot },Vert _{0,1} such that Pesin’s entropy formula holds, i.e. for any the metric entropy equals the integral, with respect to nu , of the sum of the positive Lyapunov exponents.

Adaptive partition of unity interpolation method with moving patches

The adaptive partition of unity interpolation method, introduced by Aiton and Driscoll, using Chebyshev local interpolants, is explored for interpolating functions with sharp gradients representing two-medium problems. For functions that evolve under vector fields, the partition of unity patches (covers) can be shifted and resized to follow the changing dynamics of local profiles. The method is tested for selected 1D and 2D two-medium problems with linear divergence-free vector fields. In those cases, the volume fraction in each patch contributing to volume conservation throughout the domain can be kept in high accuracy down to machine precisions. Applications that could benefit from the method include volume tracking and multiphase flow modeling.

Sharp Uncertainty Principle inequality for solenoidal fields

This paper solves the L2 version of Maz'ya's open problem [1, Section 3.9] on the sharp uncertainty principle inequality∫RN|∇u|2dx∫RN|u|2|x|2dx≥CN(∫RN|u|2dx)2 for solenoidal (namely divergence-free) vector fields u=u(x) on RN. The best value of the constant for N≥3 turns out to be CN=14(N2−4(N−3)+2)2 which exceeds the original value N2/4 for unconstrained fields. Moreover, we show the attainability of CN and specify the profiles of the extremal solenoidal fields: for N≥4, the extremals are proportional to a poloidal field that is axisymmetric and unique up to the axis of symmetry; for N=3, there exist extremal toroidal fields, in addition to extremal poloidal fields; for N=2, the extremal fields are all toroidal.

Locally-symplectic neural networks for learning volume-preserving dynamics

We propose locally-symplectic neural networks LocSympNets for learning the flow of phase volume-preserving dynamics. The construction of LocSympNets stems from the theorem of the local Hamiltonian description of the divergence-free vector field and the splitting methods based on symplectic integrators. Symplectic gradient modules of the recently proposed symplecticity-preserving neural networks SympNets are used to construct invertible locally-symplectic modules, which compositions result in volume-preserving neural networks LocSympNets. To further preserve properties of the flow of a dynamical system LocSympNets are extended to symmetric locally-symplectic neural networks SymLocSympNets, such that the inverse of SymLocSympNets is equal to the feed-forward propagation of SymLocSympNets with the negative time step, which is a general property of the flow of a dynamical system. LocSympNets and SymLocSympNets are studied numerically considering learning linear and nonlinear volume-preserving dynamics. In particular, we demonstrate learning of linear traveling wave solutions to the semi-discretized advection equation, periodic trajectories of the Euler equations of the motion of a free rigid body, and quasi-periodic solutions of the charged particle motion in an electromagnetic field. LocSympNets and SymLocSympNets can learn linear and nonlinear dynamics to a high degree of accuracy even when random noise is added to the training data. In all numerical experiments, SymLocSympNets have produced smaller errors in long-time predictions compared to the LocSympNets. When learning a single trajectory of the rigid body dynamics locally-symplectic neural networks can learn both quadratic invariants of the system with absolute relative errors below 1%. In addition, SymLocSympNets produce qualitatively good long-time predictions, when the learning of the whole system from randomly sampled data is considered. LocSympNets and SymLocSympNets can produce accurate short-time predictions of quasi-periodic solutions, which is illustrated in the example of the charged particle motion in an electromagnetic field.

On a class of interpolation inequalities on the 2D sphere

We prove estimates for the $L^p$-norms of systems of functions and divergence-free vector functions that are orthonormal in the Sobolev space $H^1$ on the 2D sphere. As a corollary, order sharp constants for the embedding $H^1\hookrightarrow L^q$, $q<\infty$, are obtained in the Gagliardo-Nirenberg interpolation inequalities. Bibliography: 25 titles.

Spectral analysis and domain truncation for Maxwell's equations

We analyse how the spectrum of the anisotropic Maxwell system with bounded conductivity σ on a Lipschitz domain Ω is approximated by domain truncation. First we prove a new non-convex enclosure for the spectrum of the Maxwell system, with weak assumptions on the geometry of Ω and none on the behaviour of the coefficients at infinity. We also establish a simple criterion for non-accumulation of eigenvalues at iR as well as resolvent estimates. For asymptotically constant coefficients, we describe the essential spectrum and show that spectral pollution may occur only in the essential numerical range We(L∞)⊂R of the quadratic pencil L∞(ω)=μ∞−1curl2−ω2ε∞, acting on divergence-free vector fields. Further, every isolated spectral point of the Maxwell system lying outside We(L∞) and outside the part of the essential spectrum on iR is approximated by spectral points of the Maxwell system on the truncated domains. Our analysis is based on two new abstract results on the (limiting) essential spectrum of polynomial pencils and triangular block operator matrices, which are of general interest. We believe our strategy of proof could be used to establish domain truncation spectral exactness for more general classes of non-self-adjoint differential operators and systems with non-constant coefficients.

Helicity, linking and the distribution of null-homologous periodic orbits for Anosov flows* *2022. This work is licensed under a CC BY license. Solly Coles was supported by the UK Engineering and Physical Sciences Research Council.

This paper concerns connections between dynamical systems, knots and helicity of vector fields. For a divergence-free vector field on a closed three-manifold that generates an Anosov flow, we show that the helicity of the vector field may be recovered as the limit of appropriately weighted averages of linking numbers of periodic orbits, regarded as knots. This complements a classical result of Arnold and Vogel that, when the manifold is a real homology three-sphere, the helicity may be obtained as the limit of the normalised linking numbers of typical pairs of long trajectories. We also obtain results on the asymptotic distribution of weighted averages of null-homologous periodic orbits.

On the advection-diffusion equation with rough coefficients: Weak solutions and vanishing viscosity

We deal with the vanishing viscosity scheme for the transport/continuity equation ∂tu+div(ub)=0 drifted by a divergence-free vector field b. Under general Sobolev assumptions on b, we show the convergence of such scheme to the unique Lagrangian solution of the transport equation. Our proof is based on the use of stochastic flows and yields quantitative rates of convergence. This offers a completely general selection criterion for the transport equation (even beyond the distributional regime) which compensates the wild non-uniqueness phenomenon for solutions with low integrability arising from convex integration constructions, as shown in recent works [8,28–30], and rules out the possibility of anomalous dissipation.

An efficient construction of divergence-free spaces in the context of exact finite element de Rham sequences

Exact finite element de Rham subcomplexes relate conforming subspaces in H1(Ω), H(curl;Ω), H(div;Ω), and L2(Ω) in a simple way by means of differential operators (gradient, curl, and divergence). The characteristics of such strong couplings are crucial for the design of stable and conservative discretizations of mixed formulations for a variety of multiphysics systems. This work explores these aspects for the construction of divergence-free vector shape functions in a robust fashion allowing stable and faster simulations of mixed formulations of incompressible porous media flows. The resulting schemes are verified by means of numerical tests with known smooth solutions and applied to a benchmark problem to confirm the expected theoretical and computational performance results.