Method of fundamental solutions with optimal regularization techniques for the Cauchy problem of the Laplace equation with singular points

The purpose of this study is to propose a high-accuracy and fast numerical method for the Cauchy problem of the Laplace equation. Our problem is directly discretized by the method of fundamental solutions (MFS). The Tikhonov regularization method stabilizes a numerical solution of the problem for given Cauchy data with high noises. The accuracy of the numerical solution depends on a regularization parameter of the Tikhonov regularization technique and some parameters of the MFS. The L-curve determines a suitable regularization parameter for obtaining an accurate solution. Numerical experiments show that such a suitable regularization parameter coincides with the optimal one. Moreover, a better choice of the parameters of the MFS is numerically observed. It is noteworthy that a problem whose solution has singular points can successfully be solved. It is concluded that the numerical method proposed in this paper is effective for a problem with an irregular domain, singular points, and the Cauchy data with high noises.     

The scattering problem in quantum mechanics as an eigenvalue problem

A method is proposed which allows the scattering problem to reduce to the eigenvalue problem. Unlike the usual method when the scattering phase is extracted from the asymptotics of solution of the Cauchy problem at a given collision energy, in the proposed method the collision energy is obtained from the solution of the Sturm-Liouville problem at a given scattering phase. The continuous analog of the Newton method is used for the numerical realization of the proposed method.     

An adaptive greedy technique for inverse boundary determination problem

In this paper, the method of fundamental solutions (MFS) is employed for determining an unknown portion of the boundary from the Cauchy data specified on parts of the boundary. We propose a new numerical method with adaptive placement of source points in the MFS to solve the inverse boundary determination problem. Since the MFS source points placement here is not trivial due to the unknown boundary, we employ an adaptive technique to choose a sub-optimal arrangement of source points on various fictitious boundaries. Afterwards, the standard Tikhonov regularization method is used to solve ill-conditional matrix equation, while the regularization parameter is chosen by the L-curve criterion. The numerical studies of both open and closed fictitious boundaries are considered. It is shown that the proposed method is effective and stable even for data with relatively high noise levels.     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

共轭梯度法求解双曲传热多宗量反问题

提出双曲传热反问题热物性参数和边界条件多宗量联合反演的一般数值求解模式,考虑了非均质和分布参数的影响,时域上采用时域精细算法进行离散,建立了便于敏度分析的有限元正演模型.由最小二乘原理建立反演模型,应用共轭梯度法进行求解.探讨了时间步长和测量误差对反演结果的影响,并进行了数值验证.     

A fast method for acoustic imaging of multiple three-dimensional objects

This paper is concerned with the inverse problem for imaging multiple three-dimensional objects using the information of the far-field pattern of the scattered wave. A spatially dependent function, which has noticeably different values inside and outside the obstacle, is derived. A numerical method based on the characterization is developed to obtain a visualization of the obstacle. The most remarkable advantage of this method is that it does not need any prior knowledge about the geometry and physical properties of the scatterer, and requires only the information of the far-field measurements for a finite number of directions of incidence and observation distributed over a limited range. Furthermore, the scheme is very simple and fast since it avoids the use of the iterative procedure and requires only the solution of a linear system. Some numerical examples with synthetic far-field data are given showing the practicality and efficiency of this scheme.     

大型粮仓温度分布反演数学模型的研究

采用有限元数值模拟方法,以响应面模型RSM(Response surface model)理论为基础对大型粮仓冬、夏两季的温度分布进行了深入分析,推导了粮仓温度反演数学模型。将采集的温度值导入RSM模型中,并结合单因素轮换法、正交实验和统计方法,通过构造一个具有明确表达形式的多项式来预测粮仓温度分布反演数学模型。在数学模型的基础上给出各因子交互作用下的设计空间、各因素对粮温影响的主效应及交互效应。实验证明该数学模型误差小,精度高,具有良好的计算性能,在大型粮仓温度监控系统中采用该模型能够解决布线少、预测空间范围广的问题。     

A Modified Kernel Method for Solving Cauchy Problem of Two-Dimensional Heat Conduction Equation

