Tetrahedral frame fields have applications to certain classes of nematic liquid crystals and frustrated media. We consider the problem of constructing a tetrahedral frame field in three dimensional domains in which the boundary normal vector is included in the frame on the boundary. To do this we identify an isomorphism between a given tetrahedral frame and a symmetric, traceless third order tensor under a particular nonlinear constraint. We then define a Ginzburg-Landau-type functional which penalizes the associated nonlinear constraint. Using gradient descent, one retrieves a globally defined limiting tensor outside of a singular set. The tetrahedral frame can then be recovered from this tensor by a determinant maximization method, developed in this work. The resulting numerically generated frame fields are smooth outside of one dimensional filaments that join together at triple junctions.

Fluid flows containing dilute or dense suspensions of thin fibers are widespread in biological and industrial processes. To describe the motion of a thin immersed fiber, or to describe the forces acting on it, it is convenient to work with one-dimensional fiber centerlines and force densities rather than two-dimensional surfaces and surface tractions. Slender body theories offer ways to model and simulate the motion of immersed fibers using only one-dimensional data. However, standard formulations can break down when the fiber surface comes close to intersecting itself or other fibers. In this paper we introduce a numerical method for a recently derived three-dimensional slender body boundary value problem that can be stated entirely in terms of a one-dimensional distribution of forces on the centerline. The method is based on a new completed single-layer potential formulation of fluid velocity which removes the nullspace associated with the unmodified single layer potential. We discretize the model and present numerical results demonstrating the good conditioning and improved performance of the method in the presence of near-intersections. To avoid the modeling and numerical choices involved with free ends, we consider closed fibers.

An $n$-cross field is a locally-defined orthogonal coordinate system invariant with respect to the cubic symmetry group. Cross fields are finding wide-spread use in mesh generation, computer graphics, and materials science among many applications. It was recently by other authors that $3$-cross fields can be embedded into the set of symmetric $4$th-order tensors. Another concurrent work further develops a relaxation of this tensor field via a certain set of varieties. In this paper, we consider the problem of generating an arbitrary $n$-cross field using a fourth-order $Q$-tensor theory that is constructed out of tensored projection matrices. We establish that by a Ginzburg-Landau relaxation towards a global projection, one can reliably generate an $n$-cross field on arbitrary Lipschitz domains. Our work provides a rigorous approach that offers several new results including porting the tensor framework to arbitrary dimensions, providing a new relaxation method that embeds the problem into a global steepest descent, and offering a relaxation scheme for aligning the cross field with the boundary. Our approach is designed to fit within the classical Ginzburg-Landau PDE theory, offering a concrete road map for the future careful study of singularities of energy minimizers.

In this note, we prove the profile decomposition for hyperbolic Schr\odinger (or mixed signature) equations on $\mathbb{R}^2$ in two cases, one mass-supercritical and one mass-critical. First, as a warm up, we show that the profile decomposition works for the ${\dot H}^{\frac12}$ critical problem, which gives a simple generalization of for instance one of the results in Fanelli-Visciglia (2013). Then, we give the derivation of the profile decomposition in the mass-critical case by proving an improved Strichartz estimate. We will use a very similar approach to that laid out in the notes of Killip-Visan (2008), but we are forced to do a double Whitney decomposition to accommodate an extra scaling symmetry that arises in the problem with mixed signature.

Slender body theory is a commonly used approximation in computational models of thin fibers in viscous fluids, especially in simulating the motion of cilia or flagella in swimming microorganisms. In Mori et al. (Commun Pure Appl Math, 2018. arXiv:1807.00178), we developed a PDE framework for analyzing the error introduced by the slender body approximation for closed-loop fibers with constant radius $$\varepsilon $$, and showed that the difference between our closed-loop PDE solution and the slender body approximation is bounded by an expression proportional to $$\varepsilon |\log \varepsilon |$$. Here we extend the slender body PDE framework to the free endpoint setting, which is more physically relevant from a modeling standpoint but more technically demanding than the closed loop analysis. The main new difficulties arising in the free endpoint setting are defining the endpoint geometry, identifying the extent of the 1D slender body force density, and determining how the well-posedness constants depend on the non-constant fiber radius. Given a slender fiber satisfying certain geometric constraints at the filament endpoints and a one-dimensional force density satisfying an endpoint decay condition, we show a bound for the difference between the solution to the slender body PDE and the slender body approximation in the free endpoint setting. The bound is a sum of the same $$\varepsilon |\log \varepsilon |$$ term appearing in the closed loop setting and an endpoint term proportional to $$\varepsilon $$, where $$\varepsilon $$ is now the maximum fiber radius.

