Fast finite difference Poisson solvers on heterogeneous architectures

In this paper we propose and evaluate a set of new strategies for the solution of three dimensional separable elliptic problems on CPU–GPU platforms. The numerical solution of the system of linear equations arising when discretizing those operators often represents the most time consuming part of larger simulation codes tackling a variety of physical situations. Incompressible fluid flows, electromagnetic problems, heat transfer and solid mechanic simulations are just a few examples of application areas that require efficient solution strategies for this class of problems. GPU computing has emerged as an attractive alternative to conventional CPUs for many scientific applications. High speedups over CPU implementations have been reported and this trend is expected to continue in the future with improved programming support and tighter CPU–GPU integration. These speedups by no means imply that CPU performance is no longer critical. The conventional CPU-control–GPU-compute pattern used in many applications wastes much of CPU’s computational power. Our proposed parallel implementation of a classical cyclic reduction algorithm to tackle the large linear systems arising from the discretized form of the elliptic problem at hand, schedules computing on both the GPU and the CPUs in a cooperative way. The experimental result demonstrates the effectiveness of this approach.     

Effects of physical refining on contents of waxes and fatty alcohols of refined olive oil

Changes in the contents of waxes and fatty alcohols during deodorization/physical refining of bleached olive oil were studied. Experiments were carried out with 1.85% acidity oil, which was physically refined in a discontinuous deodorizer of 250-kg maximum capacity using nitrogen as stripping gas instead of steam. The variables studied were load and temperature of oil in the deodorizer as well as N₂ flow. Analyses of waxes and alcohols were carried out at different operation times. The maximum content of wax was always observed when the oil reached the deodorization temperature. The variation in the wax content depended on temperature and N₂ flow. Wax decomposition started and continued during the operating time, and a progressive decrease, which was pronounced between 3 and 4 h, was observed. Small changes in waxes were observed between 4 and 5 h. Total content of fatty alcohols diminished throughout the operating time, and changes did not depend on the variables studied.     

Single frame correction of motion artifacts in PMD-based time of flight cameras

One of the leading time of flight imaging technologies for depth sensing is based on Photonic Mixer Devices (PMD). In PMD sensors each pixel samples the correlation between emitted and received light signals. Current PMD cameras compute eight correlation samples per pixel in four sequential stages to obtain depth with invariance to signal amplitude and offset variations. With motion, PMD pixels capture different depths at each stage. As a result, correlation samples are not coherent with a single depth, producing artifacts. We propose to detect and remove motion artifacts from a single frame taken by a PMD camera. The algorithm we propose is very fast, simple and can be easily included in camera hardware. We recover depth of each pixel by exploiting consistency of the correlation samples and local neighbors of the pixel. In addition, our method obtains the motion flow of occluding contours in the image from a single frame. The system has been validated in real scenes using a commercial low-cost PMD camera and high speed dynamics. In all cases our method produces accurate results and it highly reduces motion artifacts.     

A Class of Incomplete Riemann Solvers Based on Uniform Rational Approximations to the Absolute Value Function

In this paper we propose a new class of incomplete Riemann solvers, based on approximations in the \(L^\infty \) -norm to the absolute value function in \([-1,1]\) by means of rational functions, for the numerical approximation of the solution of hyperbolic systems of conservation laws. The main idea relies on the construction of a numerical approximation to the viscosity matrix by using an appropriate rational real function \(R(x)\) , that approximates the function \(|x|\) uniformly in \([-1,1]\) , evaluated at the Jacobian of the fluxes of the hyperbolic system computed at some average value (for example, Roe averages). In addition to the Jacobians of the fluxes we shall use either the maximum in absolute value of the characteristic speeds in each cell or an upper bound of them. Thus, the resulting approximate Riemann solver is incomplete in the sense that we do not use the complete spectral decomposition of the Jacobian. Moreover, the new class of Riemann solvers consists of a hierarchy of schemes running from the more dissipative to the less dissipative ones, and having as limiting case a Roe-like scheme. According to the order of the approximation of the generating rational function used, the degree of dissipation can be dosed for particular applications. We study different rational approximations: Newman-type functions, iterative generated Halley functions, and also Chebyshev polynomial approximants. We test our basic algorithms for different initial value Riemann problems for ideal gas dynamics (HD) and magnetohydrodynamics (MHD) to observe their behavior with respect to challenging scenarios in numerical simulations, including some standard numerical pathologies (e. g., heat conduction, postshock oscillations and overheating) and the formation of compound waves in ideal MHD. We also examine our proposed schemes, by computing the numerical approximation of different initial value problems for nonconservative multilayer shallow water equations, where it has been observed that intermediate waves can be properly captured for an appropriate degree of approximation of the generating rational function used. Our numerical tests indicate that the proposed schemes are robust, running stable and accurate with a satisfactory time step restriction (CFL constant), and the computational cost is more advantageous with respect to schemes that use a complete spectral decomposition of the Jacobians.     

