Singular Field Research Articles

Simple, efficient method of calculating the Detweiler-Whiting singular field to very high order

Simple, efficient method of calculating the Detweiler-Whiting singular field to very high order

A frequency-independent bound on trigonometric polynomials of Gaussians and applications

We prove a frequency-independent bound on trigonometric functions of a class of singular Gaussian random fields, which arise naturally from weak universality problems for singular stochastic PDEs. This enables us to reduce the regularity assumption on the nonlinearity of the microscopic models (for pathwise convergence) in KPZ and dynamical Φ34 in the previous works of Hairer-Xu and Furlan-Gubinelli to heuristically optimal thresholds required by PDE structures.

A derivative free projection method for the singularities of vector fields with convex constraints on Hadamard manifolds

ABSTRACT The objective of this paper is to introduce a derivative free projection method designed to find the singularities of pseudomonotone vector fields with convex constraints on Hadamard manifolds. This innovative approach combines the hyperplane projection method with a novel search direction. The global convergence of the proposed method is established under certain conditions. Our method improves some existing results in the literature on Hadamard manifolds. Additionally, illustrative numerical examples are provided to demonstrate the practical efficacy of our method.

Singularity analysis for V‐notches in a piezoelectric multimaterial and functionally graded piezoelectric materials

AbstractThis study develops a novel singularity analysis method that addresses the limitations posed by piezoelectric material quantity and the variation patterns of property functions in functionally graded piezoelectric materials (FGPMs). Given that piezoelectric composite materials consist of multiple materials, the material properties of FGPMs vary with respect to angle. By introducing the asymptotic assumption that governs the physical fields close to the notch vertex into the static equilibrium equations, a comprehensive set of characteristic equations is formulated. All singularity orders and corresponding characteristic angle functions of the notch are determined by numerically solving the established characteristic equations. The effects of the notch opening angle, piezoelectric material polarization direction, and boundary conditions on the electromechanical field singularity at the notch are assessed. The judicious selection of material variation patterns can alleviate notch singularity in FGPMs.

Experimental and numerical estimation of complex stress intensity factor for the completely debonded anti-crack embedded into a weak matrix using domain integral method

The problem of interface debonding is of fundamental importance in understanding the stress transfer mechanism and behavior of weak interfaces present due to rigid stiffeners. A domain integral method based on Betti’s reciprocal theorem was adopted to compute the complex stress intensity factor (SIF) for the completely debonded anti-crack (discontinuity) with mixed boundary conditions. The domain integral approach is based on two admissible mechanical states, which comprise of actual and auxiliary elastic fields. The singular auxiliary fields in terms of angular variations for the completely debonded anticrack were obtained using the full-field solution. The obtained angular variations were compared with the Williams eigenvalue problem for completeness. The displacement and strains were obtained using digital image correlation (DIC) experiments and numerical simulation for parallel and perpendicular orientations of the discontinuity. The obtained DIC strain contours revealed the presence of anti-symmetric and symmetric variation on the bonded side of the anti-crack. The complex SIF is estimated experimentally and numerically based on the domain integral method to seek a better comparison with the analytical estimate. The SIF estimated for the perpendicular orientation is higher than the parallel orientation.

Geometry of transcendental singularities of complex analytic functions and vector fields

Abstract On Riemann surfaces M M , there exists a canonical correspondence between a possibly multivalued function Ψ X {\Psi }_{X} whose differential is single-valued (i.e. an additively automorphic singular complex analytic function) and a vector field X X . From the point of view of vector fields, the singularities that we consider are zeros, poles, isolated essential singularities, and accumulation points of the above. The theory of singularities of the inverse function Ψ X ‒ 1 {\Psi }_{X}^{‒1} is extended from meromorphic functions to additively automorphic singular complex analytic functions. The main contribution is a complete characterization of when a singularity of Ψ X − 1 {\Psi }_{X}^{-1} is algebraic, is logarithmic, or arises from a zero with non-zero residue of X X . Relationships between analytical properties of Ψ X {\Psi }_{X} , singularities of Ψ X − 1 {\Psi }_{X}^{-1} and singularities of X X are presented. Families and sporadic examples showing the geometrical richness of vector fields on the neighbourhoods of the singularities of Ψ X − 1 {\Psi }_{X}^{-1} are studied. As applications, we have; a description of the maximal univalence regions for complex trajectory solutions of a vector field X X , a geometric characterization of the incomplete real trajectories of a vector field X X , and a description of the singularities of the vector field associated with the Riemann ξ {\rm{\xi }} -function.

Canonical lifts in multisymplectic De Donder–Weyl Hamiltonian field theories

