A comparison of in situ constant and falling head permeameter tests to assess the distribution of clogging within horizontal subsurface flow constructed wetlands

Clogging is the main operational problem associated with horizontal subsurface flow constructed wetlands (HSSF CWs). The measurement of saturated hydraulic conductivity has proven to be a suitable technique to assess clogging within HSSF CWs. The vertical and horizontal distribution of hydraulic conductivity was assessed in two full-scale HSSF CWs by using two different in situ permeameter methods (falling head (FH) and constant head (CH) methods). Horizontal hydraulic conductivity profiles showed that both methods are correlated by a power function (FH= CH 0.7821, r 2=0.76) within the recorded range of hydraulic conductivities (0-70 m/day). However, the FH method provided lower values of hydraulic conductivity than the CH method (one to three times lower). Despite discrepancies between the magnitudes of reported readings, the relative distribution of clogging obtained via both methods was similar. Therefore, both methods are useful when exploring the general distribution of clogging and, specially, the assessment of clogged areas originated from preferential flow paths within full-scale HSSF CWs. Discrepancy between methods (either in magnitude and pattern) aroused from the vertical hydraulic conductivity profiles under highly clogged conditions. It is believed this can be attributed to procedural differences between the methods, such as the method of permeameter insertion (twisting versus hammering). Results from both methods suggest that clogging develops along the shortest distance between water input and output. Results also evidence that the design and maintenance of inlet distributors and outlet collectors appear to have a great influence on the pattern of clogging, and hence the asset lifetime of HSSF CWs. © Springer Science+Business Media B.V. 2011.

Powers and Logarithms

