Model reduction for nonlinear aerodynamics and aeroelasticity using a Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of Computational Fluid Dynamics (CFD) solvers, for fast nonlinear aerodynamic and aeroelastic modeling. A nonlinear function is identified on selected interpolation points by
a discrete empirical interpolation method (DEIM). The flow field is then reconstructed using a least square approximation of the flow modes extracted
by proper orthogonal decomposition (POD). The aeroelastic reduce order
model (ROM) is completed by introducing a nonlinear mapping function
between displacements and the DEIM points. The proposed model is investigated to predict the aerodynamic forces due to forced motions using
a N ACA 0012 airfoil undergoing a prescribed pitching oscillation. To investigate aeroelastic problems at transonic conditions, a pitch/plunge airfoil
and a cropped delta wing aeroelastic models are built using linear structural models. The presence of shock-waves triggers the appearance of limit
cycle oscillations (LCO), which the model is able to predict. For all cases
tested, the new ROM shows the ability to replicate the nonlinear aerodynamic forces, structural displacements and reconstruct the complete flow
field with sufficient accuracy at a fraction of the cost of full order CFD
model.

Model reduction for nonlinear aeroelasticity using the Discrete Empirical Interpolation Method

A novel surrogate model is proposed in lieu of computational fluid dynamic (CFD) code for fast nonlinear aerodynamic modeling. First, a nonlinear function is identified on selected interpolation points defined by discrete empirical interpolation method (DEIM). The flow field is then reconstructed by a least square approximation of flow modes extracted by proper orthogonal decomposition (POD). The proposed model is applied in the prediction of limit cycle oscillation for a plunge/pitch airfoil and a delta wing with linear structural model, results are validate against a time accurate CFD-FEM code. The results show the model is able to replicate the aerodynamic forces and flow fields with sufficient accuracy while requiring a fraction of CFD cost.

Design of silicon IP cores for biorthogonal wavelet transforms

A methodology for rapid silicon design of biorthogonal wavelet transform systems has been developed. This is based on generic, scalable architectures for the forward and inverse wavelet filters. These architectures offer efficient hardware utilisation by combining the linear phase property of biorthogonal filters with decimation and interpolation. The resulting designs have been parameterised in terms of types of wavelet and wordlengths for data and coefficients. Control circuitry is embedded within these cores that allows them to be cascaded for any desired level of decomposition without any interface logic. The time to produce silicon designs for a biorthogonal wavelet system is only the time required to run synthesis and layout tools with no further design effort required. The resulting silicon cores produced are comparable in area and performance to hand-crafted designs. These designs are also portable across a range of foundries and are suitable for FPGA and PLD implementations.

Efficient modeling of novel uniplanar left-handed metamaterials

This paper presents an efficient. modeling technique for the derivation of the dispersion characteristics of novel uniplanar metallodielectric periodic structures. The analysis is based on the method of moments and an interpolation scheme, which significantly accelerates the computations. Triangular basis functions are used that allow for modeling of arbitrary shaped metallic elements. Based on this method, novel uniplanar left-handed (LH) metamaterials are proposed. Variations of the split rectangular-loop element printed on grounded dielectric substrate are demonstrated to possess LH propagation properties. Full-wave dispersion curves are presented. Based on the dual transmission-line concept, we study the distribution of the modal fields And the variation of series capacitance and shunt inductance for all the proposed elements. A verification of the left-handedness is presented by means of full-wave simulation of finite uniplanar arrays using commercial software (HFSS). The cell dimensions are a small fraction of the wavelength (approximately lambda/24) so that the structures can he considered as a homogeneous effective medium. The structures are simple, readily scalable to higher frequencies, and compatible with low-cost fabrication techniques.

Efficient modelling technique for fractal electromagnetic band-gap arrays

