Critical Exponent Proble In R2 With Neumann Boundary Conition

In this paper we prove the exostence of a solution for the following Neumann problem where is a bounded domain in R$sup:2$esup: with smooth boundary a bounded measurable function on a non–negative real number, f and g are functions of critical growth on and respectively and v is the outward unit normal to     

Finite analytic numerical method for three‐dimensional quasi‐laplace equation with conductivity in tensor form

The finite analytic numerical method for 3D quasi‐Laplace equation with conductivity in full tensor form is constructed in this article. For cubic grid system, the gradient of the potential variable will diverge when tending to the common edge joining the four grids with different conductivities. However, the potential gradient along the tangential direction is of limited value. As a consequence, the 3D quasi‐Laplace equations will behave as a quasi‐2D one. An approximate analytical solution of the 3D quasi‐Laplace equation can be found around the common edge, which is expressed as a combination of a power‐law function and a linear function. With the help of this approximate analytical solution, a 3D finite analytical numerical scheme is then constructed. Numerical examples show that the proposed numerical scheme can provide rather accurate solutions only with or subdivisions. More important, the convergent speed of the numerical scheme is independent of the conductivity heterogeneity. In contrast, when using the traditional numerical schemes, typically such as the MPFA method, the refinement ratio for the grid cell needs to increase dramatically to get an accurate result for the strong heterogeneous case.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1475–1492, 2017     

On computation of test dipoles for factorization method

In electrical impedance tomography, one tries to recover the spatial conductivity distribution inside a body from boundary measurements of current and voltage. In many important situations, the examined object has known background conductivity but is contaminated by inhomogeneities. The factorization method of Kirsch provides a tool for locating such inclusions. The computational attractiveness of the factorization technique relies heavily on efficient computation of Dirichlet boundary values of potentials created by dipole sources located inside the examined object and corresponding to the homogeneous Neumann boundary condition and to the known background conductivity. In certain simple situations, these test potentials can be written down explicitly or given with the help of suitable analytic maps, but, in general, they must be computed numerically. This work introduces an inexpensive algorithm for approximating the test potentials in the framework of real-life electrode measurements and analyzes how well this technique can be imbedded in the factorization method. The performance of the resulting fast reconstruction algorithm is tested in two spatial dimensions. The work of the second author was supported by the Academy of Finland (project 115013), the Finnish Funding Agency for Technology and Innovation (project 40084/06), the Finnish Cultural Foundation and the Finnish Foundation for Technology Promotion.     

An improved continuum skin and proximity effect model for hexagonally packed wires

This work presents approximate but closed-form expressions for “effective” complex-valued magnetic permeability and electric conductivity that represent the effects of proximity and skin effect losses in wound coil with hexagonally packed wires. Previous work is extended by providing improved accuracy versus finite element results for effective permeability and by providing an expression for effective conductivity, which was previously neglected. These material properties can then be used in 2D/axisymmetric finite element models in which the coil is modeled as a coarsely meshed, homogeneous region (i.e., removing the need for modeling each turn in the coil).     

拟共形映射和弱John域

作为John域的推广,本文定义了弱John域,并讨论了弱John域与拟圆、弱John域与拟共形 映射之间的关系,得到(1)若(?)。中的Jordan域D和它的外部 均是弱John域,则D 是拟圆;(2)R2中的弱John域是拟共不变的;(3)R2中的有界拟圆必是弱John域．最后构造例子 说明R2中的无界拟圆不一定是弱John域．     

Schrödinger Operators on Zigzag Nanotubes

We consider the Schr?dinger operator with a periodic potential on quasi-1D models of zigzag single-wall carbon nanotubes. The spectrum of this operator consists of an absolutely continuous part (intervals separated by gaps) plus an infinite number of eigenvalues with infinite multiplicity. We describe all compactly supported eigenfunctions with the same eigenvalue. We define a Lyapunov function, which is analytic on some Riemann surface. On each sheet, the Lyapunov function has the same properties as in the scalar case, but it has branch points, which we call resonances. We prove that all resonances are real. We determine the asymptotics of the periodic and antiperiodic spectrum and of the resonances at high energy. We show that there exist two types of gaps: i) stable gaps, where the endpoints are periodic and anti-periodic eigenvalues, ii) unstable (resonance) gaps, where the endpoints are resonances (i.e., real branch points of the Lyapunov function). We describe all finite gap potentials. We show that the mapping: potential all eigenvalues is a real analytic isomorphism for some class of potentials. Submitted: October 5, 2006. Accepted: December 15, 2006.     