We study quench detection in superconducting accelerator cavities cooled with He-II. A rigorous mathematical formula is derived to localize the quench position from dynamical data over a finite time interval at a second sound detector.

We consider the problem of a one-dimensional elastic filament immersed in a two-dimensional steady Stokes fluid. Immersed boundary problems in which a thin elastic structure interacts with a surrounding fluid are prevalent in science and engineering, a class of problems for which Peskin has made pioneering contributions. Using boundary integrals, we first reduce the fluid equations to an evolution equation solely for the immersed filament configuration. We then establish local well-posedness for this equation with initial data in low-regularity Holder spaces. This is accomplished by first extracting the principal linear evolution by a small-scale decomposition and then establishing precise smoothing estimates on the nonlinear remainder. Higher regularity of these solutions is established via commutator estimates with error terms generated by an explicit class of integral kernels. Furthermore, we show that the set of equilibria consists of uniformly parametrized circles and prove nonlinear stability of these equilibria with explicit exponential decay estimates, the optimality of which we verify numerically. Finally, we identify a quantity that respects the symmetries of the problem and controls global-in-time behavior of the system. (c) 2018 Wiley Periodicals, Inc.

We consider the problem of reconstructing the features of a weak anisotropic background potential by the trajectories of vortex dipoles in a nonlinear Gross-Pitaevskii equation. At leading order, the dynamics of vortex dipoles are given by a Hamiltonian system. If the background potential is sufficiently smooth and flat, the background can be reconstructed using ideas from the boundary and the lens rigidity problems. We prove that reconstructions are unique, derive an approximate reconstruction formula, and present numerical examples.

We study the dynamics of vortices in an inhomogeneous Gross-Pitaevskii equation $i \partial_t u = \Delta u + {1\over \varepsilon^2} (p_\varepsilon^2(x) - |u|^2)$. For a unique scaling regime $|p_\varepsilon(x) - 1 | = O(|\log \varepsilon|^{-1})$, it is shown that vortices can interact both with the background perturbation and with each other. Results for associated parabolic and elliptic problems are discussed.

This paper studies questions related to the dynamic transition between local and global minimizers in the Ginzburg–Landau theory of superconductivity. We derive a heuristic equation governing the dynamics of vortices that are close to the boundary, and of dipoles with small inter-vortex separation. We consider a small random perturbation of this equation and study the asymptotic regime under which vortices nucleate.

We examine the Petviashvilli method for solving the equation $ \phi - \Delta \phi = |\phi|^{p-1} \phi$ on a bounded domain $\Omega \subset \mathbb{R}^d$ with Dirichlet boundary conditions. We prove a local convergence result, using spectral analysis, akin to the result for the problem on $\mathbb{R}$ by Pelinovsky & Stepanyants, 2004. We also prove a global convergence result by generating a suite of nonlinear inequalities for the iteration sequence, and we show that the sequence has a natural energy that decreases along the sequence.

We study dynamics of vortices in solutions of the Gross-Pitaevskii equation $i \partial_t u = \Delta u + \varepsilon^{-2} u (1 - |u|^2)$ on $\mathbb{R}^2$ with nonzero degree at infinity. We prove that vortices move according to the classical Kirchhoff-Onsager ODE for a small but finite coupling parameter $\varepsilon$. By carefully tracking errors we allow for asymptotically large numbers of vortices, and this lets us connect the Gross-Pitaevskii equation on the plane to two dimensional incompressible Euler equations through the work of Schochet [21].

