Implementation and investigation of iterative solvers in the Discrete Sources Method

The implementation of iterative methods as solvers for the Discrete Sources Method (DSM) is presented. In this method, light scattering computation linear systems with dense and relative small matrices are generated. The linear systems are traditionally solved using the QR-decomposition method. For large particles or particles with extreme geometries even this commonly stable solver can fail. In these cases, we expect that iterative methods can provide a satisfying solution nevertheless.We will present our investigation in two consecutive papers. Here, we study four different iterative solvers (RGMRES, BiCGStab, BiCGStab(l), and MinRes) considering the performance and the accuracy for typical light scattering problems. Using these iterative methods we increased the quality of a solution, especially for oblate spheroids with a higher aspect ratio. Preconditioning technique is considered in the following paper.     

Liouville Integrable System and Associated Integrable Coupling

Liouville integrable discrete integrable system is derived based on discrete isospectral problem. It is shown that the hierarchy is completely integrable in the Liouville sense and possesses bi-Hamiltonian structure. Finally, integrable couplings of the obtained system is given by means of semi-direct sums of Lie algebras.     

三维位势问题边界元法中几乎奇异积分的正则化

采用一种半解析正则化算法,计算了三维位势问题边界元法中近边界点的几乎强奇异和几乎超奇异面积分.该算法适用于三角形线性等参元.对高次单元将其细分为几个三节点三角形单元即可应用该算法.由于几乎奇异性,与内点邻近的单元上的积分,采用半解析正则化积分算法计算;而远处单元的积分仍保持常规高斯积分.对三维热传导算例,计算了近边界点的温度和热流.数值结果证明了该算法的有效性和精确性.     

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

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

普遍化范德华配分函数

对纯流体及其混合物的普遍化范德华配分函数的推导作系统阐述。由该配分函数可导出一些重要的状态方程、混合规则及活度系数模型。还分析了与密度有关的局部组成模型的低密度边界条件，从而可为建立新的局部组成模型及其应用于流体相平衡计算提供理论依据     

An improved boundary element-free method （IBEFM） for two-dimensional potential problems

The interpolating moving least-squares (IMLS) method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method (BEFM), combining the boundary integral equation (BIE) method with the IMLS method, the improved boundary element-free method (IBEFM) for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker δ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker δ function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.     

Symmetry Reduction and Cauchy Problems for a Class of Fourth-Order Evolution Equations

We exploit higher-order conditional symmetry to reduce initial-value problems for evolution equations to Cauchy problems for systems of ordinary differential equations （ODEs）. We classify a class of fourth-order evolution equations which admit certain higher-order generalized conditional symmetries （GCSs） and give some examples to show the main reduction procedure. These reductions cannot be derived within the framework of the standard Lie approach, which hints that the technique presented here is something essential for the dimensional reduction of evolution equations.     

解弹性力学第二类边界积分方程的求积法与分裂外推

利用Sidi奇异求积公式,提出了解曲边多角形域上线性弹性力学第二类边界积分方程的求积法,即离散矩阵的每个元素的生成只需赋值不需计算任何奇异积分.通过估计离散矩阵的特征值和利用Anselone聚紧收敛理论,证明了近似解的收敛性;同时得到了误差的多参数渐近展开式;通过并行地解粗网格上的离散方程,利用分裂外推获得了高精度近似解和后验误差.     

Depth-adaptive regularized reconstruction for reflection diffuse optical tomography

We study reflection diffuse optical tomography using two-dimensional (2D) continuous-wave source-detector arrays on the surface of semi-infinite medium, aiming at imaging the perfusion and the hemoglobin oxygen saturation variation of human cerebral cortex with brain activation. We had previously formulated the inverse problem with Moore-Penrose inversion. When we use simple regularization in this inverse problem, the reconstruction sensitivity decreases markedly with the depth so that the signal in the deep range may be masked by an unwanted signal in the shallow range. In this paper, we propose a depth-adaptive regularized reconstruction, in which we assign a smaller regularization parameter with the depth. We demonstrate improvement of the three-dimensional (3D) reconstruction uniformity using the proposed scheme.     

二维抛物型问题的Strang型交替分段区域分裂格式

给出求解二维抛物型方程的Strang型的交替分段区域分裂格式。交替分段思想可以将区域分为一些不重叠的子区域,Strang型算子分裂技巧通过将高维问题的求解分解为几个低维问题的求解来降低其求解的复杂度。方法是无条件稳定的,理论分析了截断误差。数值算例说明格式的有效性及时空的二阶精度.