On a class of preconditioners for solving the Helmholtz equation

In 1983, a preconditioner was proposed [J. Comput. Phys. 49 (1983) 443] based on the Laplace operator for solving the discrete Helmholtz equation efficiently with CGNR. The preconditioner is especially effective for low wavenumber cases where the linear system is slightly indefinite. Laird [Preconditioned iterative solution of the 2D Helmholtz equation, First Year's Report, St. Hugh's College, Oxford, 2001] proposed a preconditioner where an extra term is added to the Laplace operator. This term is similar to the zeroth order term in the Helmholtz equation but with reversed sign. In this paper, both approaches are further generalized to a new class of preconditioners, the so-called “shifted Laplace” preconditioners of the form Δφ−k²φ with . Numerical experiments for various wavenumbers indicate the effectiveness of the preconditioner. The preconditioner is evaluated in combination with GMRES, Bi-CGSTAB, and CGNR.     

The Hamiltonian Canonical Form for Euler-Lagrange Equations

Based on the theory of calculus of variation, some suffcient conditions are given for some Euler-Lagrangcequations to be equivalently represented by finite or even infinite many Hamiltonian canonical equations. Meanwhile,some further applications for equations such as the KdV equation, MKdV equation, the general linear Euler Lagrangeequation and the cylindric shell equations are given.     

On a fast direct elliptic solver by a modified Fourier method

We describe high order numerical algorithms for the solution of second order elliptic equations in rectangular domains. These algorithms are based on the Fourier method in combination with a subtraction procedure. The singularities at the corner points, arising due to non-smoothness of the boundaries, are treated explicitly using properly constructed singular corner functions. The present algorithm is a generalization of the Fast Poisson Solver developed in our previous paper. This revised version was published online in August 2006 with corrections to the Cover Date.     

Incremental unknowns preconditioning for solving the Helmholtz equation

An efficient preconditioner is developed for solving the Helmholtz problem in both high and low frequency (wavenumber) regimes. The preconditioner is based on hierarchical unknowns on nested grids, known as incremental unknowns (IU). The motivation for the IU preconditioner is provided by an eigenvalue analysis of a simplified Helmholtz problem. The performance of our preconditioner is tested on the iterative solution of two‐dimensional electromagnetic scattering problems. When compared with other well‐known methods, our technique is shown to often provide a better numerical efficacy and, most importantly, to be more robust. Moreover, for the best performance, the number of IU levels used in the preconditioner should be designed for the coarsest grid to have roughly two points per linear wavelength. This result is consistent with the conventional sampling criteria for wave phenomena in contrast with existing IU applications for solving the Laplace/Poisson problem, where the coarsest grid comprises just one interior point. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Scattering from infinite rough tubular surfaces

We study the Helmholtz equation in the exterior of an infinite perturbed cylinder with a Dirichlet boundary condition. Existence and uniqueness of solutions are established using the variational technique introduced (SIAM J. Math. Anal. 2005; 37 (2):598–618). We also provide stability estimates with explicit dependence of the constants in terms of the frequency and the perturbed cylinder thickness. Copyright © 2006 John Wiley & Sons, Ltd.     

圆外区域求解Helmholtz方程

1引言许多科学和工程计算问题都可以归结为无界区域上的偏微分方程边值问题．而求解椭圆方程边值问题的常用技术是有限元方法,可是对于无界区域,在用有限元方法求解时,往往遇到困难．最简单的办法显然是直接略去区域的无界部分求解,但这样做或者导致过低的计算精度,或者要付出很高的计算代价．边界归化,即将求解偏微分方程边值问题转化为边界积分方程,是求解某些无界区域问题的强有力的手段．自70年代以来,有限元和     

A study on the applications of the acoustic design sensitivity analysis of vibrating bodies