Degradation of Leaf Litter Phenolics by Aquatic and Terrestrial Isopods

To investigate species-specific decomposition rates of litter from native (Quercus faginea) and introduced (Eucalyptus globulus) tree species in Portugal, we monitored changes in the phenolic signature of leaf litter during decomposition as mediated by an aquatic, Proasellus coxalis (Isopoda: Asellota), and two terrestrial, Porcellio dispar and Eluma caelatum (Isopoda: Oniscidea), detritivores. Although the litter of Eucalyptus and Quercus did not differ in overall protein precipitation capacity, we detected differences in terms of contents of particular phenolic compounds and phenol oxidation products. Accordingly, we observed food-specific consumption rates in Proasellus, but not in the terrestrial isopods. Proasellus digested Eucalyptus at significantly higher rates than Quercus, whereas the opposite was the case for Eluma, and Porcellio digested both litter types equally well. Despite slight differences in detail, effects of Proasellus on changes in the signature of litter phenolics were similar for both litter types, whereas terrestrial isopods—Porcellio and Eluma, although they differed from each other—digestively degraded phenolic compounds in Eucalyptus and Quercus litter, respectively, in different ways. Overall, however, degradation of litter phenolics was similarly effective on both litter types. From these data, we conclude that decomposition of Eucalyptus litter does not proceed more slowly than of litter from native Portuguese trees.     

Preparation and properties of high molecular weight poly(1-germa-) or poly(1-sila-cis-pent-3-enes)

High molecular weight poly(1,1-dimethyl-1-germa-cis-pent-3-ene), poly(1,1-diphenyl-1-germa-cis-pent-3-ene), poly(1,1-dimethyl-1-sila-cis-pent-3-ene), and poly(1-methyl-1-phenyl-1-sila-cis-pent-3-ene) have been prepared. The thermal stability of these polymers is found to increase with their molecular weight.     

Effect of refining on the phenolic composition of crude olive oils

By definition, virgin olive oil is consumed unrefined, although a great proportion of the olive oil produced has to be refined to render it edible. Phenolic compounds are among the substances eliminated during the refining process; in the present work these were characterized by HPLC, and their evolution during the different refining steps was studied. The complete refining process removed most polyphenols from oils, but the behavior of individual compounds at each step also was observed. o-Diphenols (hydroxytyrosol, catechol, and hydroxytyrosol acetate) and flavonoids (luteolin and apigenin) were eliminated first during the alkaline treatment. Tyrosol and 4-ethylphenol remained in the oil until the deodorization step. A large amount of phenolic compounds was discovered in the refining by-products such as soapstocks and deodorization distillates. In the latter streams, the concentrations of tyrosol and 4-ethylphenol reached up to 149 and 3720 mg/kg by-product, respectively. This high level of 4-ethylphenol and its well-known strong off-odor can interfere during further processing of the deodorization distillates, and this must be taken into account when deciding what is to become of them. Similarly, the results of this work open the possibility of recovering phenolic compounds from the “second centrifugation olive oils” by adding a new washing step prior to the refining process. By including this new step, the most polar polyphenols, hydroxytyrosol and tyrosol, will diffuse from oil to water and a concentration of up to 1400 mg/L of hydroxytyrosol may be achieved.     

Ambient and High-Temperature Stable Fracture Tests in Ceramics: Applications to Yttria-Partially-Stabilized Zirconia

A new technique that provides stable fracture tests in brittle ceramic materials at room and at high temperature is presented. This technique uses the crack mouth opening displacement signal from a laser extensometer to feed the servo control of a hydraulic testing machine. In this way it is possible to obtain—in addition to the fracture toughness—the fracture energy and the R -curve, in a single test.