On an integral equation approach for the exterior Robin problem for the Helmholtz equation

A boundary integral equation for the exterior Robin problem for Helmholtz's equation is analyzed in this paper. This integral operator is not compact. A proof based on a suitable regularization of this integral operator and the Fredholm alternative for the regularized compact operator was given by other authors. In this paper, we will give a direct existence and uniqueness proof for the boundary non-compact integral equation in the space settings C^1,λ(S) and C^0,λ(S), where S is a closed bounded smooth surface.     

Integral Equation Method for the First and Second Problems of the Stokes System

Using the integral equation method we study solutions of boundary value problems for the Stokes system in Sobolev space H ¹(G) in a bounded Lipschitz domain G with connected boundary. A solution of the second problem with the boundary condition $\partial {\bf u}/\partial {\bf n} -p{\bf n}={\bf g}$ is studied both by the indirect and the direct boundary integral equation method. It is shown that we can obtain a solution of the corresponding integral equation using the successive approximation method. Nevertheless, the integral equation is not uniquely solvable. To overcome this problem we modify this integral equation. We obtain a uniquely solvable integral equation on the boundary of the domain. If the second problem for the Stokes system is solvable then the solution of the modified integral equation is a solution of the original integral equation. Moreover, the modified integral equation has a form f?+?S f?=?g, where S is a contractive operator. So, the modified integral equation can be solved by the successive approximation. Then we study the first problem for the Stokes system by the direct integral equation method. We obtain an integral equation with an unknown ${\bf g}=\partial {\bf u}/\partial {\bf n} -p{\bf n}$ . But this integral equation is not uniquely solvable. We construct another uniquely solvable integral equation such that the solution of the new eqution is a solution of the original integral equation provided the first problem has a solution. Moreover, the new integral equation has a form ${\bf g}+\tilde S{\bf g}={\bf f}$ , where $\tilde S$ is a contractive operator, and we can solve it by the successive approximation.     

Bounds for the number of points on curves over finite fields

Let U be a bounded open subset of ?^d, d ≥ 2 and f ∈ C(?U). The Dirichlet solution f^CU of the Dirichlet problem associated with the Laplace equation with a boundary condition f is not continuous on the closure ū of U in general if U is not regular but it is always Baire-one.Let H(U) be the space of all functions continuous on the closure ū and harmonic on U and F(H(U)) be the space of uniformly bounded absolutely convergent series of functions in H(U). We prove that f^CU can be obtained as a uniform limit of a sequence of functions in F(H(U)). Thus f^CU belongs to the subclass B_1/2 of Baire-one functions studied for example in [8]. This is not only an improvement of the result obtained in [10] but it also shows that the Dirichlet solution on the closure ū can share better properties than to be only a Baire-one function. Moreover, our proof is more elementary than that in [10].A generalization to the abstract context of simplicial function space on a metrizable compact space is provided.We conclude the paper with a brief discussion on the solvability of the abstract Dirichlet problem with a boundary condition belonging to the space of differences of bounded semicontinuous functions complementing the results obtained in [17].     

Additive operators preserving rank-additivity on symmetry matrix spaces

We characterize the additive operators preserving rank-additivity on symmetry matrix spaces. LetS _n(F) be the space of alln×n symmetry matrices over a fieldF with 2,3 ∈F ^*, thenT is an additive injective operator preserving rank-additivity onS _n(F) if and only if there exists an invertible matrixU∈M _n(F) and an injective field homomorphism ? ofF to itself such thatT(X)=cUX ^?U^T, ?X=(x_ij)∈S_n(F) wherec∈F ^*,X ^?=(?(x _ij)). As applications, we determine the additive operators preserving minus-order onS _n(F) over the fieldF.     

Continuous time random walks and the Cauchy problem for the heat equation

We deal with anomalous diffusions induced by continuous time random walks - CTRW in ?ⁿ. A particle moves in ?ⁿ in such a way that the probability density function u(·, t) of finding it in region Ω of ?ⁿ is given by ∫_Ωu(x, t)dx. The dynamics of the diffusion is provided by a space time probability density J(x, t) compactly supported in {t ≥ 0}. For t large enough, u satisfies the equation

$$u\left( {x,t} \right) = \left[ {\left( {J - \delta } \right)*u} \right]\left( {x,t} \right)$$

, where δ is the Dirac delta in space-time. We give a sense to a Cauchy type problem for a given initial density distribution f. We use Banach fixed point method to solve it and prove that under parabolic rescaling of J, the equation tends weakly to the heat equation and that for particular kernels J, the solutions tend to the corresponding temperatures when the scaling parameter approaches 0.

     

Regularization of a three-element functional equation