An efficient modelling technique is proposed for the analysis of a fractal-element electromagnetic band-gap array. The modelling is based on a method of moments modal analysis in conjunction with an interpolation scheme, which significantly accelerates the computations. The plane-wave and the surface-wave responses of the structure have been studied by means of transmission coefficients and dispersion diagrams. The multiband properties and the compactness of the proposed structure are presented. The technique is general and can be applied to arbitrary-shaped element geometries.

High precision C-14 dating of Holocene peat deposits: A comparison of Bayesian calibration and wiggle-matching approaches

The chronologies of five northern European ombrotrophic peat bogs subjected to a large ANIS C-14 dating effort (32-44 dates/site) are presented here. The results of Bayesian calibration (BCal) of dates with a prior assumption of chronological ordering were compared with a Bayesian wiggle-match approach (Bpeat) which assumes constant linear accumulation over sections of the peat profile. Interpolation of BCal age estimates of dense sequences of C-14 dates showed variable patterns of peat accumulation with time, with changes in accumulation occurring at intervals ranging from 20 to 50 cm. Within these intervals, peat accumulation appeared to be relatively linear. Close analysis suggests that some of the inferred variations in accumulation rate were related to the plant macrofossil composition of the peat. The wiggle-matched age-depth models had relatively high chronological uncertainty within intervals of closely spaced 14 C dates, suggesting that the premise of constant linear accumulation over large sections of the peat profile is unrealistic. Age models based on the assumption of linear accumulation over large parts of a peat core (and therefore only effective over millennial timescales), are not compatible with studies examining environmental change during the Holocene, where variability often occurs at decadal to centennial time-scales. Ideally, future wiggle-match age models should be constrained, with boundaries between sections based on the plant macrofossil composition of the peat and physical-chemical parameters such as the degree of decomposition. Strategies for the selection of material for dating should be designed so that there should be enough C-14 dates to accurately reconstruct the peat accumulation rate of each homogeneous stratigraphic unit. (c) 2006 Elsevier Ltd. All rights reserved.

Advanced glycation end products cause increased CCN family and extracellular matrix gene expression in the diabetic rodent retina

Aims/hypothesis: Referred to as CCN, the family of growth factors consisting of cystein-rich protein 61 (CYR61, also known as CCN1), connective tissue growth factor (CTGF, also known as CCN2), nephroblastoma overexpressed gene (NOV, also known as CCN3) and WNT1-inducible signalling pathway proteins 1, 2 and 3 (WISP1, -2 and -3; also known as CCN4, -5 and -6) affects cellular growth, differentiation, adhesion and locomotion in wound repair, fibrotic disorders, inflammation and angiogenesis. AGEs formed in the diabetic milieu affect the same processes, leading to diabetic complications including diabetic retinopathy. We hypothesised that pathological effects of AGEs in the diabetic retina are a consequence of AGE-induced alterations in CCN family expression.

Materials and methods: CCN gene expression levels were studied at the mRNA and protein level in retinas of control and diabetic rats using real-time quantitative PCR, western blotting and immunohistochemistry at 6 and 12 weeks of streptozotocin-induced diabetes in the presence or absence of aminoguanidine, an AGE inhibitor. In addition, C57BL/6 mice were repeatedly injected with exogenously formed AGE to establish whether AGE modulate retinal CCN growth factors in vivo.

Results: After 6 weeks of diabetes, Cyr61 expression levels were increased more than threefold. At 12 weeks of diabetes, Ctgf expression levels were increased twofold. Treatment with aminoguanidine inhibited Cyr61 and Ctgf expression in diabetic rats, with reductions of 31 and 36%, respectively, compared with untreated animals. Western blotting showed a twofold increase in CTGF production, which was prevented by aminoguanidine treatment. In mice infused with exogenous AGE, Cyr61 expression increased fourfold and Ctgf expression increased twofold in the retina.

Conclusions/interpolation: CTGF and CYR61 are downstream effectors of AGE in the diabetic retina, implicating them as possible targets for future intervention strategies against the development of diabetic retinopathy.