We define canonical lifts of vector fields to the multisymplectic multimomentum bundles of De Donder–Weyl Hamiltonian (first-order) field theories and to the appropriate premultisymplectic embedded constraint submanifolds on which singular field theories are studied. These new canonical lifts are used to study the so-called natural Noether symmetries present in both regular and singular Hamiltonian field theories along with their associated conserved quantities obtained from Noether’s theorem. The Klein–Gordon field, the Polyakov bosonic string, and Einstein–Cartan gravity in 3+1 dimensions are analyzed in depth as applications of these concepts; as a peripheral result obtained in the analysis of the bosonic string, we provide a new geometrical interpretation of the well-known Virasoro constraint.

Remarks on geometric engineering, symmetry TFTs and anomalies

Geometric engineering is a collection of tools developed to establish dictionaries between local singularities in string theory and (supersymmetric) quantum fields. Extended operators and defects, as well as their higher quantum numbers captured by topological symmetries, can be encoded within geometric engineering dictionaries. In this paper we revisit and clarify aspects of these techniques, with special emphasis on ’t Hooft anomalies, interpreted from the SymTFT perspective as obstructions to the existence of Neumann boundary conditions. These obstructions to gauging higher symmetries are captured via higher link correlators for the SymTFT on spheres. In this work, we give the geometric engineering counterpart of this construction in terms of higher links of topological membranes. We provide a consistency check in the context of 5D SCFTs with anomalous 1-form symmetries, where we give two independent derivations of the anomaly in terms of higher links, one purely field theoretical and the other purely geometrical. Along the way, we also recover the construction of non-invertible duality defects in 4D N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 SYM from a geometric engineering perspective.

Extended legality of curved boundary integral method.

The angular spectrum method is an efficient approach for synthesizing electromagnetic beams from planar electric field distributions. The electric field definition is restricted to a plane, which can introduce inaccuracy when applying the synthesized beam to curved surface features. The angular spectrum method can also be interpreted as a pure source method defining the field symmetrically with respect to the creation plane. Recently, we generalized that symmetric field method to arbitrary source distributions, which are valid at any point on compact, regular surface Ω in R 3. We call this approach the Curved Boundary Integral method. The electromagnetic fields synthesized with this method satisfy the Helmholtz equation and are adjusted via amplitude and phase at the desired surface. The fields are obtained as a relatively simple integral. However, restrictions on where in space the synthesized field is valid were included in the mathematical proof length to avoid obscuring the main points. These restrictions can be significant depending on the shape and degree of curvature of surface Ω. In this article, we remove these restrictions so that the integral representation of the electromagnetic beam becomes valid at all points r∈R 3∖Ω, with a minor restriction. Its modification can work even on Ω. We demonstrate the importance of this extended legality with a source field parametrized into the torus surface. The electromagnetic radiation of this structure would not be valid at any point in space without this extension. Finally, we show that by changing the order of integration, the field singularity at each source point is eliminated.

Analysis of Stokes singularities using a lateral shear interferometer

Polarization and Poincaré singularities in the optical fields can be studied by analyzing the phase singularities of mathematically constructed Stokes vector fields. The wider applicability of the Stokes construction is found by exploring the generation and detection methods for various types of Stokes singularities and their analysis. Here, we detect and analyze all forms of the Stokes singularities through lateral shear interferometry. Specifically, the projections of a Stokes singularity on three pairs of orthogonal polarization basis states, defined by the eigen polarization states of Pauli’s matrices, are analyzed through unique fork patterns in the shearogram pairs. These interference patterns also provide the topological indices of the singularities. Such a self-referencing interferometric method also helps to remove the degeneracy in the Stokes index and polarization. Through both, simulations and experiments, we have analyzed specific beams represented by higher order Poincaré sphere and hybrid order Poincaré sphere topological constructs.

Singularities of Diffuse Wave Fields in Scattering Media with Refractive Index Gradients

The subject of this paper is the radiative transfer in media with gradients of the refractive index. Asymptotic expressions of the intensity distributions in the vicinity of singular directions are derived. Critical conditions for the appearance of singular intensity distributions in fields of scattered radiation in a medium are formulated. It is shown that singular radiation fields in such media can be generated, among other things, by non-singular configurations of radiation sources. The obtained results are verified by direct comparison with the results of numerical calculations using the Monte Carlo statistical modeling method.

Dissipation during crack growth in a viscoelastic material from a cohesive model for a finite specimen

