Blow-up of solution of an initial boundary value problem for a generalized Camassa-Holm equation

In this Letter, we study the following initial boundary value problem for a generalized Camassa-Holm equation

     

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.     

Symmetries of the heat equation on the lattice

Discrete versions of the heat equation on two-dimensional uniform lattices are shown to possess the same symmetry algebra as their continuum limits. Solutions with definite symmetry properties are presented.     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.     

A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.     

Solving the hypersingular boundary integral equation in three-dimensional acoustics using a regularization relationship

Regularization of the hypersingular integral in the normal derivative of the conventional Helmholtz integral equation through a double surface integral method or regularization relationship has been studied. By introducing the new concept of discretized operator matrix, evaluation of the double surface integrals is reduced to calculate the product of two discretized operator matrices. Such a treatment greatly improves the computational efficiency. As the number of frequencies to be computed increases, the computational cost of solving the composite Helmholtz integral equation is comparable to that of solving the conventional Helmholtz integral equation. In this paper, the detailed formulation of the proposed regularization method is presented. The computational efficiency and accuracy of the regularization method are demonstrated for a general class of acoustic radiation and scattering problems. The radiation of a pulsating sphere, an oscillating sphere, and a rigid sphere insonified by a plane acoustic wave are solved using the new method with curvilinear quadrilateral isoparametric elements. It is found that the numerical results rapidly converge to the corresponding analytical solutions as finer meshes are applied.     

A basic study on the sound field analysis of periodic structures by boundary integral equation method

This study deals with the development of the approximate method to analyze the sound field around equally spaced finite obstacles, using the periodic boundary condition. First, on the assumption that the equally spaced finite obstacles are the periodically arranged obstacles, the sound field is analyzed by boundary integral equation method with a Green’s function which satisfies the periodic boundary condition. Furthermore, by comparing these results and the exact solution by using the fundamental solution as Green’s function, the validity of the approximate method is also investigated. Next, in order to evaluate the applicability of the approximate method, the simple formula using some parameters, i.e., the frequency, the period, and the number of obstacles, etc., is proposed. The results of the sound field analysis applied the formula are presented.     

Numerical solution of fractional differential equations via a Volterra integral equation approach

The main focus of this paper is to present a numerical method for the solution of fractional differential equations. In this method, the properties of the Caputo derivative are used to reduce the given fractional differential equation into a Volterra integral equation. The entire domain is divided into several small domains, and by collocating the integral equation at two adjacent points a system of two algebraic equations in two unknowns is obtained. The method is applied to solve linear and nonlinear fractional differential equations. Also the error analysis is presented. Some examples are given and the numerical simulations are also provided to illustrate the effectiveness of the new method.     

Characteristic values of the integral equation satisfied by the Mathieu functions and its application to a system with chirality-pair interaction on a one-dimensional lattice

We investigate low-temperature behaviors of a system with chirality-pair interaction on a one-dimensional lattice. In the course of the investigation, we evaluate asymptotic forms of the characteristic values of the integral equation satisfied by the Mathieu functions. It turns out that the low-temperature behavior of correlation length of the chirality-pair correlation function is different from the one for the Ising model of spin ±1 but akin to the one for the Ising model of infinite spin.     

Random integral equation formulation of a generalized Langevin equation

In this paper, the generalized Langevin equation introduced by Kubo and Mori is formulated as a random integral equation. We consider (1) the existence and uniqueness of the solution, (2) moments of the solution process, (3) a comparison theorem for solution processes, and (4) the Cauchy polygonal approximation to the solution.     

Numerical study of a modified BBGY integral equation for the primitive model of an electrolyte solution

A modified Poisson-Boltzmann equation derived from the BBGY hierarchy of equations is solved numerically for the restricted primitive model of an electrolyte solution. Computations are carried out for 1 : 1 and 2 : 2 electrolytes, and results compared with 1 : 1 Monte Carlo data and results from the Poisson-Boltzmann extensions arising from the Kirkwood hierarchy. Satisfactory agreement with the Monte Carlo results are found in the 1 : 1 case, and the correct qualitative behaviour of thermodynamic properties obtained in the 2 : 2 case. At low concentrations in the 2 : 2 case the modified Poisson-Boltzmann equations from the Kirkwood and BBGY hierarchies predict the anomalous behaviour found in the thermodynamic properties of real electrolytes. The behaviour of the mean electrostatic potential changes from a damped exponential to a damped oscillatory form for 1·2 < κa < 1·3 (where κ is the Debye-Hückel constant and a is an ionic diameter) in the 1 : 1 case and at κa ～ 0·8 in the 2 : 2 case for the parameters considered.     

一维波动方程的解的长期稳定性

对有限区间上的一维波动方程,在各种边界条件下建立了位移函数与初始条件之间的不等式,并由此推论出多种边界条件下的解的长期稳定性．     

边界积分方程法求解目标的声散射T矩阵