The String Method as a Dynamical System

This paper is devoted to the theoretical analysis of the zero-temperature string method, a scheme for identifying minimum energy paths (MEPs) on a given energy landscape. By definition, MEPs are curves connecting critical points on the energy landscape which are everywhere tangent to the gradient of the potential except possibly at critical points. In practice, MEPs are mountain pass curves that play a special role, e.g., in the context of rare reactive events that occur when one considers a steepest descent dynamics on the potential perturbed by a small random noise. The string method aims to identify MEPs by moving each point of the curve by steepest descent on the energy landscape. Here we address the question of whether such a curve evolution necessarily converges to an MEP. Surprisingly, the answer is no, for an interesting reason: MEPs may not be isolated, in the sense that there may be families of them that can be continuously deformed into one another. This degeneracy is related to the presence of critical points of Morse index 2 or higher along the MEP. In this paper, we elucidate this issue and completely characterize the limit set of a curve evolving by the string method. We establish rigorously that the limit set of such a curve is again a curve when the MEPs are isolated. We also show under the same hypothesis that the string evolution converges to an MEP. However, we identify and classify situations where the limit set is not a curve and may contain higher dimensional parts. We present a collection of examples where the limit set of a path contains a 2D region, a 2D surface, or a region of an arbitrary dimension up to the dimension of the space. In some of our examples the evolving path wanders around without converging to its limit set. In other examples it fills a region, converging to its limit set, which is not an MEP.     

On the Su–Schrieffer–Heeger model of electron transport: Low-temperature optical conductivity by the Mellin transform

We describe the low-temperature optical conductivity as a function of frequency for a quantum-mechanical system of electrons that hop along a polymer chain. To this end, we invoke the Su–Schrieffer–Heeger tight-binding Hamiltonian for noninteracting spinless electrons on a one-dimensional (1D) lattice. Our goal is to show via asymptotics how the interband conductivity of this system behaves as the smallest energy bandgap tends to close. Our analytical approach includes: (i) the Kubo-type formulation for the optical conductivity with a nonzero damping due to microscopic collisions, (ii) reduction of this formulation to a 1D momentum integral over the Brillouin zone, and (iii) evaluation of this integral in terms of elementary functions via the three-dimensional Mellin transform with respect to key physical parameters and subsequent inversion in a region of the respective complex space. Our approach reveals an intimate connection of the behavior of the conductivity to particular singularities of its Mellin transform. The analytical results are found in good agreement with direct numerical computations.     

A Fourier method for computing the perturbation to a two‐dimensional electric potential by a conductor of arbitrary shape

Introducing an electric conductor into a region pervaded by an initial electric potential perturbs that potential by inducing a charge distribution on the conductor's surface, necessary to guarantee that the surface is an equipotential of the total potential. Some numerical method is required to compute the perturbation potential, when the conductor's shape does not admit a standard analytic solution. For two‐dimensional situations, a method is proposed for solving for the perturbation potential that involves expansion of the boundary perturbation potential and its normal derivative as truncated Fourier series. This boundary potential is known to within an additive constant from the requirement that its sum with the initial potential must be a constant. The standard representation theorem for the Dirichlet problem gives a consistency relation between the boundary function and its normal derivative, which here becomes a set of linear algebraic relations between Fourier series coefficients, with matrix entries found by appropriate applications of the fast Fourier transform. These are solved for the boundary derivative coefficients; at any point exterior to the conductor, the perturbation potential can then be evaluated from the two sets of Fourier coefficients, using further application of the fast Fourier transform. Examples are shown for two conductor shapes, with several initial potentials. © 2001 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 17: 673–683, 2001     

A new linear ordering of fuzzy numbers on subsets of {{\mathcal F}({\pmb{\mathbb{R}}}})