In this paper, a Cauchy problem of two-dimensional heat conduction equation is investigated. This is a severely ill-posed problem. Based on the solution of Cauchy problem of two-dimensional heat conduction equation, we propose to solve this problem by modifying the kernel, which generates a well-posed problem. Error estimates between the exact solution and the regularized solution are given. We provide a numerical experiment to illustrate the main results.     

A temperature-related boundary Cauchy–Born method for multi-scale modeling of silicon nano-structures

The surface, edge and corner effects have significant influences in the electrical and optical properties of silicon nano-structures. In this paper, a novel hierarchical temperature-related multi-scale model is presented based on the boundary Cauchy–Born method to investigate not only the surface but also the edge and corner effects in thermal properties of diamond-like structures such as silicon nano-structures at finite temperature. A combined finite element method and molecular dynamics are respectively employed in macro- and micro-scale levels. The temperature-related Cauchy–Born rule is applied using the Helmholtz free energy, as the energy density of equivalent continua relating to the Tersoff inter-atomic potential. The model employs radial quadratures at the surface, edge and corner elements as an indicator of material behavior. The capability of computational algorithm is illustrated by numerical simulation of a nano-scale cube at finite temperature and the results are compared with the atomistic model.     

Inversion of Robin coefficient by a spectral stochastic finite element approach

This paper investigates a variational approach to the nonlinear stochastic inverse problem of probabilistically calibrating the Robin coefficient from boundary measurements for the steady-state heat conduction. The problem is formulated into an optimization problem, and mathematical properties relevant to its numerical computations are investigated. The spectral stochastic finite element method using polynomial chaos is utilized for the discretization of the optimization problem, and its convergence is analyzed. The nonlinear conjugate gradient method is derived for the optimization system. Numerical results for several two-dimensional problems are presented to illustrate the accuracy and efficiency of the stochastic finite element method.     

等离子体平衡逆问题的一种数值解法

本文借助离散化算子,把滤波方法推广到非圆截面,对具有非圆截面的等离子体平衡逆问题建立了一个稳定的显式逆推滤波算法;给出了部分非圆截面的计算结果;对圆截面情形的数值解与解析解作了比较。     

Numerical solution of 3D Green's function for the dynamic system of anisotropic elasticity

In this Letter, we used homotopy perturbation method to obtain numerical solution of the 3D Green's function for the dynamic system of anisotropic elasticity. Application of homotopy perturbation method to this problem shows the rapid convergence of the sequence constructed by this method to the exact solution. The numerical results obtained from convolution of Green's function and data of the Cauchy problem are compared with the exact solution for cubic media. The results reveal that the proposed method is very effective and simple.     

Discretization and its proof for numerical solution of a Cauchy problem for Laplace equation with inaccurately given Cauchy conditions on an inaccurately defined arbitrary surface

Cauchy problem for the Laplace equation with inaccurately given Cauchy conditions on an inaccurately defined arbitrary surface is considered. Discretization was performed and proved to obtain a numerical solution. An economic algorithm is proposed. The text was submitted by the authors in English.     

Exact solutions of semiclassical non-characteristic Cauchy problems for the sine-Gordon equation

The use of the sine-Gordon equation as a model of magnetic flux propagation in Josephson junctions motivates studying the initial-value problem for this equation in the semiclassical limit in which the dispersion parameter ε tends to zero. Assuming natural initial data having the profile of a moving −2π kink at time zero, we analytically calculate the scattering data of this completely integrable Cauchy problem for all ε>0 sufficiently small, and further we invert the scattering transform to calculate the solution for a sequence of arbitrarily small ε. This sequence of exact solutions is analogous to that of the well-known N-soliton (or higher-order soliton) solutions of the focusing nonlinear Schrödinger equation. We then use plots obtained from a careful numerical implementation of the inverse-scattering algorithm for reflectionless potentials to study the asymptotic behavior of solutions in the semiclassical limit. In the limit ε↓0 one observes the appearance of nonlinear caustics, i.e. curves in space-time that are independent of ε but vary with the initial data and that separate regions in which the solution is expected to have different numbers of nonlinear phases.In the appendices, we give a self-contained account of the Cauchy problem from the perspectives of both inverse scattering and classical analysis (Picard iteration). Specifically, Appendix A contains a complete formulation of the inverse-scattering method for generic L¹-Sobolev initial data, and Appendix B establishes the well-posedness for L^p-Sobolev initial data (which in particular completely justifies the inverse-scattering analysis in Appendix A).     