In this paper we study the three-element functional equation

$(V\Phi )(z) \equiv \Phi (iz) + \Phi ( - iz) + G(z)\Phi \left( {\frac{1}{z}} \right) = g(z), z \in R,$

, subject to

$R: = \{ z:\left| z \right| < 1, \left| {\arg z} \right| < \frac{\pi }{4}\} .$

We assume that the coefficients G(z) and g(z) are holomorphic in R and their boundary values G ⁺(t) and g ⁺(t) belong to H(Γ), G(t)G(t ^?1) = 1. We seek for solutions Φ(z) in the class of functions holomorphic outside of \(\bar R\) such that they vanish at infinity and their boundary values Φ^?(t) also belong to H(Γ). Using the method of equivalent regularization, we reduce the problem to the 2nd kind integral Fredholm equation.

     

Fundamental solution for linear two-point boundary value problem

The quadratic functional minimization with differential restrictions represented by the command linear systems is considered. The optimal solution determination implies the solving of a linear problem with two points boundary values. The proposed method implies the construction of a fundamental solution S(t)—a n×n matrix—and of a vector h(t) defining an adjoint variable λ(t) depending of the state variable x(t). From the extremum necessary conditions it is obtained the Riccati matrix differential equation having the S(t) as unknown fundamental solution is obtained. The paper analyzes the existence of the Riccati equation solution S(t) and establishes as the optimal solution of the proposed optimum problem. Also a superior limit of the minimum for the considered quadratic functionals class are evaluated.     

Some properties of the operators of potential theory and their application to the investigation of the basic equation of electrostatics and magnetostatics

Some new properties are proved for the operatorB ^* of the direct value of the potential of a double layer on a closed surfaceS=?ω, in particular the existence inH ^1/2(S) of a basis of eigenfunctions. On the basis of these properties it is proved that the vector integral equation $$\alpha {\rm M}(x) + \nabla \int\limits_\Omega {M(y)} \nabla _y |x - y|dy = H(x), \alpha \geqslant 0,\Omega \subset R^3 ,$$ which is encountered in classical problems of electro- and magnetostatics, is equivalent to the well-known scalar equation with operatorB ^*. The properties of the operator on the left-hand side and of the solutions of the vector equation are investigated.     

Existence and nonuniqueness of solutions to a functional-differential equation

We examine the functional-differential equation Δu(x) — div(u(H(x))f (x)) = 0 on a torus which is a generalization of the stationary Fokker-Planck equation. Under sufficiently general assumptions on the vector field f and the map H, we prove the existence of a nontrivial solution. In some cases the subspace of solutions is established to be multidimensional.     

A multiscale Galerkin method for the hypersingular integral equation reduced by the harmonic equation

The aim of this paper is to investigate the numerical solution of the hypersingular integral equation reduced by the harmonic equation. First, we transform the hypersingular integral equation into 2π-periodic hypersingular integral equation with the map x=cot(θ/2). Second, we initiate the study of the multiscale Galerkin method for the 2π-periodic hypersingular integral equation. The trigonometric wavelets are used as trial functions. Consequently, the 2j+1 × 2j+1 stiffness matrix Kj can be partitioned j×j block matrices. Furthermore, these block matrices are zeros except main diagonal block matrices. These main diagonal block matrices are symmetrical and circulant matrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform and the inverse fast Fourier transform instead of the inverse matrix. Finally, we provide several numerical examples to demonstrate our method has good accuracy even though the exact solutions are multi-peak and almost singular.     

Integral representations and summations of the modified Struve function

It is known that the Struve function H _ν and the modified Struve function L _ν are closely connected to the Bessel function of the first kind J _ν and to the modified Bessel function of the first kind I _ν and possess representations through higher transcendental functions like the generalized hypergeometric ₁ F ₂ and the Meijer G function. Also, the NIST project and Wolfram formula collection contain a set of Kapteyn type series expansions for L _ν (x). In this paper firstly, we obtain various another type integral representation formulae for L _ν (x) using the technique developed by D. Jankov and the authors. Secondly, we present some summation results for different kind of Neumann, Kapteyn and Schlömilch series built by I _ν (x) and L _ν (x) which are connected by a Sonin–Gubler formula, and by the associated modified Struve differential equation. Finally, solving a Fredholm type convolutional integral equation of the first kind, Bromwich–Wagner line integral expressions are derived for the Bessel function of first kind J _ν and for an associated generalized Schlömilch series.     

Finite element method for space-time fractional diffusion equation

