Using Gaussian Eigenfunctions to Solve Boundary Value Problems

Kernel-based methods are popular in computer graphics, machine learning, and statistics, among other fields; because they do not require meshing of the domain under consideration, higher dimensions and complicated domains can be managed with reasonable effort. Traditionally, the high order of accuracy associated with these methods has been tempered by ill-conditioning, which arises when highly smooth kernels are used to conduct the approximation. Recent advances in representing Gaussians using eigenfunctions have proven successful at avoiding this destabilization in scattered data approximation problems. This paper will extend these techniques to the solution of boundary value problems using collocation. The method of particular solutions will also be considered for elliptic problems, using Gaussian eigenfunctions to stably produce an approximate particular solution.     

On the Initial Boundary Value Problems for the Enskog Equation in Irregular Domains

The paper is concerned with the Enskog equation with a constant high density factor for large initial data in L ¹(R ⁿ). The initial boundary value problem is investigated for bounded domains with irregular boundaries. The proof of an H-theorem for the case of general domains and boundary conditions is given. The main result guarantees the existence of global solutions of boundary value problems for large initial data with all v-moments initially finite and domains having boundary with finite Hausdorff measure and satisfying a cone condition. Existence and uniqueness are first proved for the case of bounded velocities. The solution has finite norm where q = (t ₀, x) is taken on all possible n-dimensional planes Q(v) in R ^n+l intersecting a fixed point and orthogonal to vectors (1, v), v R ⁿ.     

Boundary value problems for Yang—Mills fields

This paper investigates boundary value problems for Hermitian Yang—Mills equations over complex manifolds. The main result is the unique solubility of the Dirichlet problem for the Hermitian Yang—Mills equation. Connections with a number of topics are found, including the link with loop groups.     

A Boussinesq system for two-way propagation of interfacial waves

The theory of internal waves between two layers of immiscible fluids is important both for its applications in oceanography and engineering, and as a source of interesting mathematical model equations that exhibit nonlinearity and dispersion. A Boussinesq system for two-way propagation of interfacial waves in a rigid lid configuration is derived. In most cases, the nonlinearity is quadratic. However, when the square of the depth ratio is close to the density ratio, the coefficients of the quadratic nonlinearities become small and cubic nonlinearities must be considered. The propagation as well as the collision of solitary waves and/or fronts is studied numerically.     

一类各项异性半线性椭圆方程自然边界元与有限元耦合法

将冯康和余德浩提出的自然边界归化方法用于研究一类半线性椭圆方程外区域问题,提出一种自然边界元与有限元的耦合算法、针对某一类半线性椭圆方程,应用变分原理,研究其弱解性及Galerkin逼近,得到有限元解的误差估计及收敛阶O(h^n),最后给出相应数值例子。     

A Gerbe Obstruction to Quantization of Fermions on Odd-Dimensional Manifolds with Boundary

We consider the canonical quantization of fermions on an odd-dimensional manifold with boundary, with respect to a family of elliptic Hermitian boundary conditions for the Dirac Hamiltonian. We show that there is a topological obstruction to a smooth quantization as a function of the boundary conditions. The obstruction is given in terms of a gerbe and its Dixmier–Douady class is evaluated.     

基于等几何分析的边界元法求解Helmholtz问题

将基于一类局部双变量B样条函数的等几何分析方法和Burton-Miller方法相结合,分析三维Helmholtz问题.对于某些从二维参数域映射到三维空间具有奇异点的参数曲面,该方法可以有效地避免奇异点处大量奇异与近奇异积分的计算.数值算例表明该方法具有较好的计算精度和计算效率.复杂问题的分析表明,该方法具有良好的工程应用前景.     

Heat Kernel Asymptotics of Zaremba Boundary Value Problem

The Zaremba boundary-value problem is a boundary value problem for Laplace-type second-order partial differential operators acting on smooth sections of a vector bundle over a smooth compact Riemannian manifold with smooth boundary but with discontinuous boundary conditions, which include Dirichlet boundary conditions on one part of the boundary and Neumann boundary conditions on another part of the boundary. We study the heat kernel asymptotics of Zaremba boundary value problem. The construction of the asymptotic solution of the heat equation is described in detail and the heat kernel is computed explicitly in the leading approximation. Some of the first nontrivial coefficients of the heat kernel asymptotic expansion are computed explicitly. This revised version was published online in July 2006 with corrections to the Cover Date.     

稳态问题混合边界积分方程的高精度求积法与分裂外推

提出了求积法解稳态问题的混合边界积分方程,它拥有高精度,低复杂度.通过并行地解粗网格上的离散方程,根据误差的多参数渐近展开,应用分裂外推算法得到高精度的近似解,同时获得后验误差估计.     

Inverse Eigenvalue Problems for Exploring the Dynamics of Systems Biology Models

