Dynamic Depletion of Vortex Stretching and Non-Blowup of the 3-D Incompressible Euler Equations

We study the interplay between the local geometric properties and the non-blowup of the 3D incompressible Euler equations. We consider the interaction of two perturbed antiparallel vortex tubes using Kerr's initial condition . We use a pseudo-spectral method with resolution up to 1536 × 1024 × 3072 to resolve the nearly singular behavior of the Euler equations. Our numerical results demonstrate that the maximum vorticity does not grow faster than doubly exponential in time, up to t = 19, beyond the singularity time t = 18.7 predicted by Kerr's computations , . The velocity, the enstrophy, and the enstrophy production rate remain bounded throughout the computations. As the flow evolves, the vortex tubes are flattened severely and turned into thin vortex sheets, which roll up subsequently. The vortex lines near the region of the maximum vorticity are relatively straight. This local geometric regularity of vortex lines seems to be responsible for the dynamic depletion of vortex stretching.     

Lagrangian chaos and multiphase processes in vortex flows

We discuss experimental and numerical studies of the effects of Lagrangian chaos (chaotic advection) on the stretching of a drop of an immiscible impurity in a flow. We argue that the standard capillary number used to describe this process is inadequate since it does not account for advection of a drop between regions of the flow with varying velocity gradient. Consequently, we propose a Lagrangian-generalized capillary number C_L number based on finite-time Lyapunov exponents. We present preliminary tests of this formalism for the stretching of a single drop of oil in an oscillating vortex flow, which has been shown previously to exhibit Lagrangian chaos. Probability distribution functions (PDFs) of the stretching of this drop have features that are similar to PDFs of C_L. We also discuss on-going experiments that we have begun on drop stretching in a blinking vortex flow.     

Finite-time Collapse of Three Point Vortices in the Plane

We investigate the finite-time collapse of three point vortices in the plane utilizing the geometric formulation of three-vortexmotion from Krishnamurthy, Aref and Stremler (2018) Phys. Rev. Fluids 3, 024702. In this approach, the vortex system is described in terms of the interior angles of the triangle joining the vortices, the circle that circumscribes that triangle, and the orientation of the triangle. Symmetries in the governing geometric equations of motion for the general three-vortex problem allow us to consider a reduced parameter space in the relative vortex strengths. The well-known conditions for three-vortex collapse are reproduced in this formulation, and we show that these conditions are necessary and sufficient for the vortex motion to consist of collapsing or expanding self-similar motion. The geometric formulation enables a new perspective on the details of this motion. Relationships are determined between the interior angles of the triangle, the vortex strength ratios, the (finite) system energy, the time of collapse, and the distance traveled by the configuration prior to collapse. Several illustrative examples of both collapsing and expanding motion are given.     

Finite-time singularity formation at a magnetic neutral line in Hall magnetohydrodynamics

The formation of a current sheet in a weakly collisional plasma can be modelled as a finite-time singularity solution of magnetohydrodynamic equations. We use an exact self-similar solution to confirm and generalise a previous finding that, in sharp contrast to two-dimensional solutions in standard MHD, a finite-time collapse to a current sheet can occur in Hall MHD. We derive a criterion for the finite-time singularity in terms of initial conditions, and we use an intermediate asymptotic solution for the evolution of an axial magnetic field to obtain a general expression for the singularity formation time. We illustrate the analytical results by numerical solutions.     

Geometric depletion of vortex stretch in 3D viscous incompressible flow

New geometric constraints on vorticity are obtained which suppress possible development of finite-time singularity from the nonlinear vortex stretching mechanism. We find a new condition on the smoothness of the direction of vorticity in the vortical region which yields regularity. We also detect a regularity condition of isotropy type on vorticity in the intensive vorticity region via a new cancellation principle. This is in contrast with the one of isotropy type on the curl of vorticity obtained recently by A. Ruzmaikina and Z. Gruji? [A. Ruzmaikina, Z. Gruji?, On depletion of the vortex-stretching term in the 3D Navier-Stokes equations, Comm. Math. Phys. 247 (2004) 601-611]. We improve as well all of their results by eliminating their assumption that the initial vorticity ω₀ is required to be in L¹.     

The Effect of Surface Tension on the Moore Singularity of Vortex Sheet Dynamics

We investigate the regularization of Moore’s singularities by surface tension in the evolution of vortex sheets and its dependence on the Weber number (which is inversely proportional to surface tension coefficient). The curvature of the vortex sheet, instead of blowing up at finite time t ₀, grows exponentially fast up to a O(We) limiting value close to t ₀. We describe the analytic structure of the solutions and their self-similar features and characteristic scales in terms of the Weber number in a O(We⁻¹) neighborhood of the time at which, in absence of surface tension effects, Moore’s singularity would appear. Our arguments rely on asymptotic techniques and are supported by full numerical simulations of the PDEs describing the evolution of vortex sheets.     

Unobstructed Lagrangian deformations

