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 evolution of a truncated gaussian probability density through time—modelling animal liveweights after selection

The form of the probability density derived from the evolution in time of a previously truncated frequency distribution of animal Liveweights is of interest in animal husbandry. Truncated frequency distributions arise when the heavier animals are sold for slaughter and the lighter animals retained. The demands of modern quality assurance schemes require that, given information on animal growth, the farmer is able to estimate the number of animals that would meet the specifications at some time in the future after truncation. Assuming that animal growth can be described by a linear stochastic differential equation, we derive an explicit expression for the probability density of animal Liveweights at any time after the truncation of an initial Gaussian density. It is shown that this probability density converges rapidly to a Gaussian density, so that after about 20 days of typical growth rates for lambs, the resulting density is practically indistinguishable from Gaussian.     

The Best Lower Bound for the Permanent of a Correlation Matrix of Rank Two

If 1≤k≤n, then Cor(n,k) denotes the set of all n×n real correlation matrices of rank not exceeding k. Grone and Pierce have shown that if A∈Cor (n, n-1), then per(A)≥n/(n-1). We show that if A∈Cor(n,2), then , and that this inequality is the best possible.     

Three-dimensional structure of Saccharomyces cerevisiae inorganic pyrophosphatase complexed with cobalt and phosphate ions

Crystals of Saccharomyces cerevisiae inorganic pyrophosphatase suitable for X-ray diffraction study were grown by cocrystallization of the enzyme with cobalt chloride and imidodiphosphate. Saccharomyces cerevisiae is a metal-dependent enzyme which catalyzes hydrolysis of inorganic pyrophosphate to orthophosphate. The three-dimensional structure of this enzyme was solved by the molecular-replacement method and refined at 1.8 Å resolution to an R factor of 19.5%. Cobalt and phosphate ions were revealed in the active centers of both identical subunits (A and B) of the pyrophosphatase molecule. In subunit B, a water molecule was found between two cobalt ions. It is believed that this water molecule acts as an attacking nucleophile in the enzymatic cleavage of the pyrophosphate bond. It was demonstrated that cobalt ions and a phosphate group occupy only part of the potential binding sites (two chemically identical and crystallographically independent subunits have different binding sites). The arrangement of ligands and the structure of the nucleophile-binding site are discussed in relation to the mechanism of action of the enzyme and the nature of the metal activator.     

Retardation characteristics and birefringence of a multiple-order crystalline quartz plate

A crystalline quartz plate of thickness 1.070 mm is calibrated between 370 and 794 nm. Throughout this spectral interval, the retardance varied by 32π and the plate introduced quarterwave retardance 16 times at different wavelengths. The birefringence (n_e−n_o) of crystalline quartz was calculated as a single quantity and varied from 0.00971 at 370 nm to 0.00891 at 794 nm. All measurements were carried out at 23°C.     

Performance of a CW pumped Nd:YVO4 amplifier with kHz pulses

Until recently, simple and reliable high repetition-rate laser sources with nanosecond pulses much shorter than from conventional A-O Q-switch lasers were not available. However over the past 2 years we have developed such lasers based on proprietary fast E-O switching technology, which allows designs delivering 1 ns pulses and subnanosecond jitter for good synchronisation. The technology provides pulses with multi-kW peak power and repetition-rates to >100 kHz.Most recently, the performance of these short pulse lasers has been developed further by implementing oscillator/amplifier (master oscillator and power amplifier, MOPA) technology which increases the output to >1 W average power. Here we report on a simple model that has been used to predict the performance of the CW pumped Nd:YVO₄ amplifier used in the MOPA laser. The model is based on the well-known expressions for the saturated gain applying to laser pulses, but more usually applied to pulse-excited amplifiers. The model is shown to allow a good interpretation of the amplifier behaviour for kHz pulses and to be a useful tool for predicting the performance of the MOPA laser.     

Hutchinson's shear coefficient for anisotropic beams

Timoshenko's theory of vibrating beams requires a shear correction factor to correctly take into account the effects of shear deformation for different beam cross-sections. This correction is crucial for a precise determination of the shear modulus from the resonant frequencies. Hutchinson's beam theory is used to derive a new shear correction coefficient for anisotropic materials. A comparison is made with other shear coefficients for anisotropic materials published in the literature. Computer-simulated spectra are used to validate the new anisotropic shear correction coefficient.     

Refraction holodiagrams and Snell's law

The refraction holodiagram RHD is analyzed here with respect to the law of refraction. Particularly, we study the surface that exactly conjugates by refraction a virtual point source with a real image or conversely. By using the total optical path as a parameter we build a diagram that consists in a family of Descartes ovals of the apple type that contains the Pascal's limaçon as a particular extreme case and the spherical surface with the Weierstrass points as another. These representations permit the straightforward application of Fermat's principle in the case of arbitrary refracting surfaces and show the shape of generalized Fresnel's zones in the intersections with any surface. Snell's law is applied to rays incident on the apple type surfaces to find out the conditions for exact conjugation. Sensitivity to optical path variations is also discussed. The RHD curves family can be represented in a Cartesian way where the ovals appear as equally spaced straight lines.     

Determination of flutter derivatives via a neural network approach

A neural-network-based method is offered to determine the flutter derivatives of section models under smooth and turbulent flows. The approach uses the observed dynamic responses to train an appropriate neural network. Subsequently, the modal parameters of the model for different mean velocities of wind flow are directly estimated using weight matrices in the neural network. The flutter derivatives can then be determined accurately. The validity of the present method is verified through numerical studies. Finally, the procedure is employed to process experimental data from an inverted-U-type section model, obtained from wind tunnel tests.     

两个回归函数相等的Bootstrap检验

两个回归参数相等性检验一直是统计界感兴趣的问题之一．在这篇文章中，四个检验统计量被用于度量两曲线的差异，在原假设下统计量的分布采用向量数据的重复抽样来逼近，并给出了—些模拟结果．