A Solver for Helmholtz System Generated by the Discretization of Wave Shape Functions

An interesting discretization method for Helmholtz equations was introduced in B. Després [1]. This method is based on the ultra weak variational formulation (UWVF) and the wave shape functions, which are exact solutions of the governing Helmholtz equation. In this paper we are concerned with fast solver for the system generated by the method in [1]. We propose a new preconditioner for such system, which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in [13]. In our numerical experiments, this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations, and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.     

Boltzmann's Formulation of the Helmholtz Vorticity Theorem

A dual form of the general relativistic and covariant generalization of the first vorticity theorem of Helmholtz is proven. This dual form is the generalization of Boltzmann's formulation of the Helmholtz theorem.     

On the indefinite Helmholtz equation: Complex stretched absorbing boundary layers,iterative analysis,and preconditioning

This paper studies and analyzes a preconditioned Krylov solver for Helmholtz problems that are formulated with absorbing boundary layers based on complex coordinate stretching. The preconditioner problem is a Helmholtz problem where not only the coordinates in the absorbing layer have an imaginary part, but also the coordinates in the interior region. This results into a preconditioner problem that is invertible with a multigrid cycle. We give a numerical analysis based on the eigenvalues and evaluate the performance with several numerical experiments. The method is an alternative to the complex shifted Laplacian and it gives a comparable performance for the studied model problems.     

Data completion method for the characterization of sound sources

This paper presents a different approach to solve the inverse acoustic problem. This problem is an "ill-posed" problem since the solution is very sensitive to measurement precision. A classical way to solve this problem consists in inversing a propagation operator which relates structure quantities (acoustic pressures or gradients) to near-field quantities (acoustic pressures or gradients). This can be achieved by using near-field acoustical holography (NAH) in separable coordinate systems. In order to overcome this limitation, the inverse boundary element method (IBEM) can be implemented to recover all acoustic quantities in a three-dimensional space and on an arbitrary three-dimensional source surface. In this paper, the data completion method (DCM) is developed: the acoustic gradients and pressures are known on a surface surrounding the source, but are unknown on its structure. The solution is given by the resolution of the Helmholtz formulation applied on the empty domain between the two boundaries made by the measurements quantities and the structure of the source. The conventional method applies directly the integral formulation for the empty domain. Another way of solving this Helmholtz formulation can be achieved by splitting it in two well-posed subproblems in a Steklov-Poincare?'s formulation. The data completion method allows one to solve the problem with acoustic perturbations due to sources on the exterior domain, or due to a confined domain, without altering the results.     

Spin connection resonance in magnetic motors

A mechanism is proposed for rotation of magnetic assemblies by a torque consisting of the magnetic dipole moment of the assembly and a magnetic field generated from space–time in Einstein–Cartan–Evans (ECE) field theory. It is shown that when the magnetic assembly is stationary, the space–time is described by a Helmholtz wave equation in the tetrad as eigenfunction. This is a balance condition in which the Cartan torsion of the space–time is zero, but in which the tetrad and spin connection are non-zero. This balance may be broken by a driving current density produced by the magnetic assembly. The Helmholtz equation becomes an undamped oscillator equation. At resonance the torque on the magnetic assembly may be amplified sufficiently to cause the whole assembly to rotate, as observed experimentally in a repeatable and reproducible manner.     

Bracket formulation for irreversible classical fields

A bracket formulation for irreversible fields analogous to that for hamiltonian fields is presented. The formulation contains a bracket with symmetric and antisymmetric components and a generator of time translation. Plasma examples are given when the generator of time translation is the energy, entropy and Helmholtz free energy.     

A Novel Three-Dimensional Wide-Angle Beam Propagation Method Based on Split-Step Fast Fourier Transform

A novel three-dimensional wide-angle beam propagation method based on the split-step fast Fourier transform is developed. The formulation is based on the three-dimensional Helmholtz wave equation. Each propagation step is performed by utilizing both the FFT and split-step scheme. The solution of Helmholtz wave equation does not make the slowly varying envelope and one-way propagation approximations. To validate the efficiency and accuracy, numerical results for a propagation beam in a tilted step-index optical waveguide are compared with other beam propagation algorithms.     

用多极理论分析圆桩对称微波谐振腔

用多极理论计算具有复杂几何形状、圆柱对称微波腔的谐振频率,推导出用多极理论计算圆柱对称微波腔谐振频率的本征值方程。三个工程实例的计算结果表明,用多极理论计算圆柱对称微波腔谐振频率,不仅具有较高的计算精度,而且可以很方便地应用于复杂几何形状圆柱对称微波腔工程问题的设计与计算,多极理论是计算圆柱对称微波腔谐振频率的一种有效方法。     

Solving the hypersingular boundary integral equation in three-dimensional acoustics using a regularization relationship

Regularization of the hypersingular integral in the normal derivative of the conventional Helmholtz integral equation through a double surface integral method or regularization relationship has been studied. By introducing the new concept of discretized operator matrix, evaluation of the double surface integrals is reduced to calculate the product of two discretized operator matrices. Such a treatment greatly improves the computational efficiency. As the number of frequencies to be computed increases, the computational cost of solving the composite Helmholtz integral equation is comparable to that of solving the conventional Helmholtz integral equation. In this paper, the detailed formulation of the proposed regularization method is presented. The computational efficiency and accuracy of the regularization method are demonstrated for a general class of acoustic radiation and scattering problems. The radiation of a pulsating sphere, an oscillating sphere, and a rigid sphere insonified by a plane acoustic wave are solved using the new method with curvilinear quadrilateral isoparametric elements. It is found that the numerical results rapidly converge to the corresponding analytical solutions as finer meshes are applied.     

On the modeling of narrow gaps using the standard boundary element method

