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.     

Heavy Majorana neutrinos in e + e − → W + W−

In this paper we investigate the effects of heavy Majorana neutrinos in the reaction e ⁺ e ^? → W ⁺ W ^?. We consider neutrino masses in the 1–10 TeV region. We show that at LEP II and NLC energies it is possible to use this processes to verify indirect evidence of heavy neutral particles with mixing angles of the order sin² α = 0.01. We discuss the unitarity restrictions that can be obtained for vector singlet and fermion-mirror-fermion models.     

Vibrational spectrum of uridine obtained by means of nonresonant-background-free CARS spectroscopy

Received: 18 June 1996/Revised version: 3 January 1997     

Double-diffusive convection for a heated cylinder submerged in a salt-stratified fluid layer

Double-diffusive convection due to a cylindrical source submerged in a salt-stratified solution is numerically investigated in this study. For proper simulation of the vortex generated around the cylinder, a computational domain with irregular shape is employed. Flow conditions depend strongly on the thermal Rayleigh number, Ra _T , and the buoyancy ratio, R _ρ. There are two types of onset of instability existing in the flow field. Both types are due to either the interaction of the upward temperature gradient and downward salinity gradient or the interaction of the lateral temperature gradient and downward salinity gradient. The onset of layer instability due to plume convection is due to the former, whereas, the onset of layer instability of layers around the cylinder is due to the latter. Both types can be found in the flow field. The transport mechanism of layers at the top of the basic plume belongs to former while that due to basic plume and layer around the cylinder are the latter. The increase in Ra _T reinforces the plume convection and reduces the layer numbers generated around the cylinder for the same buoyancy ratio. For the same Ra _T , the increase of R _ρ suppresses the plume convection but reinforces the layers generated around the cylinder. The profiles of local Nusselt number reflects the heat transfer characteristics of plume convection and layered structure. The profiles of averaged Nusselt number are between the pure conduction and natural convection modes and the variation is due to the evolution of layers. Received on 13 September 1996     

Time-harmonic response of a vertical crack in plates

A technique for the characterization of an interior crack perpendicular to the surface of the plates is presented. A time-harmonic line load is applied on the upper surface of the plate. The scattered displacement field is calculated by using the strip element method. It is found that there are additional superimposed oscillations in the curve of the absolute values of the displacement in the region between the load and the crack. The mean value of the peaks of the additional oscillations increases with the crack length, and is less dependent on the distance from the load to the crack. The crack position is determined from the displacement distribution curve and the crack length is determined from the mean value of the peaks of the additional oscillations. Approximate polynomial formula of degree four for predicting the length of interior vertical cracks in isotropic homogeneous plates is also proposed. Numerical examples are given to demonstrate the present technique.     

Magnetic resonance imaging of rigid polymers at elevated temperatures with SPRITE

Magnetic resonance imaging has rarely been applied to rigid polymeric materials, due primarily to the strong dipolar coupling and short signal lifetimes inherent in these materials. SPRITE (single point ramped imaging withT ₁ enhancement) (B. J. Balcom, R. P. MacGregor, S. D. Beyea, D. P. Green, R. L. Armstrong, T. W. Bremner: J. Magn. Reson. A123, 131–134, 1996) is particularly well suited to imaging solid materials. With SPRITE, the only requirement is thatT ₂^* be long enough so that the signal can be phase-encoded. The minimum phase encoding time is limited by the maximum gradient strength available and by the instrument deadtime. At present this is usually tens of microseconds and will only improve with refinements in technology. We have used the SPRITE sequence in conjunction with raising the sample temperature to obtain images of rigid polymers that have largely frustrated conventional imaging methods. This approach provides a straightforward and reliable method for imaging a class of samples that, up until now, have been very difficult to image.     

Development of high-speed video cameras for Dynamic PIV

The most promising next generation Image Velocimtry (IV) is the high-speed Dynamic PIV. It requires the development of innovative high-speed video camera sensors. We started by specifying the required performance of these new sensors, for measurements in air and water flows. These criteria founded on the most recent developments in PIV algorithms and incorporate results from a large questionnaire survey of users of high-speed video cameras in Japan. The results suggest that the followings are required: (1) frame rate of 1,000,000 fps, (2) pixel count of 1,000,000 pixels, (3) frame storage capacity of 100–200 frames for tracing a single event and 10,000 frames for turbulent measurements, (4) gray levels of 4–8 bits for PTV; 12 bits for observation. Finally, we reviewed the state of the art of high-speed video-image sensors. Currently the standard parallel-readout sensors can operate at 1 Kfps with a pixel count of approximately 1 Kpixels. The In-situ Storage Image Sensor (ISIS) developed by the authors has recently achieved frame rates of 1 Mfps.     

The 2D Rancieite-type manganic acid and its alkali-exchanged derivatives: Part I — Chemical characterization and thermal behavior

The 2D Rancieite type manganic acid was prepared by reduction of KMnO₄ in acidic medium. Its ion exchange behavior allows to prepare alkali derivatives. All compounds were characterized with use of a combination of X-ray diffraction, chemical analyses, TGA, magnetic measurements and spectroscopic techniques. The evolution of their chemical composition versus temperature was studied between 180 and 400 °C. It shows that the dehydration process is partly reversible in these compounds whereas the weak reduction is irreversible. The 2D Rancieite-type manganic acid is readily different from a Birnessite-type phyllomanganate, as shown by several features: the interlayer distance, the ion exchange capacity, the thermal behavior, the interlayer cation content, the manganese average oxidation state, the magnetic behavior and the IR spectrum.     

Application of stochastic optimal control to the suboptimal design of thrust-mass profiles

Stochastic optimal control techniques are applied to compare the performance of identical medium-range air-to-air missiles which have different thrust-mass profiles. The measure of the performance is the probability of reaching a lock-on-point with a favorable range of guidance and flight parameters, during a fixed time interval [0,t _f ], given that, during the flight, the trajectories of the missile are subjected to a variety of constraints including dynamic pressure constraints.