Numerical investigations of linear least squares methods for derivative estimation

The results of a numerical investigation into the errors for least squares estimates of function gradients are presented. The underlying algorithm is obtained by constructing a least squares problem using a truncated Taylor expansion. An error bound associated with this method contains in its numerator terms related to the Taylor series remainder, while its denominator contains the smallest singular value of the least squares matrix. Perhaps for this reason the error bounds are often found to be pessimistic by several orders of magnitude. The circumstance under which these poor estimates arise is elucidated and an empirical correction of the theoretical error bounds is conjectured and investigated numerically. This is followed by an indication of how the conjecture is supported by a rigorous argument.

Compositional model repositories via dynamic constraint satisfaction with order-of-magnitude preferences

J. Keppens and Q. Shen. Compositional model repositories via dynamic constraint satisfaction with order-of-magnitude preferences. Journal of Artificial Intelligence Research, 21:499-550, 2004.

Inorganic arsenic in rice bran and its products are an order of magnitude higher than in bulk grain

Rice is more elevated in arsenic than all other grain crops tested to date, with whole grain (brown) rice having higher arsenic levels than polished (white). It is reported here that rice bran, both commercially purchased and specifically milled for this study, have levels of inorganic arsenic, a nonthreshold, class 1 carcinogen, reaching concentrations of approximately 1 mg/kg dry weight, around 10-20 fold higher than concentrations found in bulk grain. Although pure rice bran is used as a health food supplement, perhaps of more concern is rice bran solubles, which are marketed as a superfood and as a supplement to malnourished children in international aid programs. Five rice bran solubles products were tested, sourced from the United States and Japan, and were found to have 0.61-1.9 mg/kg inorganic arsenic. Manufactures recommend approximately 20 g servings of the rice bran solubles per day, which equates to a 0.012-0.038 mg intake of inorganic arsenic. There are no maximum concentration levels (MCLs) set for arsenic or its species in food stuffs. EU and U.S. water regulations, set at 0.01 mg/L total or inorganic arsenic, respectively, are based on the assumption that 1 L of water per day is consumed, i.e., 0.01 mg of arsenic/ day. At the manufacturers recommended rice bran solubles consumption rate, inorganic arsenic intake exceeds 0.01 mg/ day, remembering that rice bran solubles are targeted at malnourished children and that actual risk is based on mg kg(-1) day(-1) intake.

Carbohydrate-based micelle clusters which enhance hydrophobic drug bioavailability by up to 1 order of magnitude

Amphiphilic chitosan-based polymers (M-w < 20 kDa) self-assemble in aqueous media at low micromolar concentrations to give previously unknown micellar clusters of 100-300 nm in size. Micellar clusters comprise smaller 10-30 nm aggregates, and the nanopolarity/drug incorporation efficiency of their hydrophobic domains can be tailored by varying the degree of lipidic derivatization and molecular weight of the carbohydrate. The extent of drug incorporation by these novel micellar clusters is 1 order of magnitude higher than is seen with triblock copolymers, with molar polymer/drug ratios of 1:48 to 1:67. On intravenous injection, the pharmacodynamic activity of a carbohydrate propofol formulation is increased by 1 order of magnitude when compared to a commercial emulsion formulation, and on topical ocular application of a carbohydrate prednisolone formulation, initial drug aqueous humor levels are similar to those found with a 10-fold dose of prednisolone suspension.

Compositional Model Repositories via Dynamic Constraint Satisfaction with Order-of-Magnitude Preferences

The predominant knowledge-based approach to automated model construction, compositional modelling, employs a set of models of particular functional components. Its inference mechanism takes a scenario describing the constituent interacting components of a system and translates it into a useful mathematical model. This paper presents a novel compositional modelling approach aimed at building model repositories. It furthers the field in two respects. Firstly, it expands the application domain of compositional modelling to systems that can not be easily described in terms of interacting functional components, such as ecological systems. Secondly, it enables the incorporation of user preferences into the model selection process. These features are achieved by casting the compositional modelling problem as an activity-based dynamic preference constraint satisfaction problem, where the dynamic constraints describe the restrictions imposed over the composition of partial models and the preferences correspond to those of the user of the automated modeller. In addition, the preference levels are represented through the use of symbolic values that differ in orders of magnitude.

Logarithms of the values of all vulgar fractions, with numerator and denominator not exceeding 100, arranged in order of magnitude.

"Excerpt Minutes of proceedings of the Institution of Civil Engineers, vol. lxv. Session 1880-81. Part iii."

Methylmercury varies more than one order of magnitude in commercial European rice

P.M. thanks the Royal Thai Government for funding and C.C.B. thanks the School of Natural and Computing Science and PS Analytical for funding.

A mathematical and experimental investigation of microcycle spectral estimates of semiconductor flicker noise

The experimental portion of this thesis tries to estimate the density of the power spectrum of very low frequency semiconductor noise, from 10^-6.3 cps to 1. cps with a greater accuracy than that achieved in previous similar attempts: it is concluded that the spectrum is 1/f^α with α approximately 1.3 over most of the frequency range, but appearing to have a value of about 1 in the lowest decade. The noise sources are, among others, the first stage circuits of a grounded input silicon epitaxial operational amplifier. This thesis also investigates a peculiar form of stationarity which seems to distinguish flicker noise from other semiconductor noise.

In order to decrease by an order of magnitude the pernicious effects of temperature drifts, semiconductor "aging", and possible mechanical failures associated with prolonged periods of data taking, 10 independent noise sources were time-multiplexed and their spectral estimates were subsequently averaged. If the sources have similar spectra, it is demonstrated that this reduces the necessary data-taking time by a factor of 10 for a given accuracy.

