The clinical repair of teeth using direct filling materials: engineering considerations

This paper reviews the way in which teeth damaged by caries may be repaired clinically. The mechanical effects of caries are described, as are the materials available to repair the damage caused by this disease. Studies are reported which have shown that caries reduces the compressive strength of the tooth to less than 50 per cent of its original value and that, by use of appropriate materials and placement techniques, this can be restored to some 80 per cent of this value. However, very few studies have been carried out which view tooth repair from an engineering perspective. Instead, emphasis is placed on determining clinical durability of repairs. This is related to repair strength but brings in other factors, such as the oral hygiene of the patient. Despite this complication, durability studies show that modern restorative materials perform well under clinical conditions, from which it may be concluded that the repair process allows a structure to be fabricated that is essentially sound from an engineering viewpoint, even if inferior to the original tooth structure provided by nature.     

On the Efficient Evaluation of Coalescence Integrals in Population Balance Models

The solution of population balance equations is a function f(t,r,x) describing a population density of particles of the property x at time t and space r. For instance, the additional independent variable x may denote the particle size. The describing partial differential equation contains additional sink and source terms involving integral operators. Since the coordinate x adds at least one further dimension to the spatial directions and time coordinate, an efficient numerical treatment of the integral terms is crucial. One of the more involved integral terms appearing in population balance models is the coalescence integral, which is of the form ∫ ₀ ^x κ(x–y, y) f(y) f(x–y)dy. In this paper, we describe an evaluation method of this integral which needs only operations, where n is the number of degrees of freedom with respect to the variable x. This cost can also be obtained in the case of a grid geometrically refined towards x=0.     

Optimal neuronal tuning for finite stimulus spaces

The efficiency of neuronal encoding in sensory and motor systems has been proposed as a first principle governing response properties within the central nervous system. We present a continuation of a theoretical study presented by Zhang and Sejnowski, where the influence of neuronal tuning properties on encoding accuracy is analyzed using information theory. When a finite stimulus space is considered, we show that the encoding accuracy improves with narrow tuning for one- and two-dimensional stimuli. For three dimensions and higher, there is an optimal tuning width.     

Motif analysis for automatic segmentation of CT surface contours into individual surface features

This paper presents a new layer-based technique for automatic high-level segmentation of 3-D surface contours into individual surface features through motif analysis. The procedure starts from a contour-based surface model representing a composite surface area of an object. For each of the surface contours, a relative turning angle (RTA) map is derived. The RTA map usually contains noise and minor features. Algorithms based on motif analysis are applied for extracting a main profile of the RTA map free from background noise and other minor features. All feature points on the extracted profile are further identified from the extracted main profile through further motif analysis. The original contour is thus partitioned into individual segments with the identified feature points. A collection of consecutive contour segments among different layers form an individual 3-D surface feature of the original composite surface. The developed approach using motif analysis is particularly useful for the identification of smooth joins between individual surface features and for the elimination of superposed noise and unwanted minor features.     

A simple strategy for calibrating the geometry of light sources

We present a methodology for calibrating multiple light source locations in 3D from images. The procedure involves the use of a novel calibration object that consists of three spheres at known relative positions. The process uses intensity images to find the positions of the light sources. We conducted experiments to locate light sources in 51 different positions in a laboratory setting. Our data shows that the vector from a point in the scene to a light source can be measured to within 2.7±4° at α=.05 (6 percent relative) of its true direction and within 0.13±.02 m at α=.05 (9 percent relative) of its true magnitude compared to empirically measured ground truth. Finally, we demonstrate how light source information is used for color correction     

Parallel Implementation of a Class of Adaptive Signal Processing Applications