Numerical methods based on the Helmholtz integral equation are well suited for solving acoustic scattering and diffraction problems at relatively low frequencies. However, it is well known that the standard method becomes degenerate if the objects that disturb the sound field are very thin. This paper makes use of a standard axisymmetric Helmholtz integral equation formulation and its boundary element method (BEM) implementation to study the behavior of the method on two test cases: a thin rigid disk of variable thickness and two rigid cylinders separated by a gap of variable width. Both problems give rise to the same kind of degeneracy in the method, and modified formulations have been proposed to overcome this difficulty. However, such techniques are better suited for the so-called thin-body problem than for the reciprocal narrow-gap problem, and only the first is usually dealt with in the literature. A simple integration technique that can extend the range of thicknesses/widths tractable by the otherwise unmodified standard formulation is presented and tested. This technique is valid for both cases. The modeling of acoustic transducers like sound intensity probes and condenser microphones has motivated this work, although the proposed technique has a wider range of applications.     

A PROPOSAL OF QUANTIZATION OF SYSTEMS WITHOUT STANDARD LAGRANGIANS

For many classical systems which disobey Helmholtz conditions and thus possess no stan-dard Lagrangian formulation on velocity phase space TQ, ordinary quantization procedures cannot be directly applied. By shifting to an enlarged space TTQ and viewing the systems as typically constrained ones, we propose a quantization method which seems to be practicable as illustrated by the examples.     

用多极理论分析圆柱对称微波谐振腔

 用多极理论计算具有复杂几何形状、圆柱对称微波腔的谐振频率,推导出用多极理论计算圆柱对称微波腔谐振频率的本征值方程。三个工程实例的计算结果表明,用多极理论计算圆柱对称微波腔谐振频率,不仅具有较高的计算精度, 而且可以很方便地应用于复杂几何形状圆柱对称微波腔工程问题的设计与计算,多极理论是计算圆柱对称微波腔谐振频率的一种有效方法。     

Burton–Miller-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.     

A study on avoiding hyper-singular integrals in exterior acoustic radiation analysis

Burton and Miller's formulation (BMF) is well known for solving non-uniqueness problems in an exterior acoustic problem. In this paper, an IBMF (improved BMF) is presented by directly taking into account of the normal derivatives of solid angles on the surface. The introduction of the normal derivative of solid angles in building the Burton and Miller's formulation will cause the elimination of hyper-singular terms. Here lies the advantage of improved formulation. Then, reproduction of diagonal terms (RDT) is applied to further improve singular integral in IBMF. Three cases are compared for closed axi-symmetric surfaces, including SHF (conventional surface Helmholtz formulation), IBMF and IBMF-RDT. Numerical results show that the combined scheme of IBMF-RDT is valid at all wave numbers, and gives satisfied calculating accuracy, which is about 10 times higher than that of SHF.     

A finite-element model of the aperture method for determining the effective radiating area of physiotherapy treatment heads

This paper describes a theoretical study of the way in which a circular aperture placed in front of a plane-piston modifies the ultrasonic field it generates. Specifically, the radiated acoustic power transmitted by the aperture and the radiation force experienced by an absorbing target placed in the transmitted beam, are evaluated at a distance from the exit-side of the aperture. The calculations used a finite element (FE) method, in conjunction with a surface Helmholtz integral formulation to solve the fluid/structure interaction problem. The PAFEC (Program for Automatic Finite Element Computation) vibroacoustics software was used for the FE calculations and the implementation of the surface Helmholtz integral formulation method. Acoustic pressures and particle velocities were computed at required points, whilst accounting for the presence of the aperture in the medium, together with its dynamic properties when subjected to an incident sound field. This enabled the calculation of the radiation force on the absorber and its variation with the applied aperture diameter was investigated. As part of the validation process for the new FE aperture model, the ratio of radiation force to radiated acoustic power obtained using the FE method in the unapertured case, through the use of the Rayleigh integral, yielded good agreement with results obtained through an analytical solution. The study has been carried out to provide a better understanding of the factors affecting the measurement uncertainty for the aperture method of determining the effective radiating area (A(ER)) of physiotherapy ultrasound treatment heads.     

Fast multipole accelerated boundary element method for the Helmholtz equation in acoustic scattering problems

We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficultie...     

A new simulation method for the determination of forces in polymer/colloid systems

We introduce a new simulation method, which we call the contact-distribution method, for the determination of the Helmholtz potential for polymer/colloid systems from lattice Monte-Carlo simulations. This method allows one to obtain forces between finite or semi-infinite objects of any arbitrary shape and dimensions in the presence of polymer chains in solution or physisorbed or chemisorbed at interfaces. We illustrate the application of the method using two examples: (i) the interaction between the tip of an atomic force microscope (AFM) and a single, end-grafted polymer chain and (ii) the interaction between an AFM tip and a polymer brush. Numerical results for the first two cases illustrate how the method can be used to confirm and extend scaling laws for forces and Helmholtz potentials, to examine the effects of the shapes and sizes of the objects and to examine conformational transitions in the polymer chains. Received: 15 May 1998 / Revised: 11 June 1998 / Accepted: 12 June 1998     

Approximation of the parabolic equation and the wavefield of a point source in a layered random medium

This paper studies the wavefield of a source in a multidimensional randomly layered medium. They obtained asymptotical expressions of the wave statistical characteristics for different boundary conditions both in the framework of the parabolic equation approximation and the exact formulation of the boundary problem for the Helmholtz equation. It is shown that the presence of a small but finite absorption γ is most important for the statistics. The diffraction effects turn out to be like those of absorption, but γ cannot tend to zero in this problem. In an appendix they give the factorization formulae of the wave equation solution in a layered medium.