Dynamical instability in surface permeability characteristics of building sandstones in response to salt accumulation over time

This paper explores how the surface permeability of sandstone blocks changes over time in response to repeated salt weathering cycles. Surface permeability controls the amount of moisture and dissolved salt that can penetrate in and facilitate decay. Connected pores permit the movement of moisture (and hence soluble salts) into the stone interior, and where areas are more or less permeable soluble salts may migrate along preferred pathways at differential rates. Previous research has shown that salts can accumulate in the near-surface zone and lead to partial pore blocking which influences subsequent moisture ingress and causes rapid salt accumulation in the near-surface zone.

Two parallel salt weathering simulations were carried out on blocks of Peakmoor Sandstone of different volumes. Blocks were removed from simulations after 2, 5, 10, 20 and 60 cycles. Permeability measurements were taken for these blocks at a resolution of 20 mm, providing a grid of 100 permeability values for each surface. The geostatistical technique of ordinary kriging was applied to the data to produce a smoothed interpolation of permeability for these surfaces, and hence improve understanding of the evolution of permeability over time in response to repeated salt weathering cycles.

Results illustrate the different responses of the sandstone blocks of different volumes to repeated salt weathering cycles. In both cases, after an initial subtle decline in the permeability (reflecting pore blocking), the permeability starts to increase — reflected in a rise in mean, maximum and minimum values. However, between 10 and 20 cycles, there is a jump in the mean and range permeability of the group A block surfaces coinciding with the onset of meaningful debris release. After 60 cycles, the range of permeability in the group A block surface had increased markedly, suggesting the development of a secondary permeability. The concept of dynamic instability and divergent behaviour is applied at the scale of a single block surface, with initial small-scale differences across a surface having larger scale consequences as weathering progresses.

After cycle 10, group B blocks show a much smaller increase in mean permeability, and the range stays relatively steady — this may be explained by the capillary conditions set up by the smaller volume of the stone, allowing salts to migrate to the ‘back’ of the blocks and effectively relieving stress at the ‘front’ face.

Phylogeographic analysis of North American populations of the parasitic herbaceous plant Monotropa hypopitys reveals a complex history of range expansion from multiple late glacial refugia

Aim We carried out a phylogeographic study across the range of the herbaceous plant species Monotropa hypopitys L. in North America to determine whether its current disjunct distribution is due to recolonization from separate eastern and western refugia after the Last Glacial Maximum (LGM). Location North America: Pacific Northwest and north-eastern USA/south-eastern Canada. Methods Palaeodistribution modelling was carried out to determine suitable climatic regions for M. hypopitys at the LGM. We analysed between 155 and 176 individuals from 39 locations spanning the species' entire range in North America. Sequence data were obtained for the chloroplast rps2 gene (n=168) and for the nuclear ITS region (n=158). Individuals were also genotyped for eight microsatellite loci (n=176). Interpolation of diversity values was used to visualize the range-wide distribution of genetic diversity for each of the three marker classes. Minimum spanning networks were constructed showing the relationships between the rps2 and ITS haplotypes, and the geographical distributions of these haplotypes were plotted. The numbers of genetic clusters based on the microsatellite data were estimated using Bayesian clustering approaches. Results The palaeodistribution modelling indicated suitable climate envelopes for M. hypopitys at the LGM in both the Pacific Northwest and south-eastern USA. High levels of genetic diversity and endemic haplotypes were found in Oregon, the Alexander Archipelago, Wisconsin, and in the south-eastern part of the species' distribution range. Main conclusions Our results suggest a complex recolonization history for M. hypopitys in North America, involving persistence in separate eastern and western refugia. A generally high degree of congruence between the different marker classes analysed indicated the presence of multiple refugia, with at least two refugia in each area. In the west, putative refugia were identified in Oregon and the Alexander Archipelago, whereas eastern refugia may have been located in the southern part of the species' current distribution, as well as in the 'Driftless Area'. These findings are in contrast to a previous study on the related species Orthilia secunda, which has a similar disjunct distribution to M. hypopitys, but which appears to have recolonized solely from western refugia. © 2011 Blackwell Publishing Ltd.