AbstractIn the present paper, we extend results recently given by Ciavarella et al. (J Mech Phys Solids 169:105096, 2022) to show some actual calculations of the viscoelastic dissipation in a crack propagation at constant speed in a finite size specimen. It is usually believed that the cohesive models introduced by Knauss and Schapery and the dissipation-based theories introduced by de Gennes and Persson-Brener give very similar results for steady state crack propagation in viscoelastic materials, where usually only the asymptotic singular field is used for the stress. We show however that dissipation and the energy balance never reach a steady state, despite the constant propagation crack rate and stress intensity factor. Our loading protocol permits a rigorous solution, and implies a short phase with constant specimen elongation rate, but then possibly a very long phase of constant or decreasing elongation, which differs from typical experiments. For the external work we are therefore unable to use the de Gennes and Persson-Brener theories which suggested that the increase of effective fracture energy would go up to the ratio of instantaneous to relaxed modulus, at very fast rates. We show viscoelastic dissipation is in general a transient quantity, which can vary by orders of magnitude while the stress intensity factor is kept constant, and is largely affected by dissipation in the bulk rather than at the crack tip. The total work to break a specimen apart is found also to be possibly arbitrarily large for quite a large range of intermediate crack growth rates.

Initiation and arrest of cracks from corners in multi-chip semiconductor devices

A contemporary semiconductor device often contains multiple chips. Corners of the chips concentrate stress, and are principal sites to initiate failure. Here we propose to characterize the corners using a double cantilever beam, in which two silicon beams sandwich a row of chips. As the two beams are pulled open, a crack initiates at the corner of a chip, and runs unstably on the interface between the chip and a beam. The crack may or may not arrest, depending on various experimental conditions. We calculate energy release rate as a function of crack length by using a combination of finite element method and an analytical solution of the singular field around a corner. At a fixed applied displacement, the energy release rate is low for a short crack, peaks for a crack of intermediate length, and drops for a long crack. This non-monotonic behavior explains how a crack initiates, grows unstably, and possibly arrests. If the crack does arrest, as the two beams open further, the crack grows stably. We relate the initiation and arrest of the crack to machine compliance, specimen geometry, and flaw size. The force at which the crack initiates can be used to characterize the manufacturing process, whereas the stable growth of the crack can be used to measure interfacial toughness. It is hoped that this work will aid the development of multi-chip semiconductor devices.

On Darboux theorems for geometric structures induced by closed forms

This work reviews the classical Darboux theorem for symplectic, presymplectic, and cosymplectic manifolds (which are used to describe mechanical systems), as well as certain cases of multisymplectic manifolds, while extends the Darboux theorem in new ways to k-symplectic and k-cosymplectic manifolds (all these structures appear in the geometric formulation of first-order classical field theories). Moreover, we discuss the existence of Darboux theorems for classes of precosymplectic, k-presymplectic, k-precosymplectic, and premultisymplectic manifolds, which are the geometrical structures underlying some kinds of singular field theories, i.e. with locally non-invertible Legendre maps. Approaches to Darboux theorems based on flat connections associated with geometric structures are given, while new results on polarisations for (k-)(pre)(co)symplectic structures arise.

Tailoring focal plane component intensities of polarization singular fields in a tight focusing system

The scientific community studies tight focusing of radially and azimuthally-polarized vector beams as it is a versatile solution for many applications. We offer a new method to produce tight focusing that ensures a more uniform intensity profile in multiple dimensions, providing a more versatile and stable solution. We manipulate the polarization of the radially and azimuthally polarized vector beams to find an optimal operating point. We examine in detail optical fields whose polarization states lie on the equator of the relevant Poincaré spheres namely, the fundamental Poincaré sphere, the hybrid order Poincaré sphere (HyOPS), and the higher order Poincaré sphere. We find via simulation that the fields falling on these equators have focal plane intensity distributions characterized by a single rotation parameter α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} determining the individual state of polarization. The strengths of the component field distributions vary with α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} and can be tuned to achieve equal strengths of longitudinal (z) and transverse (x and y) components at the focal plane. Without control of this parameter (e.g., using α=0\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha =0$$\\end{document} in radially and α=π\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha =\\pi$$\\end{document} in azimuthally-polarized vector beams) intensity in x and y components are at 20% of the z component. In our solution with α=π/2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha =\\pi /2$$\\end{document}, all components are at 80% of the maximum possible intensity of z. In examining the impact of α\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\alpha$$\\end{document} on a tightly focused beam, we also found that a helicity inversion of HyOPS beams causes a rotation of 180 degree in the axial intensity distribution.

Singularities of 3D vector fields preserving the Martinet form