We prove that deformations of a Lagrangian singularity are unobstructed if the usual (flat) deformations are unobstructed and if a cohomological vanishing condition is satisfied. This gives another application to deformation theory of the Lagrangian de Rham complex introduced in Sevenheck and van Straten (Math. Ann. 327 (1) (2003) 79–102). To prove our theorem, we use the T¹-lifting criterion due to Ran, Kawamata and others. To cite this article: C. Sevenheck, C. R. Acad. Sci. Paris, Ser. I 338 (2004).     

Geometric conditions for Kuhn–Tucker sufficiency of global optimality in mathematical programming

We present geometric criteria for a feasible point that satisfies the Kuhn–Tucker conditions to be a global minimizer of mathematical programming problems with or without bounds on the variables. The criteria apply to multi-extremal programming problems which may have several local minimizers that are not global. We establish such criteria in terms of underestimators of the Lagrangian of the problem. The underestimators are required to satisfy certain geometric property such as the convexity (or a generalized convexity) property. We show that the biconjugate of the Lagrangian can be chosen as a convex underestimator whenever the biconjugate coincides with the Lagrangian at a point. We also show how suitable underestimators can be constructed for the Lagrangian in the case where the problem has bounds on the variables. Examples are given to illustrate our results.     

Resilient and robust finite-time H∞ control for uncertain discrete-time jump nonlinear systems

This paper studies the robust and resilient finite-time H_∞ control problem for uncertain discrete-time nonlinear systems with Markovian jump parameters. With the help of linear matrix inequalities and stochastic analysis techniques, the criteria concerning stochastic finite-time boundedness and stochastic H_∞ finite-time boundedness are initially established for the nonlinear stochastic model. We then turn to stochastic finite-time controller analysis and design to guarantee that the stochastic model is stochastically H_∞ finite-time bounded by employing matrix decomposition method. Applying resilient control schemes, the resilient and robust finite-time controllers are further designed to ensure stochastic H_∞ finite-time boundedness of the derived stochastic nonlinear systems. Moreover, the results concerning stochastic finite-time stability and stochastic finite-time boundedness are addressed. All derived criteria are expressed in terms of linear matrix inequalities, which can be solved by utilizing the available convex optimal method. Finally, the validity of obtained methods is illustrated by numerical examples.     

Finite-time uniform stability of functional differential equations with applications in network synchronization control

In this paper, we investigate finite-time uniform stability of functional differential equations with applications in network synchronization control. First, a Razumikhin-type theorem is derived to ensure finite-time uniform stability of functional differential equations. Based on the theoretical results, finite-time uniform synchronization is proposed for a class of delayed neural networks and delayed complex dynamical networks by designing nontrivial and simple control strategies and some novel criteria are established. Especially, a feasible region of the control parameters for each neuron is derived for the realization of finite-time uniform synchronization of the addressed neural networks, which provide a great convenience for the application of the theoretical results. Finally, two numerical examples with numerical simulations are provided to show the effectiveness and feasibility of the theoretical results.     

Uniqueness of unbounded solutions of the Lagrangian mean curvature flow equation for graphs

We observe that the comparison result of Barles–Biton–Ley for viscosity solutions of a class of nonlinear parabolic equations can be applied to a geometric fully nonlinear parabolic equation which arises from the graphic solutions for the Lagrangian mean curvature flow. To cite this article: J. Chen, C. Pang, C. R. Acad. Sci. Paris, Ser. I 347 (2009).     

Singularity formation during relaxation of jets and fronts toward the state of geostrophic equilibrium

Relaxation of rectilinear jets and fronts towards the state of geostrophic equilibrium (geostrophic adjustment) is studied in three models of increasing complexity. In the one-dimensional rotating shallow water model we perform high-resolution numerical simulations of adjustment and show that wave-breaking and shock formation in the jet core are ubiquitous. We show that adjusted final states are attained in the cases where their existence is not guaranteed by the theory.In the two-layer rotating shallow water model, we show how the baroclinic effects change the scenario of adjustment. We establish criteria of existence and uniqueness of the adjusted states and demonstrate that the adjustment process can lead to a state where the jet traps the fast oscillations. Symmetric instability is shown to appear for strong enough jets.In the full continuously stratified primitive equations model we analyze the final states of the relaxation process by using the Lagrangian variables. We demonstrate that existence of the adjusted state is further restricted as compared to the previous models and show that singularity formation on route to the geostrophic equilibrium is generic.     

Singularities of Lagrangian mean curvature flow: zero-Maslov class case