We introduce a novel linear order on every family of fuzzy numbers which satisfies the assumption that their modal values must be all different and must form a compact subset of . A distinct new feature is that our linear determined procedure employs the corresponding order of a class interval associated with a confidence measure which seems intuitively anticipated. It is worthy noting that although we start from an entirely different rationale, we introduce a fuzzy ordering which initially coincides with the one established earlier by Ramik and Rimanek. However, this fuzzy ordering does not apply when the supports of the fuzzy numbers overlap. In order to cover such cases we extent this initial fuzzy ordering to the “extended fuzzy order” (XFO). This new XFO method includes a possibility and a necessity measure which are compared with the widely accepted PD and NSD indices of D. Dubois and H. Prade. The comparison shows that our possibility and necessity measures comply better with our intuition.     

The approximate and the Clarke subdifferentials can be different everywhere

We prove that, for a Lipschitz function on , n2, the approximate and the Clarke subdifferentials can differ everywhere. This completely answers a question by A.D. Ioffe, which was partially answered by G. Katriel.     

On strong convergence to 3D steady vortex sheets

In this paper, we first establish a strong convergence criterion of approximate solutions for the 3D steady incompressible Euler equations. For axisymmetric flows, under the assumption that the vorticity is of one sign and uniformly bounded in L¹ space, we obtain a sufficient and necessary condition for the strong convergence in of approximate solutions. Furthermore, for one-sign and L¹-bounded vorticity, it is shown that if a sequence of approximate solutions concentrates at an isolated point in (r,z)-plane, then the concentration point can appear neither in the region near the axis (including the symmetry axis itself) nor in the region far away from the axis. Finally, we present an example of approximates solutions which converge strongly in by using Hill's spherical vortex.     

The lift on a small sphere in a linear shear flow near the interface of two immiscible fluids

Analytical expressions have been derived which predict, to lowest order, the inertial lift and the lateral migration velocity of a rigid sphere translating and rotating in a linear shear flow field near the flat interface of two immiscible fluids. This asymptotic analysis is primarily based on the assumption that the two Reynolds numbers defined by the gap width between the interface and the sphere, the shear rate and the translational slip velocity with which the spherical particle moves parallel to the interface are small. Furthermore, the radius of the sphere is assumed to be small compared to the gap width. To leading order in this creeping flow regime, the linear Stokes equations are obtained and a symmetry argument can be used to show that the Stokes solution does not predict any lift force. The transverse force experienced by the sphere and its migration velocity are due to the small but finite inertial terms in the Navier-Stokes equations, which can be studied by perturbation techniques. By applying a Green's function approach and matched asymptotic methods, which also incorporate the effects of the outer Oseen-like flow regime, the three components comprising the lift velocity have been calculated in closed form: the one induced by the shear rate only, the purely slip induced one and the one due to the interaction of the slip velocity with the shear flow field. The thus obtained expressions for the case of two immiscible fluids with arbitrary density and viscosity ratios extend the results that already exist in the literature for other flow configurations, such as an unbounded shear flow field [1] or a wall-bounded one, where the wall lies either within the leading order Stokes region [2] or in the outer Oseen region [3]. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An asymptotic formula for the voltage potential in the case of a near-surface conductivity inclusion

We rigorously derive an asymptotic expansion of the steady-state voltage potentials in the presence of a conductivity inclusion of small volume that is close to a planar surface. This new formula is motivated by the practically important inverse problem of imaging a conductivity inclusion near a planar interface. Partly supported by ACI Jeunes Chercheurs (0693) from the Ministry of Education and Scientific Research, France. Partially supported by grant R02-2003-000-10012-0 from the Korea Science and Engineering Foundation. Received: December 2, 2003; revised: January 18, 2004     

Determining conductivity by boundary measurements

In [2], A. P. Caldéron posed the following question: can one determine the heat conductivity of an object from static temperature and heat flux measurements at the boundary? We show that such measurements uniquely determine the conductivity and all of its derivatives at the boundary.     

Schrödinger potentials solvable in terms of the confluent Heun functions

We show that if the potential is proportional to an energy-independent continuous parameter, then there exist 15 choices for the coordinate transformation that provide energy-independent potentials whose shape is independent of that parameter and for which the one-dimensional stationary Schrödinger equation is solvable in terms of the confluent Heun functions. All these potentials are also energy-independent and are determined by seven parameters. Because the confluent Heun equation is symmetric under transposition of its regular singularities, only nine of these potentials are independent. Five of the independent potentials are different generalizations of either a hypergeometric or a confluent hypergeometric classical potential, one potential as special cases includes potentials of two hypergeometric types (the Morse confluent hypergeometric and the Eckart hypergeometric potentials), and the remaining three potentials include five-parameter conditionally integrable confluent hypergeometric potentials. Not one of the confluent Heun potentials, generally speaking, can be transformed into any other by a parameter choice.     

