A fictitious domain decomposition method for the solution of partially axisymmetric acoustic scattering problems. Part 2: Neumann boundary conditions

We present a fictitious domain decomposition method for the fast solution of acoustic scattering problems characterized by a partially axisymmetric sound‐hard scatterer. We apply this method to the solution of a mock‐up submarine problem, and highlight its computational advantages and intrinsic parallelism. A key component of our method is an original idea for addressing a Neumann boundary condition in the general framework of a fictitious domain method. This idea is applicable to many other linear partial differential equations besides the Helmholtz equation. Copyright © 2003 John Wiley & Sons, Ltd.     

Compact difference method for solving the fractional reaction–subdiffusion equation with Neumann boundary value condition

In this paper, we derive a high-order compact finite difference scheme for solving the reaction–subdiffusion equation with Neumann boundary value condition. The L1 method is used to approximate the temporal Caputo derivative, and the compact difference operator is applied for spatial discretization. We prove that the compact finite difference method is unconditionally stable and convergent with order O(τ^2?α+h⁴) in L₂ norm, where τ, α, and h are the temporal step size, the order of time fractional derivative and the spatial step size, respectively. Finally, some numerical experiments are carried out to show the effectiveness of the proposed difference scheme.     

Quintic nonpolynomial spline method for the solution of a second-order boundary-value problem with engineering applications

Nonpolynomial quintic spline functions are used to develop a numerical algorithm for computing an approximation to the solution of a system of second order boundary value problems associated with heat transfer. We show that the approximate solutions obtained by our algorithm are better than those produced by other spline and domain decomposition methods. A comparison of our algorithm with nonpolynomial quadratic spline method is discussed with the help of two numerical examples.     

无限维两体量子系统可分态的一个必要充分条件

构造所考虑的量子系统与辅助量子系统张量积空间中的经典态,利用经典态在辅助复合系统的偏迹运算,研究了无限维两体复合量子系统态纠缠的识别问题,得出了无限维两体复合量子系统态可分的一个必要充分条件,推广了Li Nan和Luo Shunlong的有限维复合量子系统的相关结果.对于识别无限维量子态的纠缠或可分是一个补充.     

零边界条件下二维元胞自动机矩阵可逆性分析

该文介绍了二维元胞自动机的典型邻域结构,对多状态二维元胞自动机的矩阵描述进行了拓展.通过对零边界条件下元胞自动机的状态转移矩阵的块矩阵的分析,得出矩阵的可逆性与元胞自动机冯诺依曼型邻域之间的联系.     

FETI‐DP,BDDC, and block Cholesky methods

The FETI‐DP and BDDC algorithms are reformulated using Block Cholesky factorizations, an approach which can provide a useful framework for the design of domain decomposition algorithms for solving symmetric positive definite linear system of equations. Instead of introducing Lagrange multipliers to enforce the coarse level, primal continuity constraints in these algorithms, a change of variables is used such that each primal constraint corresponds to an explicit degree of freedom. With the new formulation of these algorithms, a simplified proof is provided that the spectra of a pair of FETI‐DP and BDDC algorithms, with the same set of primal constraints, are essentially the same. Numerical experiments for a two‐dimensional Laplace's equation also confirm this result. Copyright © 2005 John Wiley & Sons, Ltd.     

A Gaussian RBFs method with regularization for the numerical solution of inverse heat conduction problems

In this paper, a recursion numerical technique is considered to solve the inverse heat conduction problems, with an unknown time-dependent heat source and the Neumann boundary conditions. The numerical solutions of the heat diffusion equations are constructed using the Gaussian radial basis functions. The details of algorithms in the one-dimensional and two-dimensional cases, involving the global or partial initial conditions, are proposed, respectively. The Tikhonov regularization method, with the generalized cross-validation criterion, is used to obtain more stable numerical results, since the linear systems are badly ill-conditioned. Moreover, we propose some results of the condition number estimates to a class of positive define matrices constructed by the Gaussian radial basis functions. Some numerical experiments are given to show that the presented schemes are favourably accurate and effective.     

基于Neumann随机有限元法的复合材料压力容器的应力分析

利用Neumann级数展开的Monte-Carlo随机有限元法，分析了纤维缠绕壳体在内压及材料工程常数随机变化作用下的应力响应，绘出了最大应力响应分布曲线，为该类结构可靠性评估及可靠性设计提供了应力响应数字依据。     

Equilibrium on line method (ELM) for imposition of Neumann boundary conditions in the finite point method (FPM)

Equilibrium on line method (ELM) for imposition of Neumann boundary conditions in the finite point method (FPM) is presented. In contrary to weak‐form‐based methods, strong‐form‐based methods such as the FPM are often unstable and less accurate, especially for problems governed by partial differential equations with Neumann (derivative) boundary conditions. In this paper, a truly meshless approach for imposition of Neumann boundary conditions in the FPM is proposed and adopted for 2D elasticity analyses. In the proposed method, equilibrium on lines on the Neumann boundary conditions is satisfied as Neumann boundary condition equations. Numerical studies show that this method for imposition of Neumann boundary is simple to implement and computationally efficient and also leads to more stable and accurate results. Copyright © 2006 John Wiley & Sons, Ltd.     

Nuclear polarization and entanglement in spin systems

We investigated relationships between entanglement measures and the order parameter (nuclear polarization) in nuclear spin systems controlled by the nuclear magnetic resonance (NMR) technique. Since spin polarization can be easily manipulated by the NMR technique, experimentalists are presented with an opportunity to study the dynamic properties of entanglement, i.e., the creation and evolution of entangled states. Our approach may constitute the basis for researching the relations between the entanglement measures and measurable parameters of order in other quantum systems. An erratum to this article can be found at