We study singularities of Lagrangian mean curvature flow in ℂ ⁿ when the initial condition is a zero-Maslov class Lagrangian. We start by showing that, in this setting, singularities are unavoidable. More precisely, we construct Lagrangians with arbitrarily small Lagrangian angle and Lagrangians which are Hamiltonian isotopic to a plane that, nevertheless, develop finite time singularities under mean curvature flow. We then prove two theorems regarding the tangent flow at a singularity when the initial condition is a zero-Maslov class Lagrangian. The first one (Theorem A) states that that the rescaled flow at a singularity converges weakly to a finite union of area-minimizing Lagrangian cones. The second theorem (Theorem B) states that, under the additional assumptions that the initial condition is an almost-calibrated and rational Lagrangian, connected components of the rescaled flow converges to a single area-minimizing Lagrangian cone, as opposed to a possible non-area-minimizing union of area-minimizing Lagrangian cones. The latter condition is dense for Lagrangians with finitely generated H ₁(L,ℤ).     

On the evolution of sharp fronts for the quasi‐geostrophic equation

We consider the problem of the evolution of sharp fronts for the surface quasi‐geostrophic (QG) equation. This problem is the analogue to the vortex patch problem for the two‐dimensional Euler equation. The special interest of the quasi‐geostrophic equation lies in its strong similarities with the three‐dimensional Euler equation, while being a two‐dimen‐sional model. In particular, an analogue of the problem considered here, the evolution of sharp fronts for QG, is the evolution of a vortex line for the three‐dimensional Euler equation. The rigorous derivation of an equation for the evolution of a vortex line is still an open problem. The influence of the singularity appearing in the velocity when using the Biot‐Savart law still needs to be understood. We present two derivations for the evolution of a periodic sharp front. The first one, heuristic, shows the presence of a logarithmic singularity in the velocity, while the second, making use of weak solutions, obtains a rigorous equation for the evolution explaining the influence of that term in the evolution of the curve. Finally, using a Nash‐Moser argument as the main tool, we obtain local existence and uniqueness of a solution for the derived equation in the C^∞ case. © 2004 Wiley Periodicals, Inc.     

A Geometric Heat-Flow Theory of Lagrangian Coherent Structures

We consider Lagrangian coherent structures (LCSs) as the boundaries of material subsets whose advective evolution is metastable under weak diffusion. For their detection, we first transform the Eulerian advection–diffusion equation to Lagrangian coordinates, in which it takes the form of a time-dependent diffusion or heat equation. By this coordinate transformation, the reversible effects of advection are separated from the irreversible joint effects of advection and diffusion. In this framework, LCSs express themselves as (boundaries of) metastable sets under the Lagrangian diffusion process. In the case of spatially homogeneous isotropic diffusion, averaging the time-dependent family of Lagrangian diffusion operators yields Froyland’s dynamic Laplacian. In the associated geometric heat equation, the distribution of heat is governed by the dynamically induced intrinsic geometry on the material manifold, to which we refer as the geometry of mixing. We study and visualize this geometry in detail, and discuss connections between geometric features and LCSs viewed as diffusion barriers in two numerical examples. Our approach facilitates the discovery of connections between some prominent methods for coherent structure detection: the dynamic isoperimetry methodology, the variational geometric approaches to elliptic LCSs, a class of graph Laplacian-based methods and the effective diffusivity framework used in physical oceanography.

     

Reuse,Recycle, Reduce (3R) – strategies for the calculation of transient magnetic fields

The discretization of transient magneto-dynamic field problems with geometric discretization schemes such as the Finite Integration Technique or the Finite-Element Method based on Whitney form functions results in nonlinear differential-algebraic systems of equations of index 1. Their time integration with embedded s-stage singly diagonal implicit Runge–Kutta methods requires the solution of s nonlinear systems within one time step. Accelerated solution of these schemes is achieved with techniques following so-called 3R-strategies (“reuse, recycle, reduce”). This involves e.g. the solution of the linear(-ized) equations in each time step where the solution process of the iterative preconditioned conjugate gradient method reuses and recycles spectral information of linear systems from previous stages. Additionally, a combination of an error controlled spatial adaptivity and an error controlled implicit Runge–Kutta scheme is used to reduce the number of unknowns for the algebraic problems effectively and to avoid unnecessary fine grid resolutions both in space and time. First numerical results for 2D nonlinear magneto-dynamic problems validate the presented approach and its implementation. The space discretization in the numerical examples is done by Lagrangian nodal finite elements but the presented algorithms also work in combination with other discretization schemes for the Maxwell equations such as the Whitney vector finite elements.     

Geometric measure-type regularity criteria for the 3D magnetohydrodynamical system

Several formulations of a local geometric measure-type condition are imposed on super-level sets of mild solutions to the homogeneous incompressible 3D magnetohydrodynamical system with bounded initial data to prevent finite-time singularity formation. Supporting this, results regarding the existence, uniqueness, and real analyticity of mild solutions are established as is a sharp lower bound on the radius of analyticity.     

Numerical stability and numerical viscosity in certain meshless vortex methods as applied to the Navier-Stokes and heat equations

The stability and numerical viscosity of certain Lagrangian methods for flow simulation are analyzed in numerical experiments. Most attention is given to the viscous vortex domains method and its generalization, the viscous vortex thermal domains method.