This paper is devoted to the problem of recovering a potential q in a domain in ℝ^d for d ≥ 3 from the Dirichlet to Neumann map. This problem is related to the inverse Calderón conductivity problem via the Liouville transformation. It is known from the work of Haberman and Tataru [11] and Nachman and Lavine [17] that uniqueness holds for the class of conductivities of one derivative and the class of W^2,d/2 conductivities respectively. The proof of Haberman and Tataru is based on the construction of complex geometrical optics (CGO) solutions initially suggested by Sylvester and Uhlmann [22], in functional spaces introduced by Bourgain [2]. The proof of the second result, in the work of Ferreira et al. [10], is based on the construction of CGO solutions via Carleman estimates. The main goal of the paper is to understand whether or not an approach which is based on the construction of CGO solutions in the spirit of Sylvester and Uhlmann and involves only standard Sobolev spaces can be used to obtain these results. In fact, we are able to obtain a new proof of uniqueness for the Calderón problem for 1) a slightly different class as the one in [11], and for 2) the class of W^2,d/2 conductivities. The proof of statement 1) is based on a new estimate for CGO solutions and some averaging estimates in the same spirit as in [11]. The proof of statement 2) is on the one hand based on a generalized Sobolev inequality due to Kenig et al. [14] and on another hand, only involves standard estimates for CGO solutions [22]. We are also able to prove the uniqueness of a potential for 3) the class of W^s, 3/s (⫌ W^2, 3/2) conductivities with 3/2 < s < 2 in three dimensions. As far as we know, statement 3) is new.

We study the Landau–Lifshitz–Gilbert equation for the dynamics of a magnetic vortex system. We present a PDE-based method for proving vortex dynamics that does not rely on strong well-preparedness of the initial data and allows for instantaneous changes in the strength of the gyrovector force due to bubbling events. The main tools are estimates of the Hodge decomposition of the supercurrent and an analysis of the defect measure of weak convergence of the stress energy tensor. Ginzburg–Landau equations with mixed dynamics in the presence of excess energy are also discussed.

We establish vortex dynamics for the time-dependent Ginzburg-Landau equation for asymptotically large numbers of vortices for the problem without a gauge field and either Dirichlet or Neumann boundary conditions. As our main tool, we establish quantitative bounds on several fundamental quantities, including the kinetic energy, that lead to explicit convergence rates. For dilute vortex liquids we prove that sequences of solutions converge to the hydrodynamic limit.

We prove via explicitly constructed initial data that solutions to the gravity-capillary wave system in R3 representing a 2d air–water interface immediately fail to be C3 with respect to the initial data if the initial data (h0,ψ0)∈Hs+12⊗Hs for s<3. Similar results hold in R2 domains with a 1d interface. Furthermore, we discuss the related threshold for the pure gravity water wave system.

We consider the linearization of the three-dimensional water waves equation with surface tension about a flat interface. Using oscillatory integral methods, we prove that solutions of this equation demonstrate dispersive decay at the somewhat surprising rate of t(-5/6). This rate is due to competition between surface tension and gravitation at O(1) wave numbers and is connected to the fact that, in the presence of surface tension, there is a so-called "slowest wave". Additionally, we combine our dispersive estimates with L-2 type energy bounds to prove a family of Strichartz estimates.

We study condensate solutions of a nonlinear elliptic equation in ℝ2, which models a W-boson with a cosmic string background. The existence of condensate solutions and an energy identity are discussed, based on which the refined asymptotic behavior of condensate solutions is established by studying the corresponding evolution dynamical system. Applying the “shrinking-sphere” method, we also prove the symmetry under inversions of condensate solutions for some special cases.

A simplified model for the energy of the magnetization of a thin ferromagnetic film gives rise to a version of the theory of Ginzburg–Landau vortices for sphere-valued maps. In particular, we have the development of vortices as a certain parameter tends to 0. The dynamics of the magnetization are ruled by the Landau–Lifshitz–Gilbert equation, which combines characteristic properties of a nonlinear Schrödinger equation and a gradient flow. This paper studies the motion of the vortex centers under this evolution equation.