Recently, High Performance Computing (HPC) platforms have been employed to realize many computationally demanding applications in signal and image processing. These applications require real-time performance constraints to be met. These constraints include latency as well as throughput. In order to meet these performance requirements, efficient parallel algorithms are needed. These algorithms must be engineered to exploit the computational characteristics of such applications. In this paper we present a methodology for mapping a class of adaptive signal processing applications onto HPC platforms such that the throughput performance is optimized. We first define a new task model using the salient computational characteristics of a class of adaptive signal processing applications. Based on this task model, we propose a new execution model. In the earlier linear pipelined execution model, the task mapping choices were restricted. The new model permits flexible task mapping choices, leading to improved throughput performance compared with the previous model. Using the new model, a three-step task mapping methodology is developed. It consists of (1) a data remapping step, (2) a coarse resource allocation step, and (3) a fine performance tuning step. The methodology is demonstrated by designing parallel algorithms for modern radar and sonar signal processing applications. These are implemented on IBM SP2 and Cray T3E, state-of-the-art HPC platforms, to show the effectiveness of our approach. Experimental results show significant performance improvement over those obtained by previous approaches. Our code is written using C and the Message Passing Interface (MPI). Thus, it is portable across various HPC platforms. Received April 8, 1998; revised February 2, 1999.     

Acceleration and filtering in the generalized Landweber iteration using a variable shaping matrix

The generalized Landweber iteration with a variable shaping matrix is used to solve the large linear system of equations arising in the image reconstruction problem of emission tomography. The method is based on the property that once a spatial frequency image component is almost recovered within in in the generalized Landweber iteration, this component will still stay within in during subsequent iterations with a different shaping matrix, as long as this shaping matrix satisfies the convergence criterion for the component. Two different shaping matrices are used: the first recovers low-frequency image components; and the second may be used either to accelerate the reconstruction of high-frequency image components, or to attenuate these components to filter the image. The variable shaping matrix gives results similar to truncated inverse filtering, but requires much less computation and memory, since it does not rely on the singular value decomposition.     

Resonance Raman and optical transient studies on the light-induced proton pump of bacteriorhodopsin reveal parallel photocycles

The photocycle of bacteriorhodopsin (bR) was studied at ambient temperature in aqueous suspensions of purple membranes using time-resolved resonance Raman (RR) and optical transient spectroscopy (OTS). The samples were photolyzed, and the fractional concentrations of the retinylidene chromophore in its parent state, BR570, and in the intermediate states L550, M412, N560, and O640 were determined in the time domain 20 microseconds-1 s and in the pH range 4-10.5. Two kinetically different L components could be identified. At pH 7 one fraction of L (approximately 65%) decays in 80 microseconds to M (deprotonation of the Schiff base), whereas the residual part is converted in approximately 0.5 ms to N. The RR spectra reveal only minor structural changes of the chromophore in the L-->N transition. These were attributed to a conformational change of the protein backbone [Ormos, P., Chu, K., & Mourant, J. (1992) Biochemistry 31, 6933]. With decreasing pH the L-->N transition is delayed to > 2 ms following a titration-like function with pKa approximately 6.2. The decay of M412 monitored by OTS can be fitted for each pH value by two different amplitudes and time constants (Mf, tau f; Ms, tau s; f = fast, s = slow). Both Mf and Ms consist of subcomponents which can be distinguished by their different reaction pathways (but not by OTS). Mf occurs in the reaction sequences L-->Mf-->N-->BR and L-->Mf-->O-->BR. The population of the first sequence, in which N is formed with the time constant tau f (approximately 2-4 ms, pH 6-10.5), increases with pH. Ms is also found in two different reaction sequences of the form L-->Ms-->BR. The quantitative analysis reveals that each "titration effect" can be related to a certain fraction of bR. It is proposed that each fraction can be identified with a "subspecies" of bR which undergoes an independent and individual cyclic reaction. A complete reaction scheme is set up which represents the manifold of observed phenomena. It is concluded from the pH dependence of the lifetimes of Ms and N that the reconstitution of BR570 in the reaction steps Ms-->BR and N-->BR requires the uptake of a proton from the external phase. It is argued that this proton catalyzes the reisomerization of retinal, whereas the Schiff base is internally reprotonated from Asp-85. A model for proton pumping is proposed in which the proton taken up from the external phase to catalyze the reisomerization of retinal is the one which is pumped through the membrane during the photocycle of bR.     

Fluid dynamics of a quantized vortex filament in a hole

The behavior of nucleated vortex loops and of remanent vortex filaments in idealized circular and slit-like orifices has been investigated by direct computation. It is found that such vortices can be stretched by the diverging flow on the exit side of the orifice. The energy needed to stretch the vortex is abstracted from the flow field and observed as dissipation. This occurs in the form of discrete phase-slip events in the case of nucleated vortex loops, and in the form of multiple-phase-slip cascades when a remanent vortex is involved.