Reducing Uncertainty in Aeroelastic Flutter Boundaries Using Experimental Data

Flutter prediction as currently practiced is usually deterministic, with a single structural model used to represent an aircraft. By using interval analysis to take into account structural variability, recent work has demonstrated that small changes in the structure can lead to very large changes in the altitude at which
utter occurs (Marques, Badcock, et al., J. Aircraft, 2010). In this follow-up work we examine the same phenomenon using probabilistic collocation (PC), an uncertainty quantification technique which can eficiently propagate multivariate stochastic input through a simulation code,
in this case an eigenvalue-based fluid-structure stability code. The resulting analysis predicts the consequences of an uncertain structure on incidence of
utter in probabilistic terms { information that could be useful in planning
flight-tests and assessing the risk of structural failure. The uncertainty in
utter altitude is confirmed to be substantial. Assuming that the structural uncertainty represents a epistemic uncertainty regarding the
structure, it may be reduced with the availability of additional information { for example aeroelastic response data from a flight-test. Such data is used to update the structural uncertainty using Bayes' theorem. The consequent
utter uncertainty is significantly reduced across the entire Mach number range.

Electron impact excitation of Kr XXVIII

Collision strengths (Ω ) are calculated for all 6328 transitions among the lowest 113 levels belonging to the 2s22p5,2s2p6,2s22p43ℓ,2s2p53ℓ, and 2p63ℓ configurations of fluorine-like krypton, Kr XXVIII, using the Dirac Atomic R -matrix Code. All partial waves with angular momentum J⩽40 are included, sufficient for the convergence of Ω for forbidden transitions. For allowed transitions a top-up is employed to obtain converged values of Ω up to an energy of 400 Ryd. Resonances in the thresholds region are resolved on a narrow energy mesh, and results for effective collision strengths (ϒ) are obtained after averaging the values of Ω over a Maxwellian distribution of electron velocities. Values of ϒ are reported over a wide temperature range below View the MathML source, and the accuracy of the results is assessed. In addition, effective collision strengths are listed for the temperature range View the MathML source, obtained from non-resonant collision strengths generated with the FAC code.

Optical source model for the 23.2-23.6 nm radiation from the multielement germanium soft X-ray laser

Distributions of source intensity in two dimensions (designated the source model), averaged over a single laser pulse, based on experimental measurements of spatial coherence, are considered for radiation from the unresolved 23.2/23.6 nm spectral lines from the germanium collisional X-ray laser. The model derives from measurements of the visibility of Young slit interference fringes determined by a method based on the Wiener-Khinchin theorem. Output from amplifiers comprising three and four target elements have similar coherence properties in directions within the horizontal plane corresponding to strong plasma refraction effects and fitting the coherence data shows source dimensions (FWHM) are similar to 26 mu m (horizontal), significantly smaller than expected by direct imaging, and similar to 125 mu m (vertical: equivalent to the height of the driver excitation). (C) 1999 Elsevier Science B.V. All rights reserved.

Coherence and bandwidth measurements of harmonics generated from solid surfaces irradiated by intense picosecond laser pulses

We present measurements of the transverse and longitudinal coherence lengths of the fourth harmonic of a 1053-nm, 2.5-ps laser generated during high-intensity (up to 10(19) W cm(-2)) interactions with a solid target. Coherence lengths were measured by use of a Young's double-slit interferometer. The effective source size, as defined by the Van Cittert-Zernicke theorem, was found to be 10-12 mu m, and the coherence time was observed to be in the range 0.02-0.4 ps.

Development of an Arcellacean (testate lobose amoebae) based transfer function for sedimentary phosphorous in lakes.

Arcellacea (testate lobose amoebae) communities were assessed from 73 sediment-water interface samples collected from 33 lakes in urban and rural settings within the Greater Toronto Area (GTA), Ontario, Canada, as well as from forested control areas in the Lake Simcoe area, Algonquin Park and eastern Ontario. The results were used to: (1) develop a statistically rigorous arcellacean-based training set for sedimentary phosphorus (Olsen P (OP)) loading; and (2) derive a transfer function to reconstruct OP levels during the post-European settlement era (AD1870s onward) using a chronologically well-constrained core from Haynes Lake on the environmentally sensitive Oak Ridges Moraine, within the GTA. Ordination analysis indicated that OP most influenced arcellacean assemblages, explaining 6.5% (p < 0.005) of total variance. An improved training set where the influence of other important environmental variables (e.g. total organic carbon, total nitrogen, Mg) was reduced, comprised 40 samples from 31 lakes, and was used to construct a transfer function for lacustrine arcellaceans for sedimentary phosphorus (Olsen P) using tolerance downweighted weighted averaging (WA-Tol) with inverse deshrinking (RMSEPjack-77pp; r2jack = 0.68). The inferred reconstruction indicates that OP levels remained near pre-settlement background levels from settlement in the late AD 1970s through to the early AD 1970s. Since OP runoff from both forests and pasture is minimal, early agricultural land use within the lake catchment was as most likely pasture and/or was used to grow perennial crops such as Timothy-grass for hay. A significant increase in inferred OP concentration beginning ~ AD 1972 may have been related to a change in crops (e.g. corn production) in the catchment resulting in more runoff, and the introduction of chemical fertilizers. A dramatic decline in OP after ~ AD 1985 probably corresponds to a reduction in chemical fertilizer use related to advances in agronomy, which permitted a more precise control over required fertilizer application. Another significant increase in OP levels after ~ AD 1995 may have been related to the construction of a large golf course upslope and immediately to the north of Haynes Lake in AD 1993, where significant fertilizer use is required to maintain the fairways. These results demonstrate that arcellaceans have great potential for reconstructing lake water geochemistry and will complement other proxies (e.g. diatoms) in paleolimnological research.

Enhancement of positron annihilation on molecules due to vibrational Feshbach resonances

The zero-range potential model is used to investigate positron collisions and annihilation with molecules. The Kr dimer is considered as an example. It is shown that (i) although positrons do not bind to individual Kr atoms, they do form bound states with Kr. (ii) A sequence of vibrationally excited states of the positron-molecule complex extends into the positron continuum, where it manifests as vibrational Feshbach resonances. (iii) These resonances give a very large contribution to the positron annihilation rate. Even after averaging over the thermal positron energy distribution, the contribution of the lowest Feshbach resonance exceeds that of the non-resonant background by an order of magnitude. © 2002 Elsevier Science B.V. All rights reserved.

A sequential algorithm for sparse support vector classifiers

Support vector machines (SVMs), though accurate, are not preferred in applications requiring high classification speed or when deployed in systems of limited computational resources, due to the large number of support vectors involved in the model. To overcome this problem we have devised a primal SVM method with the following properties: (1) it solves for the SVM representation without the need to invoke the representer theorem, (2) forward and backward selections are combined to approach the final globally optimal solution, and (3) a criterion is introduced for identification of support vectors leading to a much reduced support vector set. In addition to introducing this method the paper analyzes the complexity of the algorithm and presents test results on three public benchmark problems and a human activity recognition application. These applications demonstrate the effectiveness and efficiency of the proposed algorithm.


--------------------------------------------------------------------------------

Market Basket Survey Shows Elevated Levels of As in South Central U.S. Processed Rice Compared to California: Consequences for Human Dietary Exposure:consequences for human dietary exposure

We report the largest market basket survey of arsenic (As) in U.S. rice to date. Our findings show differences in transitional-metal levels between polished and unpolished rice and geographical variation in As and selenium (Se) between rice processed in California and the South Central U.S. The mean and median As grain levels for the South Central U.S. were 0.30 and 0.27 µg As g-1, respectively, for 107 samples. Levels for California were 41% lower than the South Central U.S., with a mean of 0.17 µg As g-1 and a median of 0.16 µg As g-1 for 27 samples. The mean and median Se grain levels for the South Central U.S. were 0.19 µg Se g-1. Californian rice levels were lower, averaging only 0.08 and 0.06 µg Se g-1 for mean and median values, respectively. The difference between the two regions was found to be significant for As and Se (General Linear Model (GLM):? As p < 0.001; Se p < 0.001). No statistically significant differences were observed in As or Se levels between polished and unpolished rice (GLM:? As p = 0.213; Se p = 0.113). No significant differences in grain levels of manganese (Mn), cobalt (Co), copper (Cu), or zinc (Zn) were observed between California and the South Central U.S. Modeling arsenic intake for the U.S. population based on this survey shows that for certain groups (namely Hispanics, Asians, sufferers of Celiac disease, and infants) dietary exposure to inorganic As from elevated levels in rice potentially exceeds the maximum intake of As from drinking water (based on consumption of 1 L of 0.01 mg L-1 In. As) and Californian state exposure limits. Further studies on the transformation of As in soil, grain As bioavailability in the human gastrointestinal tract, and grain elemental speciation trends are critical.

Hardware elliptic curve cryptographic processor over GF(p)

A novel hardware architecture for elliptic curve cryptography (ECC) over GF(p) is introduced. This can perform the main prime field arithmetic functions needed in these cryptosystems including modular inversion and multiplication. This is based on a new unified modular inversion algorithm that offers considerable improvement over previous ECC techniques that use Fermat's Little Theorem for this operation. The processor described uses a full-word multiplier which requires much fewer clock cycles than previous methods, while still maintaining a competitive critical path delay. The benefits of the approach have been demonstrated by utilizing these techniques to create a field-programmable gate array (FPGA) design. This can perform a 256-bit prime field scalar point multiplication in 3.86 ms, the fastest FPGA time reported to date. The ECC architecture described can also perform four different types of modular inversion, making it suitable for use in many different ECC applications. © 2006 IEEE.

Stellar Surface Magneto-Convection as a Source of Astrophysical Noise. I. Multi-component Parameterisation of Absorption Line Profiles

We outline our techniques to characterise photospheric granulation as an astrophysical noise source. A four component parameterisation of granulation is developed that can be used to reconstruct stellar line asymmetries and radial velocity shifts due to photospheric convective motions. The four components are made up of absorption line profiles calculated for granules, magnetic intergranular lanes, non-magnetic intergranular lanes, and magnetic bright points at disc centre. These components are constructed by averaging Fe I $6302 \mathrm{\AA}$ magnetically sensitive absorption line profiles output from detailed radiative transport calculations of the solar photosphere. Each of the four categories adopted are based on magnetic field and continuum intensity limits determined from examining three-dimensional magnetohydrodynamic simulations with an average magnetic flux of $200 \mathrm{G}$. Using these four component line profiles we accurately reconstruct granulation profiles, produced from modelling 12 x 12 Mm$^2$ areas on the solar surface, to within $\sim \pm$ 20 cm s$^{-1}$ on a $\sim$ 100 m s$^{-1}$ granulation signal. We have also successfully reconstructed granulation profiles from a $50 \mathrm{G}$ simulation using the parameterised line profiles from the $200 \mathrm{G}$ average magnetic field simulation. This test demonstrates applicability of the characterisation to a range of magnetic stellar activity levels.

A multi-fibre multi-layer representative volume element (M2RVE) for prediction of matrix and interfacial damage in composite laminates

Multiscale micro-mechanics theory is extensively used for the prediction of the material response and damage analysis of unidirectional lamina using a representative volume element (RVE). Th is paper presents a RVE-based approach to characterize the materi al response of a multi-fibre cross-ply laminate considering the effect of matrix damage and fibre-matrix interfacial strength. The framework of the homogenization theory for periodic media has been used for the analysis of a 'multi-fibre multi-layer representative volume element' (M2 RVE) representing cross-ply laminate. The non-homogeneous stress-strain fields within the M2RVE are related to the average stresses and strains by using Gauss theorem and the Hill-Mandal strain energy equivalence principle. The interfacial bonding strength affects the in-plane shear stress-strain response significantl y. The material response predicted by M2 RVE is in good agreement with the experimental results available in the literature. The maximum difference between the shear stress predicted using M2 RVE and the experimental results is ~15% for the bonding strength of 30MPa at the strain value of 1.1%

On the required shape corrections to the local density and generalized gradient approximations to the Kohn-Sham potentials for molecular response calculations of (hyper)polarizabilities and excitation energies

It is well known that shape corrections have to be applied to the local-density (LDA) and generalized gradient (GGA) approximations to the Kohn-Sham exchange-correlation potential in order to obtain reliable response properties in time dependent density functional theory calculations. Here we demonstrate that it is an oversimplified view that these shape corrections concern primarily the asymptotic part of the potential, and that they affect only Rydberg type transitions. The performance is assessed of two shape-corrected Kohn-Sham potentials, the gradient-regulated asymptotic connection procedure applied to the Becke-Perdew potential (BP-GRAC) and the statistical averaging of (model) orbital potentials (SAOP), versus LDA and GGA potentials, in molecular response calculations of the static average polarizability alpha, the Cauchy coefficient S-4, and the static average hyperpolarizability beta. The nature of the distortions of the LDA/GGA potentials is highlighted and it is shown that they introduce many spurious excited states at too low energy which may mix with valence excited states, resulting in wrong excited state compositions. They also lead to wrong oscillator strengths and thus to a wrong spectral structure of properties like the polarizability. LDA, Becke-Lee-Yang-Parr (BLYP), and Becke-Perdew (BP) characteristically underestimate contributions to alpha and S-4 from bound Rydberg-type states and overestimate those from the continuum. Cancellation of the errors in these contributions occasionally produces fortuitously good results. The distortions of the LDA, BLYP, and BP spectra are related to the deficiencies of the LDA/GGA potentials in both the bulk and outer molecular regions. In contrast, both SAOP and BP-GRAC potentials produce high quality polarizabilities for 21 molecules and also reliable Cauchy moments and hyperpolarizabilities for the selected molecules. The analysis for the N-2 molecule shows, that both SAOP and BP-GRAC yield reliable energies omega(i) and oscillator strengths f(i) of individual excitations, so that they reproduce well the spectral structure of alpha and S-4.(C) 2002 American Institute of Physics.

Exchange potential from the common energy denominator approximation for the Kohn-Sham Green's function: Application to (hyper)polarizabilities of molecular chains

An approximate Kohn-Sham (KS) exchange potential v(xsigma)(CEDA) is developed, based on the common energy denominator approximation (CEDA) for the static orbital Green's function, which preserves the essential structure of the density response function. v(xsigma)(CEDA) is an explicit functional of the occupied KS orbitals, which has the Slater v(Ssigma) and response v(respsigma)(CEDA) potentials as its components. The latter exhibits the characteristic step structure with "diagonal" contributions from the orbital densities \psi(isigma)\(2), as well as "off-diagonal" ones from the occupied-occupied orbital products psi(isigma)psi(j(not equal1)sigma). Comparison of the results of atomic and molecular ground-state CEDA calculations with those of the Krieger-Li-Iafrate (KLI), exact exchange (EXX), and Hartree-Fock (HF) methods show, that both KLI and CEDA potentials can be considered as very good analytical "closure approximations" to the exact KS exchange potential. The total CEDA and KLI energies nearly coincide with the EXX ones and the corresponding orbital energies epsilon(isigma) are rather close to each other for the light atoms and small molecules considered. The CEDA, KLI, EXX-epsilon(isigma) values provide the qualitatively correct order of ionizations and they give an estimate of VIPs comparable to that of the HF Koopmans' theorem. However, the additional off-diagonal orbital structure of v(xsigma)(CEDA) appears to be essential for the calculated response properties of molecular chains. KLI already considerably improves the calculated (hyper)polarizabilities of the prototype hydrogen chains H-n over local density approximation (LDA) and standard generalized gradient approximations (GGAs), while the CEDA results are definitely an improvement over the KLI ones. The reasons of this success are the specific orbital structures of the CEDA and KLI response potentials, which produce in an external field an ultranonlocal field-counteracting exchange potential. (C) 2002 American Institute of Physics.

Shape corrections to exchange-correlation potentials by gradient-regulated seamless connection of model potentials for inner and outer region

Shape corrections to the standard approximate Kohn-Sham exchange-correlation (xc) potentials are considered with the aim to improve the excitation energies (especially for higher excitations) calculated with time-dependent density functional perturbation theory. A scheme of gradient-regulated connection (GRAC) of inner to outer parts of a model potential is developed. Asymptotic corrections based either on the potential of Fermi and Amaldi or van Leeuwen and Baerends (LB) are seamlessly connected to the (shifted) xc potential of Becke and Perdew (BP) with the GRAC procedure, and are employed to calculate the vertical excitation energies of the prototype molecules N-2, CO, CH2O, C2H4, C5NH5, C6H6, Li-2, Na-2, K-2. The results are compared with those of the alternative interpolation scheme of Tozer and Handy as well as with the results of the potential obtained with the statistical averaging of (model) orbital potentials. Various asymptotically corrected potentials produce high quality excitation energies, which in quite a few cases approach the benchmark accuracy of 0.1 eV for the electronic spectra. Based on these results, the potential BP-GRAC-LB is proposed for molecular response calculations, which is a smooth potential and a genuine "local" density functional with an analytical representation. (C) 2001 American Institute of Physics.

Model categories for orthogonal calculus

We restate the notion of orthogonal calculus in terms of model categories. This provides a cleaner set of results and makes the role of O(n)-equivariance clearer. Thus we develop model structures for the category of n-polynomial and n-homogeneous functors, along with Quillen pairs relating them. We then classify n-homogeneous functors, via a zig-zag of Quillen equivalences, in terms of spectra with an O(n)-action. This improves upon the classification theorem of Weiss. As an application, we develop a variant of orthogonal calculus by replacing topological spaces with orthogonal spectra.