基于同伦方法反演非饱和土中镉离子传输参数

依据实验室土柱实验结果，引入同伦方法对镉离子在不同非饱和土样中的传输参数进行了反演计算，得到了相应的反演结果．     

Iterated Tikhonov regularization with a general penalty term

Tikhonov regularization is one of the most popular approaches to solving linear discrete ill‐posed problems. The choice of the regularization matrix may significantly affect the quality of the computed solution. When the regularization matrix is the identity, iterated Tikhonov regularization can yield computed approximate solutions of higher quality than (standard) Tikhonov regularization. This paper provides an analysis of iterated Tikhonov regularization with a regularization matrix different from the identity. Computed examples illustrate the performance of this method.     

A coupled complex boundary method for an inverse conductivity problem with one measurement

We recently proposed in [Cheng, XL et al. A novel coupled complex boundary method for inverse source problems Inverse Problem 2014 30 055002] a coupled complex boundary method (CCBM) for inverse source problems. In this paper, we apply the CCBM to inverse conductivity problems (ICPs) with one measurement. In the ICP, the diffusion coefficient q is to be determined from both Dirichlet and Neumann boundary data. With the CCBM, q is sought such that the imaginary part of the solution of a forward Robin boundary value problem vanishes in the problem domain. This brings in advantages on robustness and computation in reconstruction. Based on the complex forward problem, the Tikhonov regularization is used for a stable reconstruction. Some theoretical analysis is given on the optimization models. Several numerical examples are provided to show the feasibility and usefulness of the CCBM for the ICP. It is illustrated that as long as all the subdomains share some portion of the boundary, our CCBM-based Tikhonov regularization method can reconstruct the diffusion parameters stably and effectively.     

时间分数次扩散方程反演源项问题的迭代正则化方法

时间分数次扩散方程中反演源项问题是一类经典不适定问题.本文构造了一种新的迭代格式作为正则化方法,给出了先验和后验参数选取下相应的收敛性分析.数值算例验证该方法的有效性.     

Multiresolution Symplectic Scheme for wave propagation in complex media

A fast adaptive symplectic algorithm named Multiresolution Symplectic Scheme (MSS) was first presented to solve the problem of the wave propagation (WP) in complex media, using the symplectic scheme and Daubechies‘ compactly supported orthogonal wavelet transform to respectively discretise the time and space dimension of wave equation. The problem was solved in multiresolution symplectic geometry space under the conservative Hamiltonian system rather than the traditional Lagrange system. Due to the fascinating properties of the wavelets and symplectic scheme, MSS is a promising method because of little computational burden, robustness and reality of long-time simulation.     

EIGENVALUE PROBLEM OF A LARGE SCALE INDEFINITE GYROSCOPIC DYNAMIC SYSTEM

Gyroscopic dynamic system can be introduced to Hamiltonian system.Based on an adjoint symplectic subspace iteration method of Hamiltonian gyroscopic system, an adjoint symplectic subspace iteration method of indefinite Hamiltonian function gy- roscopic system was proposed to solve the eigenvalue problem of indefinite Hamiltonian function gyroscopic system.The character that the eigenvalues of Hamiltonian gyroscopic system are only pure imaginary or zero was used.The eigenvalues that Hamiltonian function is negative can be separated so that the eigenvalue problem of positive definite Hamiltonian function system was presented,and an adjoint symplectic subspace iteration method of positive definite Hamiltonian function system was used to solve the separated eigenvalue problem.Therefore,the eigenvalue problem of indefinite Hamiltonian function gyroscopic system was solved,and two numerical examples were given to demonstrate that the eigensolutions converge exactly.     

大型辛矩阵本征问题的逆迭代法

基于共轭辛子空间迭代法,求解了大型辛矩阵的主要本征解。随着迭代的进行,可以无限地逼近其精确解。     

Determination of the molecular weight distribution of linear polymers by inversion of a blending law on complex viscosities

The method presented in this paper allows to calculate the molecular weight distribution (MWD) of linear homopolymer melts from the complex shear modulus data measured in a wide frequency domain. An empirical blending law on complex viscosities is first developed; as a consequence, the variations of the storage and loss modulus as a function of MWD are presented. This simulation demonstrates also the role of the shape of the MWD itself, and shows that one should not postulate a priori the shape of the MWD. An efficient numerical approach based on a Tikhonov regularization method with constraint is used to solve this ill-posed problem; the MWD is hence derived without any assumption on its shape. This method is first applied on simulated data to prove its numerical efficiency. Then the inversion method is applied on complex moduli data of various commercial polymers (polypropylene, polyethylene and polystyrene) and on an artificial mixture of polystyrene that have been presented in the literature. For amorphous polymers, the coupling of the terminal relaxation domains with the transition region at higher frequency leads to errors in the low molecular weight tail: one way to solve this problem is to cut off the experimental data at the high frequencies. This general method needs only a few physical parameters, namely the scaling law for the Newtonian viscosity η₀=f(M _w ) and the plateau modulus G _N ⁰, and leads to reasonable results with respect to the simplicity of the viscoelastic model used. Received: 27 October 1997 Accepted: 24 February 1998     

四角点支承四边自由矩形薄板屈曲问题的新解析解

角点支承矩形薄板的屈曲问题是板壳力学的一类重要课题，控制方程和边界条件的复杂性导致寻求该类问题的解析解十分困难。虽然各类近似/数值方法可用于解决此类难题，但作为基准的精确解析解在公开文献中鲜有报道。本文基于近年来提出的辛叠加方法，解析求解了四角点支承四边自由矩形薄板的屈曲问题。首先将问题拆分为两个子问题，接着利用分离变量与辛本征展开推导出子问题的解析解，最后通过叠加获得原问题的解。由于求解过程从基本控制方程出发，逐步严格推导，无需假定解的形式，因此本文解法是一种理性的解析方法。数值算例给出了不同长宽比和不同面内载荷比情况下，四角点支承四边自由矩形薄板的屈曲载荷和典型屈曲模态，并经有限元方法验证，确认了解析解的正确性。     

四边任意支承条件下弹性矩形薄板弯曲问题的解析解

利用辛几何法推导出了四边为任意支承条件下矩形薄板弯曲的解析解。在分析过程中首先把矩形薄板弯曲问题表示成Hamilton正则方程，然后利用辛几何方法对全状态相变量进行分离变量，求出其本征值后，再按本征函数展开的方法求出四边为任意支承条件下矩形薄板弯曲的解析解。由于在求解过程中并不需要人为的事先选取挠度函数，而是从弹性矩形薄板弯曲的基本方程出发，直接利用数学的方法求出问题的解析解，使得这类问题的求解更加理论化和合理化。文中的最后还给出了计算实例来验证本文方法的正确性。