Transonic Aeroelastic Stability Analysis Using a Kriging-Based Schur Complement Formulation

A method is described to allow searches for transonic aeroelastic instability of realistically sized aircraft models in multidimensional parameter spaces when computational fluid dynamics are used to model the aerodynamics. Aeroelastic instability is predicted from a small nonlinear eigenvalue problem. The approximation of the computationally expensive interaction term modeling the fluid response is formulated to allow the automated and blind search for aeroelastic instability. The approximation uses a kriging interpolation of exact numerical samples covering the parameter space. The approach, demonstrated for the Goland wing and the multidisciplinary optimization transport wing, results in stability analyses over whole flight envelopes at an equivalent cost of several steady-state simulations.

Technical Note: Assessing a 24/7 solution for monitoring water quality loads in small river catchments

Quantifying nutrient and sediment loads in catchments is dif?cult owing to diffuse controls related to storm hydrology. Coarse sampling and interpolation methods are prone to very high uncertainties due to under-representation of high discharge, short duration events. Additionally, important low-?ow processes such as diurnal signals linked to point source impacts are missed. Here we demonstrate a solution based on a time-integrated approach to sampling with a standard 24 bottle autosampler con?gured to take a sample every 7 h over a week according to a Plynlimon design. This is evaluated with a number of other sampling strategies using a two-year dataset of sub-hourly discharge and phosphorus concentration data. The 24/7 solution is shown to be among the least uncertain in estimating load (inter-quartile range: 96% to 110% of actual load in year 1 and 97% to 104% in year 2) due to the increased frequency raising the probability of sampling storm events and point source signals. The 24/7 solution would appear to be most parsimonious in terms of data coverage and certainty, process signal representation, potential laboratory commitment, technology requirements and the ability to be widely deployed in complex catchments.

A lower mass for the exoplanet WASP-21b

We present high-precision transit observations of the exoplanet WASP-21b, obtained with the Rapid Imager to Search for Exoplanets instrument mounted on the 2.0-m Liverpool Telescope. A transit model is fitted, coupled with a Markov chain Monte Carlo routine, to derive accurate system parameters. The two new high-precision transits allow us to estimate the stellar density directly from the light curve. Our analysis suggests that WASP-21 is evolving off the main sequence which led to a previous overestimation of the stellar density. Using isochrone interpolation, we find a stellar mass of 0.86 ± 0.04 Msun, which is significantly lower than previously reported (1.01 ± 0.03 Msun). Consequently, we find a lower planetary mass of 0.27 ± 0.01 MJup. A lower inclination (87?4 ± 0?3) is also found for the system than previously reported, resulting in a slightly larger stellar (R*= 1.10 ± 0.03 Rsun) and planetary radius (Rp= 1.14 ± 0.04 RJup). The planet radius suggests a hydrogen/helium composition with no core which strengthens the correlation between planetary density and host star metallicity. A new ephemeris is determined for the system, i.e. T0= 245 5084.519 74 ± 0.000 20 (HJD) and P= 4.322 5060 ± 0.000 0031 d. We found no transit timing variations in WASP-21b.

Electronic properties of zircon and hafnon from many-body perturbation theory

The electronic properties of zircon and hafnon, two wide-gap high-kappa materials, are investigated using many-body perturbation theory (MBPT) combined with the Wannier interpolation technique. For both materials, the calculated band structures differ from those obtained within density-functional theory and MBPT by (i) a slight displacement of the highest valence-band maximum from the Gamma point and (ii) an opening of the indirect band gap to 7.6 and 8.0 eV for zircon and hafnon, respectively. The introduction of vertex corrections in the many-body self-energy does not modify the results except for a global rigid shift of the many-body corrections.

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.