In view of the measured high temperature sensitivity of the noise sources, it was necessary to combine the passive attenuation of a special-material container with active control. The noise sources were placed in a copper-epoxy container of high heat capacity and medium heat conductivity, and that container was immersed in a temperature controlled circulating ethylene-glycol bath.

Other spectra of interest, estimated from data taken concurrently with the semiconductor noise data were the spectra of the bath's controlled temperature, the semiconductor surface temperature, and the power supply voltage amplitude fluctuations. A brief description of the equipment constructed to obtain the aforementioned data is included.

The analytical portion of this work is concerned with the following questions: what is the best final spectral density estimate given 10 statistically independent ones of varying quality and magnitude? How can the Blackman and Tukey algorithm which is used for spectral estimation in this work be improved upon? How can non-equidistant sampling reduce data processing cost? Should one try to remove common trands shared by supposedly statistically independent noise sources and, if so, what are the mathematical difficulties involved? What is a physically plausible mathematical model that can account for flicker noise and what are the mathematical implications on its statistical properties? Finally, the variance of the spectral estimate obtained through the Blackman/Tukey algorithm is analyzed in greater detail; the variance is shown to diverge for α ≥ 1 in an assumed power spectrum of k/|f|^α, unless the assumed spectrum is "truncated".

Two orders of magnitude increase in metal piezoresistor sensitivity through nanoscale inhomogenization

Metal-based piezoresistive sensing devices could find a much wider applicability if their sensitivity to mechanical strain could be substantially improved. Here, we report a simple method to enhance the strain sensitivity of metal films by over two orders of magnitude and demonstrate it on specially designed microcantilevers. By locally inhomogenizing thin gold films using controlled electromigration, we have achieved a logarithmic divergence in the strain sensitivity with progressive microstructural modification. The enhancement in strain sensitivity could be explained using non-universal tunneling-percolation transport. We find that the Johnson noise limited signal-to-noise ratio is an order of magnitude better than silicon piezoresistors. This method creates a robust platform for engineering low resistance, high gauge factor metallic piezoresistors that may have profound impact on micro and nanoscale self-sensing technology. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4761817]

Application of a laminar flamelet model to confined explosion hazards

The majority of computational studies of confined explosion hazards apply simple and inaccurate combustion models, requiring ad hoc corrections to obtain realistic flame shapes and often predicting an order of magnitude error in the overpressures. This work describes the application of a laminar flamelet model to a series of two-dimensional test cases. The model is computationally efficient applying an algebraic expression to calculate the flame surface area, an empirical correlation for the laminar flame speed and a novel unstructured, solution adaptive numerical grid system which allows important features of the solution to be resolved close to the flame. Accurate flame shapes are predicted, the correct burning rate is predicted near the walls, and an improvement in the predicted overpressures is obtained. However, in these fully turbulent calculations the overpressures are still too high and the flame arrival times too low, indicating the need for a model for the early laminar burning phase. Due to the computational expense, it is unrealistic to model a laminar flame in the complex geometries involved and therefore a pragmatic approach is employed which constrains the flame to propagate at the laminar flame speed. Transition to turbulent burning occurs at a specified turbulent Reynolds number. With the laminar phase model included, the predicted flame arrival times increase significantly, but are still too low. However, this has no significant effect on the overpressures, which are predicted accurately for a baffled channel test case where rapid transition occurs once the flame reaches the first pair of baffles. In a channel with obstacles on the centreline, transition is more gradual and the accuracy of the predicted overpressures is reduced. However, although the accuracy is still less than desirable in some cases, it is much better than the order of magnitude error previously expected.

Energy inequalities and error estimates for axisymmetric torsion of thin elastic shells of revolution

The problem motivating this investigation is that of pure axisymmetric torsion of an elastic shell of revolution. The analysis is carried out within the framework of the three-dimensional linear theory of elastic equilibrium for homogeneous, isotropic solids. The objective is the rigorous estimation of errors involved in the use of approximations based on thin shell theory.

The underlying boundary value problem is one of Neumann type for a second order elliptic operator. A systematic procedure for constructing pointwise estimates for the solution and its first derivatives is given for a general class of second-order elliptic boundary-value problems which includes the torsion problem as a special case.

The method used here rests on the construction of “energy inequalities” and on the subsequent deduction of pointwise estimates from the energy inequalities. This method removes certain drawbacks characteristic of pointwise estimates derived in some investigations of related areas.

Special interest is directed towards thin shells of constant thickness. The method enables us to estimate the error involved in a stress analysis in which the exact solution is replaced by an approximate one, and thus provides us with a means of assessing the quality of approximate solutions for axisymmetric torsion of thin shells.

Finally, the results of the present study are applied to the stress analysis of a circular cylindrical shell, and the quality of stress estimates derived here and those from a previous related publication are discussed.

Error estimates for a fully discrete spectral scheme for a class of nonlinear, nonlocal dispersive wave equations

We analyze a fully discrete spectral method for the numerical solution of the initial- and periodic boundary-value problem for two nonlinear, nonlocal, dispersive wave equations, the Benjamin–Ono and the Intermediate Long Wave equations. The equations are discretized in space by the standard Fourier–Galerkin spectral method and in time by the explicit leap-frog scheme. For the resulting fully discrete, conditionally stable scheme we prove an L2-error bound of spectral accuracy in space and of second-order accuracy in time.