In this paper, we consider two types of space-time fractional diffusion equations(STFDE) on a finite domain. The equation can be obtained from the standard diffusion equation by replacing the second order space derivative by a Riemann-Liouville fractional derivative of order β (1 < β ≤ 2), and the first order time derivative by a Caputo fractional derivative of order γ (0 < γ ≤ 1). For the 0 < γ < 1 case, we present two schemes to approximate the time derivative and finite element methods for the space derivative, the optimal convergence rate can be reached O(τ^2?γ + h²) and O(τ² + h²), respectively, in which τ is the time step size and h is the space step size. And for the case γ = 1, we use the Crank-Nicolson scheme to approximate the time derivative and obtain the optimal convergence rate O(τ² + h²) as well. Some numerical examples are given and the numerical results are in good agreement with the theoretical analysis.     

Homogenization of a nonstationary model equation of electrodynamics

In L ₂(?₃;?₃), we consider a self-adjoint operator ? _ε , ε > 0, generated by the differential expression curl η(x/ε)^?1 curl??ν(x/ε) div. Here the matrix function η(x) with real entries and the real function ν(x) are periodic with respect to some lattice, are positive definite, and are bounded. We study the behavior of the operators cos(τ? _ε ^1/2 ) and ? _ε ^?1/2 sin(τ? _ε ^1/2 ) for τ ∈ ? and small ε. It is shown that these operators converge to cos(τ(?⁰)^1/2) and (?₀)^?1/2 sin(τ(?₀)^1/2), respectively, in the norm of the operators acting from the Sobolev space H ^s (with a suitable s) to ?₂. Here ?⁰ is an effective operator with constant coefficients. Error estimates are obtained and the sharpness of the result with respect to the type of operator norm is studied. The results are used for homogenizing the Cauchy problem for the model hyperbolic equation ? _τ ² v _ε = ?? _ε v _ε , div v _ε = 0, appearing in electrodynamics. We study the application to a nonstationary Maxwell system for the case in which the magnetic permeability is equal to 1 and the dielectric permittivity is given by the matrix η(x/ε).     

Laudal-type theorems in P

In this paper we classify all integral, non-degenerate, locally Cohen-Macaulay subvarieties in P^N, whose general complementary section is a complete intersection set of points: they are either complete intersections or curves on a quadric surface in P³ or degree 4 arithmetically Buchsbaum surfaces in P⁴ (i.e. the Veronese surface or a degeneration of it). As a consequence we show that every locally Cohen-Macaulay threefold in P^S of degree 4 is a complete intersection.Moreover, we obtain a generalization of Laudal's Lemma to threefolds in P⁵ and fourfolds in P⁶, which gives a bound on the degree of a codimension 2, integral subvariety X in P^N, depending both on N and a non-lifting level s of X.     

Galerkin finite element method and error analysis for the fractional cable equation

The cable equation is one of the most fundamental equations for modeling neuronal dynamics. These equations can be derived from the Nernst-Planck equation for electro-diffusion in smooth homogeneous cylinders. Fractional cable equations are introduced to model electrotonic properties of spiny neuronal dendrites. In this paper, a Galerkin finite element method(GFEM) is presented for the numerical simulation of the fractional cable equation(FCE) involving two integro-differential operators. The proposed method is based on a semi-discrete finite difference approximation in time and Galerkin finite element method in space. We prove that the numerical solution converges to the exact solution with order O(τ+h^l+1) for the lth-order finite element method. Further, a novel Galerkin finite element approximation for improving the order of convergence is also proposed. Finally, some numerical results are given to demonstrate the theoretical analysis. The results show that the numerical solution obtained by the improved Galerkin finite element approximation converges to the exact solution with order O(τ²+h^l+1).     

The existence and uniqueness of the local solution for a Camassa-Holm type equation

A shallow water equation of Camassa-Holm type, containing nonlinear dissipative effect, is investigated. Using the techniques of the pseudoparabolic regularization and some prior estimates derived from the equation itself, we establish the existence and uniqueness of its local solution in Sobolev space H^s(R) with . Meanwhile, a new lemma and a sufficient condition which guarantee the existence of solutions of the equation in lower order Sobolev space H^s with are presented.     

Classification of solutions for an integral equation

Let n be a positive integer and let 0 < α < n. Consider the integral equation We prove that every positive regular solution u(x) is radially symmetric and monotone about some point and therefore assumes the form with some constant c = c(n, α) and for some t > 0 and x₀ ? ?ⁿ. This solves an open problem posed by Lieb 12 . The technique we use is the method of moving planes in an integral form, which is quite different from those for differential equations. From the point of view of general methodology, this is another interesting part of the paper. Moreover, we show that the family of well‐known semilinear partial differential equations is equivalent to our integral equation (0.1), and we thus classify all the solutions of the PDEs. © 2005 Wiley Periodicals, Inc.