A modified Lax-Wendroff correction for wave propagation in media described by Zener elements

A modified Lax-Wendroff correction for wave propagation in attenuating and dispersive media described by Zener elements is presented. As opposed to the full correction, this new technique is explicit and offers large computational savings. The technique may be applied to a wide variety of hyperbolic problems. Here, the concept is illustrated for wave propagation in visco-acoustic medi

Correcting an acoustic wavefield for elastic effects

Finite-difference simulations are an important tool for studying elastic and acoustic wave propagation, but remain computationally challenging for elastic waves in three dimensions. Computations for acoustic waves are significantly simpler as they require less memory and operations per grid cell, and more significantly can be performed with coarser grids, both in space and time. In this paper, we present a procedure for correcting acoustic simulations for some of the effects of elasticity, at a cost considerably less than full elastic simulations. Two models are considered: the full elastic model and an equivalent acoustic model with the same P velocity and density. In this paper, although the basic theory is presented for anisotropic elasticity, the specific examples are for an isotropic model. The simulations are performed using the finite-difference method, but the basic method could be applied to other numerical techniques. A simulation in the acoustic model is performed and treated as an approximate solution of the wave propagation in the elastic model. As the acoustic solution is known, the error to the elastic wave equations can be calculated. If extra sources equal to this error were introduced into the elastic model, then the acoustic solution would be an exact solution of the elastic wave equations. Instead, the negative of these sources is introduced into a second acoustic simulation that is used to correct the first acoustic simulation. The corrected acoustic simulation contains some of the effects of elasticity without the full cost of an elastic simulation. It does not contain any shear waves, but amplitudes of reflected P waves are approximately corrected. We expect the corrected acoustic solution to be useful in regions of space and time around a P-wave source, but to deteriorate in some regions, for example, wider angles, and later in time, or after propagation through many interfaces. In this paper, we outline the theory of the correction method, and present results for simulations in a 2-D model with a plane interface. Reflections from a plane interface are simple enough that an analytic analysis is possible, and for plane waves, we give the correction to the acoustic reflection and transmission coefficients. Finally, finite-difference calculations for plane waves are used to confirm the analytic results. Results for wave propagation in more complicated, realistic models will be presented elsewher

The Canadian Joint Replacement Registry—what have we learned?

The Canadian Joint Replacement Registry (CJRR) was launched in 2000 through the collaborative efforts of the Canadian Orthopedic Association and the Canadian Institutes for Health Information. Participation is voluntary, and data collected by participating surgeons in the operating room is linked to hospital stay information from administrative databases to compile yearly reports. In the fiscal year 2006–2007, there were 62,196 hospitalizations for hip and knee replacements in Canada, excluding Quebec. This represents a 10-year increase of 101% and a 1-year increase of 6%. Compared to men, Canadian women have higher age-adjusted rates per 105 for both TKA (148 vs. 110) and THA (86 vs. 76). There also exist substantial inter-provincial variations in both age-adjusted rates of arthroplasty and implant utilization that cannot be explained entirely on the basis of differing patient demographics. The reasons for these variations are unclear, but probably represent such factors as differences in provincial health expenditure, efforts to reduce waiting lists, and surgeon preference. The main challenge currently facing the CJRR is to increase procedure capture to > 90%. This is being pursued through a combination of efforts including simplification of the consent process, streamlining of the data collection form, and the production of customized reports with information that has direct clinical relevance for surgeons and administrators. As the CJRR continues to mature, we are optimistic that it will provide clinically important information on the wide range of factors that affect arthroplasty outcome

Free and smooth boundaries in 2-D finite-difference schemes for transient elastic waves

A method is proposed for accurately describing arbitrary-shaped free boundaries in single-grid finite-difference schemes for elastodynamics, in a time-domain velocity-stress framework. The basic idea is as follows: fictitious values of the solution are built in vacuum, and injected into the numerical integration scheme near boundaries. The most original feature of this method is the way in which these fictitious values are calculated. They are based on boundary conditions and compatibility conditions satisfied by the successive spatial derivatives of the solution, up to a given order that depends on the spatial accuracy of the integration scheme adopted. Since the work is mostly done during the preprocessing step, the extra computational cost is negligible. Stress-free conditions can be designed at any arbitrary order without any numerical instability, as numerically checked. Using 10 grid nodes per minimal S-wavelength with a propagation distance of 50 wavelengths yields highly accurate results. With 5 grid nodes per minimal S-wavelength, the solution is less accurate but still acceptable. A subcell resolution of the boundary inside the Cartesian meshing is obtained, and the spurious diffractions induced by staircase descriptions of boundaries are avoided. Contrary to what occurs with the vacuum method, the quality of the numerical solution obtained with this method is almost independent of the angle between the free boundary and the Cartesian meshing.Comment: accepted and to be published in Geophys. J. In

Extending an industrial root controller : implementation and applications of a fast open sensor interface

An overview is given of the design and implementation of a platform for fast external sensor integration in an industrial robot system called ABB S4CPlus. As an application and motivating example, the implementation of force-controlled grinding and deburring within the AUTOFETT-project is discussed. Experiences from industrial usage of the fully developed prototype confirms the appropriateness of the design choices, thus also confirming the fact that control and software need to be tightly integrated. The new sensor can be used for the prototyping and development of a wide variety of new application