Central potential reconstruction from scattering data by the variable boundary method with the use of Bargmann potentials

An effective numerical method for solving the inverse quantum scattering problem with a fixed angular momentum is suggested. The method is based on combining the variable boundary method suggested earlier and the method of Bargmann potentials. Test calculations are evidence of high method effectiveness in various particular cases, for instance, when there are bound states and resonances.     

Comparison of solvers for the generalized Riemann problem for hyperbolic systems with source terms

We compare four different approximate solvers for the generalized Riemann problem (GRP) for non-linear systems of hyperbolic equations with source terms. The GRP is a special Cauchy problem for a hyperbolic system with source terms whose initial condition is piecewise smooth. We briefly review the four solvers currently available and carry out a systematic assessment of these in terms of accuracy and computational efficiency. These solvers are the building block for constructing high-order numerical schemes of the ADER type for solving the general initial-boundary value problem for inhomogeneous systems in multiple space dimensions, in the frameworks of finite volume and discontinuous Galerkin finite element methods.     

An inverse scattering problem from an impedance obstacle

In this paper we consider an inverse scattering problem from an obstacle with impedance boundary condition. Our aim is to recover the unknown scatterer from the far field pattern iteratively assuming the impedance function. Our method, while remaining in the framework of Newton’s method, based on a system of two nonlinear integral equations which is equivalent to the original inverse problem, avoids the need of calculating a direct problem at each iteration. Because of the ill-posedness of this problem, regularization method for example, Tikhonov regularization, is incorporated in our solution scheme. Several numerical examples with only one incident wave are given at the end of the paper to show the feasibility of our method.     

Propagator of a Charged Particle with a Spin in Uniform Magnetic and Perpendicular Electric Fields

We construct an explicit solution of the Cauchy initial value problem for the time-dependent Schrödinger equation for a charged particle with a spin moving in a uniform magnetic field and a perpendicular electric field varying with time. The corresponding Green function (propagator) is given in terms of elementary functions and certain integrals of the fields with a characteristic function, which should be found as an analytic or numerical solution of the equation of motion for the classical oscillator with a time-dependent frequency. We discuss a particular solution of a related nonlinear Schrödinger equation and some special and limiting cases are outlined.     

间断解问题的有限体积法

有限体积法（FVM）是集有限差分法和有限元方法之优点而发展起来的一种新的数值方法,首先利用守恒问题给出FVM的基本概念和基本过程,接下来就间断解问题的求解进一步介绍了FVM的一些比较活跃的构造方法,并给出了用其中几个方法模拟二维溃坝问题的数值结果,然后按问题分类介绍了FVM数值方法和实际应用的进展,最后给出了FVM理论分析的一些进展。     

Inverse radiation problem in one-dimensional slab by time-resolved reflected and transmitted signals

A time-domain inverse approach is proposed for estimating the distribution of absorbing and scattering coefficients in one-dimensional inhomogeneous media. The temporal reflected and transmitted signals are detected when an ultra-short pulse irradiates on the boundary of semi-transparent scattering media. Forward computation and inverse algorithm employ the least-squares finite element method and conjugate gradient method, respectively. As the prevalent diffusion approximation is not employed in our model, the present approach can be extended to more comprehensive application. The investigation about detected signals indicates that the reflected signals play a significant role in reconstructing optical properties; the signals in early sampling time are more important than those at long-time logarithm slope, and so, more attention should be paid to the early signals in the solution of inverse radiation problem. Three different inverse radiation problems are investigated to show the ability of the present approach to deal with the two-layer, three-layer and continuous inhomogeneous media. The effect of measured errors on the accuracy of reconstruction is investigated by adding artificial random errors. The results indicate that accurate reconstruction depends on not only precise numerical simulation but also quality of detected data.