This paper describes inverse eigenvalue problems that arise in studying qualitative dynamics in systems biology models. An algorithm based on lift-and-project iterations is proposed, where the lifting step entails solving a constrained matrix inverse eigenvalue problem. In particular, prior to carrying out the iterative steps, $a$-$priori$ bounds on the entries of the Jacobian matrix are computed by relying on the reaction network structure as well as the form of the rate law expressions for the model under consideration. Numerical results on a number of models show that the proposed algorithm can be used to computationally explore the possible dynamical scenarios while identifying the important mechanisms via the use of sparsity-promoting regularization.     

Fourth Order Compact Boundary Value Method for Option Pricing with Jumps

We consider pricing options in a jump-diffusion model which requires solving a partial integro-differential equation. Discretizing the spatial direction with a fourth order compact scheme leads to a linear system of ordinary differential equations. For the temporal direction, we utilize the favorable boundary value methods owing to their advantageous stability properties. In addition, the resulting large sparse system can be solved rapidly by the GMRES method with a circulant Strang-type preconditioner. Numerical results demonstrate the high order accuracy of our scheme and the efficiency of the preconditioned GMRES method.     

二维Helmholtz边界超奇异积分方程解析研究

基于常规边界元法及超奇异边界积分方程复线性耦合的Burton-Miller方法应用于无限域声学问题的最大难点在于处理超奇异积分（二维问题）.目前,此类超奇异积分主要使用各种弱奇异/正则化方法求解,而这些弱奇异/正则化方法具有时间消耗大等弱点.基于围道积分定理,本文给出一种使用常值单元的二维Helmholtz边界超奇异积分的解析表达式.在有限部分积分意义下,所有的奇异和超奇异积分可以解析表达.数值算例表明该解析表达式是有效的.     

The Numerical Methods for Oscillating Singularities in Elliptic Boundary Value Problems

The singularities near the crack tips of homogeneous materials are monotone of type r^α and r^α log^δr (depending on the boundary conditions along nonsmooth domains). However, the singularities around the interfacial cracks of the heterogeneous bimaterials are oscillatory of type r^α sin( log r). The method of auxiliary mapping (MAM), introduced by Babu ka and Oh, was proven to be successful in dealing with r^α type singularities. However, the effectiveness of MAM is reduced in handling oscillating singularities. This paper deals with oscillating singularities as well as the monotone singularities by extending MAM through introducing the power auxiliary mapping and the exponential auxiliary mapping.     

Boundary element-free method for elastodynamics

The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.     

Lower-Dimensional Volumes and Kastler-Kalau-Walze Type Theorem for Manifolds with Boundary

In this paper, we define lower-dimensional volumes of spin manifolds with boundary. We compute the lower-dimensional volume Vol（2＇2） for 5-dimensional and 6-dimensional spin manifolds with boundary and we also get the Kastler-Kalau-Walze type theorem in this case.     

大规模声学边界元法的GPU并行计算

提出一种大规模声学边界元法的高效率、高精度GPU并行计算方法.基于Burton-Miller边界积分方程,推导适于GPU的并行计算格式并实现了传统边界元法的GPU加速算法.为提高原型算法的效率,研究GPU数据缓存优化方法.由于GPU的双精度浮点运算能力较低,为了降低数值误差,研究基于单精度浮点运算实现的doublesingle精度算法.数值算例表明,改进的算法实现了最高89.8%的GPU使用效率,且数值精度与直接使用双精度数相当,而计算时间仅为其1/28,显存消耗也仅为其一半.该方法可在普通PC机(8GB内存,NVIDIA Ge Force 660 Ti显卡)上快速完成自由度超过300万的大规模声学边界元分析,计算速度和内存消耗均优于快速边界元法.     

Boundary element-free method for elastodynamics

1 Introduction In recent years, more and more attention has been paid to researches on the meshless (or meshfree) method, which makes it a hot direction of computational mechanics[1,2]. The meshless method is the approximation based on nodes, then the large deformation and crack growth problems can be simulated with the method without the re-meshing technique. And the meshless method has some advantages over the traditional computa- tional methods, such as finite element method (FEM) and boun…     

Boundary value problems for the steady Boltzmann equation

We discuss steady boundary value problems for the Boltzmann equation with inflow and diffusive boundary conditions in one, two, and three dimensions, with suitable truncations of the collision kernel. General existence and uniqueness results are obtained if the domain is sufficiently small. In one dimension, the existence of solutions on general intervals is obtained by abstract fixed-point theory.     

The Solvability of the Discrete Boundary Value Problem on the Half-Line

This paper provides conditions for the existence of a solution to the second-order nonlinear boundary value problem on the half-line of the form with αx(0)+βa(0)Δx(0)=0,x(∞)=d, where d,α,β∈R, α2+β2>0. To achieve our goal, we use Schauder’s fixed-point theorem and the perturbation technique for a Fredholm operator of index 0. Moreover, we construct the necessary condition for the existence of a solution to the considered problem.