提出了计算任意表面形状刚性边界目标散射的基于边界积分方程的T矩阵方法(TMM-BIE).利用Helmholtz积分方程法(HIEM)计算目标表面声场,替代扩展边界法(EBCM)计算中对目标表面声场的近似处理,解决了扩展边界法不能计算任意形状目标的散射T矩阵问题.文中计算了刚性边界的球目标、有限长圆柱目标以及非对称的三维散射体-猫眼(cat's-eye)模型的散射指向性和T矩阵.通过与解析解和HIEM结果比较,证明该方法的有效性.     

带有分数阶热流条件的时间分数阶热波方程及其参数估计问题

研究了Caputo导数定义下带有分数阶热流条件的一维时间分数阶热波方程及其参数估计问题.首先,对正问题给出了解析解;其次,基于参数敏感性分析,利用最小二乘算法同时对分数阶阶数α和热松弛时间τ进行参数估计;最后对不同的热流分布函数所构成的两个初边值问题,分别进行参数估计仿真实验,分析温度真实值和估计值的拟合程度.实验结果表明,最小二乘算法在求解时间分数阶热波方程的两参数估计问题中是有效的.本文为分数阶热波模型的参数估计提供了一种有效的方法.     

An approximation for the boundary optimal control problem of a heat equation defined in a variable domain

In this paper, we consider a numerical approximation for the boundary optimal control problem with the control constraint governed by a heat equation defined in a variable domain. For this variable domain problem, the boundary of the domain is moving and the shape of theboundary is defined by a known time-dependent function. By making use of the Galerkin finite element method, we first project the original optimal control problem into a semi-discrete optimal control problem governed by a system of ordinary differential equations. Then, based on the aforementioned semi-discrete problem, we apply the control parameterization method to obtain an optimal parameter selection problem governed by a lumped parameter system, which can be solved as a nonlinear optimization problem by a Sequential Quadratic Programming （SQP） algorithm. The numerical simulation is given to illustrate the effectiveness of our numerical approximation for the variable domain problem with the finite element method and the control parameterization method.     

Admissibility of transition from the Smolukhovskii integral equation to the Einstein-Fokker-Planck differential equation

For Markov processes we consider limits for the growth rate of the probability density w(t_o, x_o; t, x) of the transition from state x_o at time to to states for which x lies between x and x + dx at time t. These limits follow from the condition of admissibility of transition from the Smolukhovskii integral equation to the Einstein-Fokker-Planck differential equation (EFP). We examine the influence of these limits on the form of the coefficients in EFT. The indicated examination is performed for the prelimit region as well, that is, the region at the end of which we have w(t_o, x_o; t, x) 0.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 10, pp. 99–104, October, 1986.The authors are deeply grateful to D. M. Sedrakyan and A. O. Melikyan for their interest in this work and useful comments.     

Eigenvalues of the one-dimensional Smoluchowski equation

The eigenvalues and eigenfunctions of the Smoluchowski equation are investigated for the case of potentials withN deep wells. The small parameter =kT/V, which measures the ratio of the thermal energy to a typical well depth, is used in connection with the method of matched asymptotic expansion to obtained asymptotic approximations to all the eigenvalues and eigenfunctions. It is found that the eigensolutions fall into two classes, namely (i) the top-of-the-well and (ii) the bottom-of-the-well eigensolutions. The eigenvalues for both classes of solutions are integer multiples of the squqres of the frequencies at the top or bottom of the various wells. The eigenfunctions are, in general, localized to the top or bottom of the corresponding well. The very small eigenvalues require special consideration because the asymptotic analysis is incapable of distinguishing them from the zero eigenvalue with multiplicityN. Another approximation reveals that, in addition to the true zero eigenvalue, there areN-1 eigenvalues of order exp(–). The case of other possible multiple eigenvalues is also examined.     

Numerical solutions of Green''s integral equation for the diffraction of femtosecond laser pulses through a subwavelength aperture

In this letter,we propose a method for the numerical calculations of the femtosecond laser pulse passed through a subwavelength aperture.The time-dependent laser pulse is decomposed into a series of monochromatic simple harmoniv waves.For the light field of the harmonic wave with a single frequency,the numerical calculation is made based on the solution of the Green's integral equation set of the electro-magnetic waves.Such numerical solution is iterated for all the waves with different frequencies,and all thenumerical solutions are transformed into the light fields in the time domain by inverse Fourier transform.The light intensity distributions transmitted the subwavelength aperture are calculated and the resultsshow the propagation of the light field is along the direction of the medium interface.     

On a one-dimensional model for the three-dimensional vorticity equation

The one-dimensional model for the three-dimensional vorticity equation proposed by Constantin, Lax, and Majda is discussed. Some unsatisfactory points are examined, especially when the viscosity is introduced. A different model is suggested, which, while less solvable than the previous one, can be more strictly connected with the three-dimensional vorticity behavior. The study is of interest for the numerical treatment of the three-dimensional vorticity equation.