A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

P-determinant regularization method for elliptic boundary problems

An expression for thep-determinant of the quotient of two differential elliptic operators with boundary conditions is given in terms of the boundary values of their solutions. Applications to physical examples are considered.     

A spectral boundary integral method for flowing blood cells

A spectral boundary integral method for simulating large numbers of blood cells flowing in complex geometries is developed and demonstrated. The blood cells are modeled as finite-deformation elastic membranes containing a higher viscosity fluid than the surrounding plasma, but the solver itself is independent of the particular constitutive model employed for the cell membranes. The surface integrals developed for solving the viscous flow, and thereby the motion of the massless membrane, are evaluated using an O(NlogN)$O(N\log N)$ particle-mesh Ewald (PME) approach. The cell shapes, which can become highly distorted under physiologic conditions, are discretized with spherical harmonics. The resolution of these global basis functions is, of course, excellent, but more importantly they facilitate an approximate de-aliasing procedure that stabilizes the simulations without adding any numerical dissipation or further restricting the permissible numerical time step. Complex geometry no-slip boundaries are included using a constraint method that is coupled into an implicit system that is solved as part of the time advancement routine. The implementation is verified against solutions for axisymmetric flows reported in the literature, and its accuracy is demonstrated by comparison against exact solutions for relaxing surface deformations. It is also used to simulate flow of blood cells at 30% volume fraction in tubes between 4.9 and 16.9 μm in diameter. For these, it is shown to reproduce the well-known non-monotonic dependence of the effective viscosity on the tube diameter.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

Pseudo-differential projections and the topology of certain spaces of elliptic boundary value problems

We calculate the homotopy groups of the space of elliptic boundary value problems for an elliptic differential operatorA of a first order and of the space of elliptic self-adjoint boundary value problems whenA is a formally self-adjoint. In particular we show that the spectral flow of anS ¹ family of self-adjoint elliptic boundary value problems is well defined. This provides some information on spectral properties along the lines of the Vafa-Witten approach to spectral inequalities.     

Homotopy deform method for reproducing kernel space for nonlinear boundary value problems

In this paper, the combination of homotopy deform method (HDM) and simplified reproducing kernel method (SRKM) is introduced for solving the boundary value problems (BVPs) of nonlinear differential equations. The solution methodology is based on Adomian decomposition and reproducing kernel method (RKM). By the HDM, the nonlinear equations can be converted into a series of linear BVPs. After that, the simplified reproducing kernel method, which not only facilitates the reproducing kernel but also avoids the time-consuming Schmidt orthogonalization process, is proposed to solve linear equations. Some numerical test problems including ordinary differential equations and partial differential equations are analysed to illustrate the procedure and confirm the performance of the proposed method. The results faithfully reveal that our algorithm is considerably accurate and effective as expected.     

A variational principle for boundary value problems in kinetic theory

A variational principle which applies directly to the integrodifferential form of the linearized Boltzmann equation is introduced. Extremely general boundary conditions and collision terms are allowed. For a class of interesting problems, the value of the functional to be varied is shown to be closely related to quantities of great physical interest. The formalism is applied to the treatment of plane Couette flow for different forms of the collision term (BGK model, rigid spheres, Maxwell's molecules).Research sponsored by the Air Force Office of Scientific Research under contract F 61(052)-68-C-0020, through the European Office of Aerospace Research, OAR, United States Air Force.     

Extremal points for fractional boundary value problems

This article is concerned with characterizing the first extremal point, b₀, for a Riemann–Liouville fractional boundary value problem, D^α₀₊y + p(t)y = 0, 0 < t < b, y(0) = y′(0) = y″(b) = 0, 2 < α ≤ 3, by applying the theory of u₀-positive operators with respect to a suitable cone in a Banach space. The key argument is that a mapping, which maps a linear, compact operator, depending on b to its spectral radius, is continuous and strictly increasing as a function of b. Furthermore, an application to a nonlinear case is given.     

An iterative method to solve acoustic scattering problems using a boundary integral equation

In this work, a simple iterative method to solve the acoustic scattering/radiation problems using the boundary integral equation (BIE) formulation is presented. The operator equation obtained in the BIE formulation is converted into a matrix equation using the well-known method of moments solution procedure. The present method requires much fewer mathematical operations per iteration when compared to other available iterative methods. Further, the present iterative method can easily handle multiple incident fields, a highly desirable feature not available in any other iterative method, much the same way as direct solution techniques. Several numerical examples are presented to illustrate the efficiency and accuracy of the method.     

研究椭圆积分物理问题的一种方法

在一些力学问题的研究中，严格的求解必须进行椭圆积分，这使问题变得复杂．但如果运用被积函数的某些性质，通过简单的数学运算，就可以化难为易，得到问题定性的、大致的结果．     

Solution of singular two-point boundary value problems using differential transformation method

In this Letter, we implemented relatively new, exact series method of solution known as the differential transform method for solving singular two-point boundary value problems. Several illustrative examples are given to demonstrate the effectiveness of the present method.     

Integral equation methods for elliptic problems with boundary conditions of mixed type

Laplace’s equation with mixed boundary conditions, that is, Dirichlet conditions on parts of the boundary and Neumann conditions on the remaining contiguous parts, is solved on an interior planar domain using an integral equation method. Rapid execution and high accuracy is obtained by combining equations which are of Fredholm’s second kind with compact operators on almost the entire boundary with a recursive compressed inverse preconditioning technique. Then an elastic problem with mixed boundary conditions is formulated and solved in an analogous manner and with similar results. This opens up for the rapid and accurate solution of several elliptic problems of mixed type.     

Electrostatic potential: A new approach for some mixed boundary value problems

A new approach to the solution of problems of electrostatics, some of them with mixed boundary conditions, is presented. The proposed scheme can be used in cases were we have a formal solution in the form of a series in Legendre polynomials and the boundary or matching conditions are given not on the whole interval (0, π) of the polar variable, θ, but only over the interval (0, π/2) or (π/2, π). Truncation of the series after the Nth term and the projection on the subspace generated by the set of the first N even (or odd) Legendre polynomials allows us to determine the unknown coefficients of the approximate solution. The results show rapid convergence toward the exact values as we increase the number of terms, N, included in the approximate solutions. The procedure allows to solve approximately some problems whose exact solutions, we believe, are not yet known.     

A boundary piecewise constant level set method for boundary control of eigenvalue optimization problems

This paper investigates optimization of the least eigenvalue of ?Δ with the constraint of one-dimension Hausdorff measure of Dirichlet boundary. We propose the boundary piecewise constant level set (BPCLS) method based on the regularity technique to combine two types of boundary conditions into a single Robin boundary condition. We derive the first variation of the least eigenvalue w.r.t. the BPCLS function and propose a penalty BPCLS algorithm and an augmented Lagrangian BPCLS algorithm. Numerical results are reported for experiments on ellipse and L-shape domains.     

A path integral formula for certain fourth-order elliptic operators

We define a discretized path integral formula for the operator –²–V. This formula is the generalization of the Feynman-Kac formula for +–V.