The determination of the sensitivity of the acoustical characteristics of vibrating systems with respect to the variation of the design parameters predicting these characteristics is a necessary and important step of the acoustic design and optimization process. Acoustic design sensitivity analysis includes the computation and evaluation of the sensitivity information required for this procedure. In this study, a boundary element code performing the sensitivity analysis of the acoustic pressure by using the matrix sensitivities with respect to different design variables has been developed. The effect of the precision of boundary element discretization on the acoustic pressure sensitivity is examined via this code. The formulation is applied to a multi-source system and the dimension sensitivity analysis of near field pressures of two-dilating-spherical source is performed. The last application is devoted to a real sound source: a washing machine sitting on the floor. Sensitivity of the field pressures to the machine’s dimensions (size), surface velocity and frequency is examined on the bases of the boundary element model of the machine and half-space condition. The impacts of these variables are compared; and a limiting speed for the machine responding both the acoustical and operational requirements is determined.     

Equilibrium theory of two-dimensional simple liquids

Using the Wigner-Kirkwood expansion and bare Lennard-Jones (LJ) (12-6) potential, an effectiveLJ potential is derived, which includes the quantum effects through the expressions of the effective diameter∂(T, λ) and well-depth (T, λ). We use theWCA perturbation theory to calculate the free energy and pressure for theLJ and effectiveLJ potentials. Simple analytic expressions are given for the reference system and the first order correction calculated. The results are quite good at high density. The quantum effects on the free energy and pressure are also discussed.     

ELASTIC BEHAVIOR OF PROTEIN-LIKE SINGLE CHAIN

The conformational properties and elastic behaviors of protein-like single chains in the process of tensile elongation were investigated by means of Monte Carlo method. The sequences of protein-like single chains contain two types of residues: hydrophobic (H) and hydrophilic (P). The average conformations and thermodynamics statistical properties of protein-like single chains with various elongation ratio λ were calculated. It was found that the mean-square end-to-end distance r increases with elongation ratio,λ. The tensor eigenvalues ratio of : decreases with elongation ratio λ for short (HP)x protein-like polymers, however, the ratio of : increases with elongation ratioλ,especially for long (H)x sequence. Average energy per bond increases with elongation ratioλ, especially for(H)x protein-like single chains. Helmholtz free energy per bond also increases with elongation ratioλ. Elastic force (f), energy contribution to force (fU) and entropy contribution to force (fs) for different protein-like single chains were also calculated.These investigations may provide some insights into elastic behaviors of proteins.     

Influence of the tangential velocity on the compressible Kelvin−Helmholtz instability with nonequilibrium effects

Kelvin−Helmholtz (KH) instability is a fundamental fluid instability that widely exists in nature and engineering. To better understand the dynamic process of the KH instability, the influence of the tangential velocity on the compressible KH instability is investigated by using the discrete Boltzmann method based on the nonequilibrium statistical physics. Both hydrodynamic and thermodynamic nonequilibrium (TNE) effects are probed and analyzed. It is found that, on the whole, the global density gradients, the TNE strength and area firstly increase and decrease afterwards. Both the global density gradient and heat flux intensity in the vertical direction are almost constant in the initial stage before a vortex forms. Moreover, with the increase of the tangential velocity, the KH instability evolves faster, hence the global density gradients, the TNE strength and area increase in the initial stage and achieve their peak earlier, and their maxima are higher for a larger tangential velocity. Physically, there are several competitive mechanisms in the evolution of the KH instability. (i) The physical gradients increase and the TNE effects are strengthened as the interface is elongated. The local physical gradients decrease and the local TNE intensity is weakened on account of the dissipation and/or diffusion. (ii) The global heat flux intensity is promoted when the physical gradients increase. As the contact area expands, the heat exchange is enhanced and the global heat flux intensity increases. (iii) The global TNE intensity reduces with the decreasing of physical gradients and increase with the increasing of TNE area. (iv) The nonequilibrium area increases as the fluid interface is elongated and is widened because of the dissipation and/or diffusion.