Solving dynamical equations in general homogeneous isotropic cosmologies with a scalaron

We consider gauge-dependent dynamical equations describing homogeneous isotropic cosmologies coupled to a scalar field ψ (scalaron). For flat cosmologies (k = 0), we analyze the gauge-independent equation describing the differential χ(α) ≡ ψ (a) of the map of the metric a to the scalaron field ψ, which is the main mathematical characteristic of a cosmology and locally defines its portrait in the so-called a version. In the more customary ψ version, the similar equation for the differential of the inverse map \(\bar \chi (\psi ) \equiv \chi ^{ - 1} (\alpha )\) is solved in an asymptotic approximation for arbitrary potentials v(ψ). In the flat case, \(\bar \chi (\psi )\) and χ^?1(α) satisfy first-order differential equations depending only on the logarithmic derivative of the potential, v(ψ)/v(ψ). If an analytic solution for one of the χ functions is known, then we can find all characteristics of the cosmological model. In the α version, the full dynamical system is explicitly integrable for k ≠ 0 with any potential v(α) ≡ v[ψ(α)] replacing v(ψ). Until one of the maps, which themselves depend on the potentials, is calculated, no sort of analytic relation between these potentials can be found. Nevertheless, such relations can be found in asymptotic regions or by perturbation theory. If instead of a potential we specify a cosmological portrait, then we can reconstruct the corresponding potential. The main subject here is the mathematical structure of isotropic cosmologies. We also briefly present basic applications to a more rigorous treatment of inflation models in the framework of the α version of the isotropic scalaron cosmology. In particular, we construct an inflationary perturbation expansion for χ. If the conditions for inflation to arise are satisfied, i.e., if v > 0, k = 0, χ² < 6, and χ(α) satisfies a certain boundary condition as α→-∞, then the expansion is invariant under scaling the potential, and its first terms give the standard inflationary parameters with higher-order corrections.     

Stationary distribution of a two-dimensional SRBM: geometric views and boundary measures

We present three sets of results for the stationary distribution of a two-dimensional semimartingale-reflecting Brownian motion (SRBM) that lives in the non-negative quadrant. The SRBM data can equivalently be specified by three geometric objects, an ellipse and two lines, in the two-dimensional Euclidean space. First, we revisit the variational problem (VP) associated with the SRBM. Building on Avram et al. (Queueing Syst. 37: 259–289, 2001), we show that the value of the VP at a point in the quadrant is equal to the optimal value of a linear function over a convex domain. Depending on the location of the point, the convex domain is either $\mathcal{D}^{(1)}$ or $\mathcal{D}^{(2)}$ or $\mathcal{D}^{(1)}\cap \mathcal{D}^{(2)},$ where each $\mathcal{D}^{(i)},$ $i=1, 2,$ can easily be described by the three geometric objects. Our results provide a geometric interpretation for the value function of the VP and allow one to see geometrically when one edge of the quadrant has influence on the optimal path traveling from the origin to a destination point. Second, we provide a geometric condition that characterizes the existence of a product form stationary distribution. Third, we establish exact tail asymptotics of two boundary measures that are associated with the stationary distribution; a key step in our proof is to sharpen two asymptotic inversion lemmas in Dai and Miyazawa (Stoch. Syst. 1:146–208, 2011) which allow one to infer the exact tail asymptotic of a boundary measure from the singularity of its moment-generating function.     

Rigidity of planar tilings

Summary Aperturbation of a tiling of a region inR ⁿ is a set of isometries, one applied to each tile, so that the images of the tiles tile the same region.We show that a locally finite tiling of an open region inR ² with tiles which are closures of their interiors isrigid in the following sense: any sufficiently small perturbation of the tiling must have only earthquake-type discontinuities, that is, the discontinuity set consists of straight lines and arcs of circles, and the perturbation near such a curve shifts points along the direction of that curve.We give an example to show that this type of rigidity does not hold inR ⁿ , forn>2.Using rigidity in the plane we show that any tiling problem with a finite number of tile shapes (which are topological disks) is equivalent to a polygonal tiling problem, i.e. there is a set of polygonal shapes with equivalent tiling combinatorics.Oblatum 19-III-1991