skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Polynomial inclusions: Definitions, applications, and open problems
Predictive modeling in physical science and engineering is mostly based on solving certain partial differential equations where the complexity of solutions is dictated by the geometry of the domain. Motivated by the broad applications of explicit solutions for spherical and ellipsoidal domains, in particular, the Eshelby’s solution in elasticity, we propose a generalization of ellipsoidal shapes called polynomial inclusions. A polynomial inclusion (or -inclusion for brevity) of degree is defined as a smooth, connected and bounded body whose Newtonian potential is a polynomial of degree inside the body. From this viewpoint, ellipsoids are identified as the only -inclusions of degree two; many fundamental problems in various physical settings admit simple closed-form solutions for general -inclusions as for ellipsoids. Therefore, we anticipate that -inclusions will be useful for applications including predictive materials models, optimal designs, and inverse problems. However, the existence of p-inclusions beyond degree two is not obvious, not to mention their explicit algebraic parameterizations. In this work, we explore alternative definitions and properties of p-inclusions in the context of potential theory. Based on the theory of variational inequalities, we show that -inclusions do exist for certain polynomials, though a complete characterization remains open. We reformulate the determination of surfaces of -inclusions as nonlocal geometric flows which are convenient for numerical simulations and studying geometric properties of -inclusions. In two dimensions, by the method of conformal mapping we find an explicit algebraic parameterization of p-inclusions. We also propose a few open problems whose solution will deepen our understanding of relations between domain geometry, Newtonian potentials, and solutions to general partial differential equations. We conclude by presenting examples of applications of -inclusions in the context of Eshelby inclusion problems and magnet designs.  more » « less
Award ID(s):
2306254
PAR ID:
10500897
Author(s) / Creator(s):
;
Publisher / Repository:
SccienceDirect
Date Published:
Journal Name:
Journal of the Mechanics and Physics of Solids
Volume:
181
Issue:
C
ISSN:
0022-5096
Page Range / eLocation ID:
105440
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We introduce Geometric Neural Operators (GNPs) for data-driven deep learning of geometric features for tasks in non-euclidean settings. We present a formulation for accounting for geometric contributions along with practical neural network architectures and factorizations for training. We then demonstrate how GNPs can be used (i) to estimate geometric properties, such as the metric and curvatures of surfaces, (ii) to approximate solutions of geometric partial differential equations on manifolds, and (iii) to solve Bayesian inverse problems for identifying manifold shapes. These results show a few ways GNPs can be used for incorporating the roles of geometry in the data-driven learning of operators. 
    more » « less
  2. This paper presents a scalable tensor-based approach to computing controllability and observability-type energy functions for nonlinear dynamical systems with polynomial drift and linear input and output maps. Using Kronecker product polynomial expansions, we convert the Hamilton- Jacobi-Bellman partial differential equations for the energy functions into a series of algebraic equations for the coefficients of the energy functions. We derive the specific tensor structure that arises from the Kronecker product representation and analyze the computational complexity to efficiently solve these equations. The convergence and scalability of the proposed energy function computation approach is demonstrated on a nonlinear reaction-diffusion model with cubic drift nonlinearity, for which we compute degree 3 energy function approximations in n = 1023 dimensions and degree 4 energy function approximations in n = 127 dimensions. 
    more » « less
  3. Abstract We study the structure of the solution sets in universal differential fields of certain differential equations of order two, the Poizat equations, which are particular cases of Liénard equations. We give a necessary and sufficient condition for strong minimality for equations in this class and a complete classification of the algebraic relations for solutions of strongly minimal Poizat equations. We also give an analysis of the non-strongly minimal cases as well as applications concerning the Liouvillian and Pfaffian solutions of some Liénard equations. 
    more » « less
  4. Lump solutions are analytical rational function solutions localized in all directions in space. We analyze a class of lump solutions, generated from quadratic functions, to nonlinear partial differential equations. The basis of success is the Hirota bilinear formulation and the primary object is the class of positive multivariate quadratic functions. A complete determination of quadratic functions positive in space and time is given, and positive quadratic functions are characterized as sums of squares of linear functions. Necessary and sufficient conditions for positive quadratic functions to solve Hirota bilinear equations are presented, and such polynomial solutions yield lump solutions to nonlinear partial differential equations under the dependent variable logarithmic derivative transformations. Applications are made for a few generalized KP and BKP equations. 
    more » « less
  5. It is well known that the Painlevé equations can formally degenerate to autonomous differential equations with elliptic function solutions in suitable scaling limits. A way to make this degeneration rigorous is to apply Deift-Zhou steepest-descent techniques to a Riemann-Hilbert representation of a family of solutions. This method leads to an explicit approximation formula in terms of theta functions and related algebro-geometric ingredients that is difficult to directly link to the expected limiting differential equation. However, the approximation arises from an outer parametrix that satisfies relatively simple conditions. By applying a method that we learned from Alexander Its, it is possible to use these simple conditions to directly obtain the limiting differential equation, bypassing the details of the algebro-geometric solution of the outer parametrix problem. In this paper, we illustrate the use of this method to relate an approximation of the algebraic solutions of the Painlevé-III (D$$_7$$) equation valid in the part of the complex plane where the poles and zeros of the solutions asymptotically reside to a form of the Weierstraß equation. 
    more » « less