There are applied power mappings in algebras with logarithms induced by a given linear operator D in order to study particular properties of powers of logarithms. Main results of this paper will be concerned with the case when an algebra under consideration is commutative and has a unit and the operator D satisfies the Leibniz condition, i.e. D(xy) = xDy + yDx for x, y ∈ dom D. Note that in the Number Theory there are well-known several formulae expressed by means of some combinations of powers of logarithmic and antilogarithmic mappings or powers of logarithms and antilogarithms (cf. for instance, the survey of Schinzel S[1].

The Direct and Inverse Spectral Problems for some Banded Matrices

2000 Mathematics Subject Classification: 15A29.

A New Method for Estimating the Moisture Content and Flexibility of Polymerised Bentonite Clay Mat

This paper presents a scientific development to address the current absence of a convenient technique to identify the ductile to brittle transition of bentonite clay mats. The instrumented indentation and 3-point bending tests were performed on different liquid polymer hydrated bentonite clay mats at varying moisture content. Properties measured include modified Brinell Hardness Number (BHN) and elastic structural stiffness (EI). The dependence of flexural stiffness on moisture content is demonstrated to conform to a best power function variation. The ductile to brittle transition of clay mat is affected primarily by the change in the moisture content and for the clay mat to remain flexible, critical moisture content of 1.7 times of its plastic limit is required. Results also indicate that a strong correlation between indentation hardness and the structural stiffness. The subsequent outcome in the development of a portable quality control device to monitor the acceptable moisture content level to ensure flexibility of the clay mats was also described in this paper.

GPU accelerated algorithms for computing matrix function vector products with applications to exponential integrators and fractional diffusion

The efficient computation of matrix function vector products has become an important area of research in recent times, driven in particular by two important applications: the numerical solution of fractional partial differential equations and the integration of large systems of ordinary differential equations. In this work we consider a problem that combines these two applications, in the form of a numerical solution algorithm for fractional reaction diffusion equations that after spatial discretisation, is advanced in time using the exponential Euler method. We focus on the efficient implementation of the algorithm on Graphics Processing Units (GPU), as we wish to make use of the increased computational power available with this hardware. We compute the matrix function vector products using the contour integration method in [N. Hale, N. Higham, and L. Trefethen. Computing Aα, log(A), and related matrix functions by contour integrals. SIAM J. Numer. Anal., 46(5):2505–2523, 2008]. Multiple levels of preconditioning are applied to reduce the GPU memory footprint and to further accelerate convergence. We also derive an error bound for the convergence of the contour integral method that allows us to pre-determine the appropriate number of quadrature points. Results are presented that demonstrate the effectiveness of the method for large two-dimensional problems, showing a speedup of more than an order of magnitude compared to a CPU-only implementation.

Optimal reactive power flow via the modified barrier Lagrangian function approach

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Optimal reactive power flow via the modified barrier Lagrangian function approach

A new approach called the Modified Barrier Lagrangian Function (MBLF) to solve the Optimal Reactive Power Flow problem is presented. In this approach, the inequality constraints are treated by the Modified Barrier Function (MBF) method, which has a finite convergence property: i.e. the optimal solution in the MBF method can actually be in the bound of the feasible set. Hence, the inequality constraints can be precisely equal to zero. Another property of the MBF method is that the barrier parameter does not need to be driven to zero to attain the solution. Therefore, the conditioning of the involved Hessian matrix is greatly enhanced. In order to show this, a comparative analysis of the numeric conditioning of the Hessian matrix of the MBLF approach, by the decomposition in singular values, is carried out. The feasibility of the proposed approach is also demonstrated with comparative tests to Interior Point Method (IPM) using various IEEE test systems and two networks derived from Brazilian generation/transmission system. The results show that the MBLF method is computationally more attractive than the IPM in terms of speed, number of iterations and numerical conditioning. (C) 2011 Elsevier B.V. All rights reserved.

A SiC-based matrix converter topology for inductive power transfer system

Typical inductive power transfer (IPT) systems employ two power conversion stages to generate a high-frequency primary current from low-frequency utility supply. This paper proposes a matrix-converter-based IPT system, which employs high-speed SiC devices to facilitate the generation of high-frequency current through a single power conversion stage. The proposed matrix converter topology transforms a three-phase low-frequency voltage system to a high-frequency single-phase voltage, which, in turn, powers a series compensated IPT system. A comprehensive mathematical model is developed and power losses are evaluated to investigate the efficiency of the proposed converter topology. Theoretical results are presented with simulations, which are performed in MATLAB/Simulink, in comparison to a conventional two-stage converter. Experimental evident of a prototype IPT system is also presented to demonstrate the applicability of the proposed concept.

A matrix converter based Inductive Power Transfer system

Typical Inductive Power Transfer (IPT) systems employ two power conversion stages to generate a high frequency current from low frequency utility supply. This paper proposes a matrix converter based IPT system that facilitates the generation of high frequency current through a single power conversion stage. The proposed matrix converter topology transforms a 3-phase low frequency voltage system to a high frequency single phase voltage which in turn powers a series compensated IPT system. A comprehensive mathematical model is developed to investigate the behavior of the proposed IPT topology. Theoretical results are presented in comparison to simulations, which are performed in Matlab/ Simulink, to demonstrate the applicability of the proposed concept and the validity of the developed model.

Time evolution of graphene growth on SiC as a function of annealing temperature

We followed by X-ray Photoelectron Spectroscopy (XPS) the time evolution of graphene layers obtained by annealing 3C SiC(111)/Si(111) crystals at different temperatures. The intensity of the carbon signal provides a quantification of the graphene thickness as a function of the annealing time, which follows a power law with exponent 0.5. We show that a kinetic model, based on a bottom-up growth mechanism, provides a full explanation to the evolution of the graphene thickness as a function of time, allowing to calculate the effective activation energy of the process and the energy barriers, in excellent agreement with previous theoretical results. Our study provides a complete and exhaustive picture of Si diffusion into the SiC matrix, establishing the conditions for a perfect control of the graphene growth by Si sublimation.

Processing maps for hot-working of powder metallurgy 1100 Al-10 vol % SiC-particulate metal-matrix composite

The hot-working characteristics of the metal-matrix composite (MMC) Al-10 vol % SiC-particulate (SiCp) powder metallurgy compacts in as-sintered and in hot-extruded conditions were studied using hot compression testing. On the basis of the stress-strain data as a function of temperature and strain rate, processing maps depicting the variation in the efficiency of power dissipation, given by eegr = 2m/(m+1), where m is the strain rate sensitivity of flow stress, have been established and are interpreted on the basis of the dynamic materials model. The as-sintered MMC exhibited a domain of dynamic recrystallization (DRX) with a peak efficiency of about 30% at a temperature of about 500°C and a strain rate of 0.01 s�1. At temperatures below 350°C and in the strain rate range 0.001�0.01 s�1 the MMC exhibited dynamic recovery. The as-sintered MMC was extruded at 500°C using a ram speed of 3 mm s�1 and an extrusion ratio of 10ratio1. A processing map was established on the extruded product, and this map showed that the DRX domain had shifted to lower temperature (450°C) and higher strain rate (1 s�1). The optimum temperature and strain rate combination for powder metallurgy billet conditioning are 500°C and 0.01 s�1, and the secondary metal-working on the extruded product may be done at a higher strain rate of 1 s�1 and a lower temperature of 425°C.

Characteristics of superplasticity domain in the processing map for hot working of an Al alloy 2014---20vol.%Al2O3 metal matrix composite

The processing map for hot working of Al alloy 2014-20vol.%Al2O3 particulate-reinforced cast-plus-extruded composite material has been generated covering the temperature range 300-500 degrees C and the strain rate range 0.001-10 s(-1) based on the dynamic materials model. The efficiency eta of power dissipation given by 2m/(m + 1), where m is the strain rate sensitivity, is plotted as a function of temperature and strain rate to obtain a processing map. A domain of superplasticity has been identified, with a peak efficiency of 62% occurring at 500 degrees C and 0.001 s(-1). The characteristics of this domain have been studied with the help of microstructural evaluation and hot-ductility measurements. Microstructural instability is predicted at higher strain rates above (ls(-1)) and lower temperatures (less than 350 degrees C).

Critical earthquake input power spectral density function models for engineering structures

The problem of determining optimal power spectral density models for earthquake excitation which satisfy constraints on total average power, zero crossing rate and which produce the highest response variance in a given linear system is considered. The solution to this problem is obtained using linear programming methods. The resulting solutions are shown to display a highly deterministic structure and, therefore, fail to capture the stochastic nature of the input. A modification to the definition of critical excitation is proposed which takes into account the entropy rate as a measure of uncertainty in the earthquake loads. The resulting problem is solved using calculus of variations and also within linear programming framework. Illustrative examples on specifying seismic inputs for a nuclear power plant and a tall earth dam are considered and the resulting solutions are shown to be realistic.

Universal power law in crossover from integrability to quantum chaos

We study models of interacting fermions in one dimension to investigate the crossover from integrability to nonintegrability, i.e., quantum chaos, as a function of system size. Using exact diagonalization of finite-sized systems, we study this crossover by obtaining the energy level statistics and Drude weight associated with transport. Our results reinforce the idea that for system size L -> infinity nonintegrability sets in for an arbitrarily small integrability-breaking perturbation. The crossover value of the perturbation scales as a power law similar to L-eta when the integrable system is gapless. The exponent eta approximate to 3 appears to be robust to microscopic details and the precise form of the perturbation. We conjecture that the exponent in the power law is characteristic of the random matrix ensemble describing the nonintegrable system. For systems with a gap, the crossover scaling appears to be faster than a power law.