Изучается локальная структура векторных полей, заданных в пространстве $\mathbb{R}^3$, которые сохраняют $1$-форму Мартине $\alpha=(1+x)dy\pm z dz$. Произведена классификация их особенностей с точностью до диффеоморфизмов, сохраняющих форму $\alpha$, а также их трансверсальных разверток. В результате появляется возможность создать довольно полный список бифуркаций, которым подвергаются такие векторные поля.

On inherent hyperelastic crease

We present analyses of the inherent hyperelastic crease that entails an unstable degenerative singular perturbation field over a uniformly compressed hyperelastic half-space. The analyses reveal asymptotic far-field characteristics of the singular field that bring out admissible incremental elastic deformation mechanisms of the inherent creasing instability typically observed at a compressive strain of ∼0.354 in neo-Hookean solids. The admissible singular perturbation field has an asymptotic leading-order incremental-stress singularity decaying 1/r2 as the distance from the crease tip r→∞. In general, asymptotic incremental-deformation fields of 1/r2 stress singularity can be introduced by surface flaws. We found that the singular field creates a far-field eigenmode of three energy-release angular sectors separated by two energy-elevating sectors of incremental deformation. The far-field eigenmode braces the near-tip energy-release field of the surface flaw against the transition to an unstable self-similar expansion field of creasing. Two asymptotic far-field parameters can characterize the braced-incremental-deformation (bid) and crease fields. One is the Burgers vector of a projected subsurface dislocation of the asymptotically singular far field. The other is a dimensionless shape factor representing the ratio of the subsurface depth of the dislocation to the Burgers vector. Our analyses show that the shape factor gauges the transition from the stable bid field to the unstable inherent crease field at the crease limit point as the compression increases. Two tangential-manifold conservation integrals are revealed to identify the two parameters. At the crease-limit point, matched asymptotes lead to the ratio between the two separate scaling parameters of the asymptotic solutions of the crease-tip folding field and the leading-order far field. To this end, we introduced a novel finite element method (FEM) for simulating the crease field nearly in a hyperelastic half-space with a finite-size simulation domain. Furthermore, we uncovered with the Gent model (1996) that strain-stiffening in hyperelasticity alters the dependence of the shape factor on the compressive strain, raising crease resistance. The new findings in hyperelastic crease mechanisms will help study ruga mechanics of self-organization and design soft-material structures and skin conditions for high crease resistance.

Maxwell-Schrödinger equations in singular electromagnetic field

Maxwell-Schrödinger equations in singular electromagnetic field

Electrostatic body forces in cracked dielectrics and their implication on Maxwell stress tensors

In solid mechanics, Maxwell stresses are known to be induced if a body is exposed to magnetic and, in the case of dielectrics, electric fields. Acting as tractions at outer or inner surfaces as well as volume forces, they are superimposed with tractions and stresses due to mechanical loads and provide a more or less significant contribution, depending on loading, material properties and geometric aspects. The Maxwell stress tensor, constituting the physical and mathematical basis, however, is controversially discussed to date. Several formulations are known, most of them having been suggested more than 100 years ago. Being equivalent in vacuum, they differ qualitatively just as quantitatively in solid or fluidic matter. In particular, the dissimilar effect of body forces, emanating from a choice of established Maxwell stress tensor approaches, on crack tip loading in dielectric solids is investigated theoretically in this paper. Due to the singularity of fields involved, their impact is basically non-negligible compared to external mechanical loading. The findings obtained indicate that fracture mechanics could be the basis of an experimental validation of Maxwell stress tensors.

Nonlinear Vibration of Cracked Porous FG-GPL RC Cylindrical Panels Using a Phase-Field Crack Model

This study is concerned with the nonlinear free vibration of a cracked functionally graded porous cylindrical panel reinforced with graphene platelets by introducing a phase-field crack model. Conventional crack modeling by separating the grid nodes lying on the crack line is not only painstaking but also suffers from numerical instability. To overcome this problem, the internal crack is modeled by adopting the phase-field formulation and a virtual geometry rotation. The nonlinear numerical method is developed based on the first-order shear deformation theory incorporated with the von Kármán geometry nonlinearity in the framework of the 2-D extended natural element method, a recently introduced mesh-free method. The crack-induced singular field is represented by adopting the crack-tip singular functions, and the troublesome numerical locking is restrained by combining the MITC3+ shell concept and the shear stabilization factor. The curved shell surface is mapped to a 2-D rectangular NEM grid to avoid difficulty in defining the interpolation functions. The developed numerical method is verified through a comparison with the reference solutions, and the large-amplitude free vibration of porous cracked functionally graded grapheme platelet-reinforced cylindrical panels is profoundly examined by changing the major parameters.