An algorithm for isentropic flow

An efficient algorithm is presented for the solution of the equations of isentropic gas dynamics with a general convex gas law. The scheme is based on solving linearized Riemann problems approximately, and in more than one dimension incorporates operator splitting. In particular, only two function evaluations in each computational cell are required. The scheme is applied to a standard test problem in gas dynamics for a polytropic gas

An efficient flux difference splitting algorithm for unsteady duct flows of a real gas

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a real gas in a single spatial coordinate. This includes flow in a duct of variable cross-section, as well as flow with slab, cylindrical or spherical symmetry, as well as the case of an ideal gas, and can be useful when testing codes for the two-dimensional equations governing compressible flow of a real gas. The resulting scheme requires an average of the flow variables across the interface between cells, and this average is chosen to be the arithmetic mean for computational efficiency, which is in contrast to the usual “square root” averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and for a number of equations of state. The results compare favourably with the results from other schemes.

An efficient shock capturing algorithm for compressible flows in a duct of variable cross-section

A numerical scheme is presented for the solution of the Euler equations of compressible flow of a gas in a single spatial co-ordinate. This includes flow in a duct of variable cross-section as well as flow with slab, cylindrical or spherical symmetry and can prove useful when testing codes for the two-dimensional equations governing compressible flow of a gas. The resulting scheme requires an average of the flow variables across the interface between cells and for computational efficiency this average is chosen to be the arithmetic mean, which is in contrast to the usual ‘square root’ averages found in this type of scheme. The scheme is applied with success to five problems with either slab or cylindrical symmetry and a comparison is made in the cylindrical case with results from a two-dimensional problem with no sources.

An algorithm for the shallow water equations with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the two-dimensional shallow water equations in a generalised coordinate system. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme has good jump capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for flow past a circular obstruction.

An isentropic flow algorithm with body fitted meshes

An efficient algorithm based on flux difference splitting is presented for the solution of the three-dimensional equations of isentropic flow in a generalised coordinate system, and with a general convex gas law. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The algorithm requires only one function evaluation of the gas law in each computational cell. The scheme has good shock capturing properties and the advantage of using body-fitted meshes. Numerical results are shown for Mach 3 flow of air past a circular cylinder. Furthermore, the algorithm also applies to shallow water flows by employing the familiar gas dynamics analogy.

Are differences in life history parameters of the pine beauty moth Panolis flammea modified by host plant quality or gender?

Pine beauty moth (Panolis flammea D&S, Lepidoptera: Noctuidae) were reared individually from egg hatch to pupation on one of three host plants, Pinus sylvestris (native host plant), Pinus contorta (Central Interior seed origin - good quality introduced host) and P. contorta (Alaskan seed origin - poor quality introduced host). After emerging from the pupae the adult moths were confined to a Skeena River seed origin of P. contorta. Female pupal weight and adult life span were significantly higher on P. sylvestris than on the two lodgepole pine seed origins. Development time was, however, not significantly different between treatments, but larval mean relative growth rate was found to be negatively correlated with birth weight and positively correlated with pupal weight. The time to emerge from the pupa was also not significantly different between treatments. However, there were marked differences between the genders. Male moths lost a significantly greater proportion of their weight over the pupal stage but lived significantly longer as adults than the females. Female moths emerged from the pupal stage significantly sooner than male moths. There was no apparent advantage of lai-ge birth size when looked at in terms of subsequent performance. These results are discussed in light of current life history theory.

Self-assembly of PEGylated peptide conjugates containing a modified amyloid β-peptide fragment

The self-assembly of PEGylated peptides containing a modified sequence from the amyloid beta peptide, FEK LVFF, has been studied in aqueous solution. PEG molar masses PEG1k, PEG2k, and PEG10k were used in the conjugates. It is shown that the three FFK LVFF-PEG hybrids form fibrils comprising a FFKLVFF core and a PEG corona. The beta-sheet secondary structure of the peptide is retained in the FFK LVFF fibril core. At sufficiently high concentrations, FEK LVFF-PEG1k and FEK LVFF-PEG2k form a nema tic phase, while PEG10k-FEK LVFF exhibits a hexagonal columnar phase. Simultaneous small angle neutron scattering/shear flow experiments were performed to study the shear flow alignment of the nematic and hexagonal liquid crystal phases. On drying, PEG crystallization occurs without disruption of the FFK LVFF beta-sheet structure leading to characteristic peaks in the X-ray diffraction pattern and FTIR spectra. The stability of beta-sheet structures was also studied in blends of FFKLVFF-PEG conjugates with poly(acrylic acid) (PAA). While PEG crystallization is only observed up to 25% PAA content in the blends, the FFK LVFF beta-sheet structure is retained up to 75% PAA.

A β-amino acid modified heptapeptide containing a designed recognition element disrupts fibrillization of the amyloid β-peptide

We study the complex formation of a peptide betaAbetaAKLVFF, previously developed by our group, with Abeta(1–42) in aqueous solution. Circular dichroism spectroscopy is used to probe the interactions between betaAbetaAKLVFF and Abeta(1–42), and to study the secondary structure of the species in solution. Thioflavin T fluorescence spectroscopy shows that the population of fibers is higher in betaAbetaAKLVFF/Abeta(1–42) mixtures compared to pure Abeta(1–42) solutions. TEM and cryo-TEM demonstrate that co-incubation of betaAbetaAKLVFF with Abeta(1–42) causes the formation of extended dense networks of branched fibrils, very different from the straight fibrils observed for Abeta(1–42) alone. Neurotoxicity assays show that although betaAbetaAKLVFF alters the fibrillization of Abeta(1–42), it does not decrease the neurotoxicity, which suggests that toxic oligomeric Abeta(1–42) species are still present in the betaAbetaAKLVFF/Abeta(1–42) mixtures. Our results show that our designed peptide binds to Abeta(1–42) and changes the amyloid fibril morphology. This is shown to not necessarily translate into reduced toxicity.

Self-assembly of a modified amyloid peptide fragment: pH-responsiveness and nematic phase formation

The self-assembly of peptide YYKLVFFC based on a fragment of the amyloid beta (A) peptide, A beta 16-20, KLVFF has been studied in aqueous solution. The peptide is designed with multiple functional residues to examine the interplay between aromatic interactions and charge on the self-assembly, as well as specific transformations such as the pH-induced phenol-phenolate transition of the tyrosine residue. Circular dichroism (CD) and Fourier-transform infrared (FTIR) spectroscopies are used to investigate the conditions for beta-sheet self-assembly and the role of aromatic interactions in the CD spectrum as a function of pH and concentration. The formation of well-defined fibrils at pH 4.7 is confirmed by cryo-TEM (transmission electron microscope) and negative stain TEM. The morphology changes at higher pH, and aggregates of short twisted fibrils are observed at pH 11. Polarized optical microscopy shows birefringence at a low concentration (1 wt.-%) of YYKLVFFC in aqueous solution, and small-angle X-ray scattering was used to probe nematic phase formation in more detail. A pH-induced transition from nematic to isotropic phases is observed on increasing pH that appears to be correlated to a reduction in aggregate anisotropy upon increasing pH.

Influence of the solvent on the self-assembly of a modified amyloid beta peptide fragment. II. NMR and computer simulation investigation

The conformation of a model peptide AAKLVFF based on a fragment of the amyloid beta peptide A beta 16-20, KLVFF, is investigated in methanol and water via solution NMR experiments and Molecular dynamics computer simulations. In previous work, we have shown that AAKLVFF forms peptide nanotubes in methanol and twisted fibrils in water. Chemical shift measurements were used to investigate the solubility of the peptide as a function of concentration in methanol and water. This enabled the determination of critical aggregation concentrations, The Solubility was lower in water. In dilute solution, diffusion coefficients revealed the presence of intermediate aggregates in concentrated solution, coexisting with NMR-silent larger aggregates, presumed to be beta-sheets. In water, diffusion coefficients did not change appreciably with concentration, indicating the presence mainly of monomers, coexisting with larger aggregates in more concentrated solution. Concentration-dependent chemical shift measurements indicated a folded conformation for the monomers/intermediate aggregates in dilute methanol, with unfolding at higher concentration. In water, an antiparallel arrangement of strands was indicated by certain ROESY peak correlations. The temperature-dependent solubility of AAKLVFF in methanol was well described by a van't Hoff analysis, providing a solubilization enthalpy and entropy. This pointed to the importance of solvophobic interactions in the self-assembly process. Molecular dynamics Simulations constrained by NOE values from NMR suggested disordered reverse turn structures for the monomer, with an antiparallel twisted conformation for dimers. To model the beta-sheet structures formed at higher concentration, possible model arrangements of strands into beta-sheets with parallel and antiparallel configurations and different stacking sequences were used as the basis for MD simulations; two particular arrangements of antiparallel beta-sheets were found to be stable, one being linear and twisted and the other twisted in two directions. These structures Were used to simulate Circular dichroism spectra. The roles of aromatic stacking interactions and charge transfer effects were also examined. Simulated spectra were found to be similar to those observed experimentally.(in water or methanol) which show a maximum at 215 or 218 nm due to pi-pi* interactions, when allowance is made for a 15-18 nm red-shift that may be due to light scattering effects.

Fibrillisation of hydrophobically modified amyloid peptide fragments in an organic solvent

The self-assembly of a hydrophobically modified fragment of the amyloid beta(A beta) peptide has been studied in methanol. The peptide FFKLVFF is based on A beta(16-20) extended at the N terminus by two phenylalanine residues. The formation of amyloid-type fibrils is confirmed by Congo Red staining, thioflavin T fluorescence and circular dichroism experiments. FTIR points to the formation of beta-sheet structures in solution and in dried films and suggests that aggregation occurs at low concentration and is not strongly affected by further increase in concentration, i.e. the peptide is a strong fibril-former in methanol. UV fluorescence experiments on unstained peptide and CD point to the importance of aromatic interactions between phenylalanine groups in driving aggregation into beta-sheets. The CD spectrum differs from that usually observed for beta-sheet assemblies formed by larger peptides or proteins and this is discussed for solutions in methanol and also trifluoroethanol. The fibril structure is imaged by transmission electron microscopy and scanning electron microscopy on dried samples and is confirmed by small-angle X-ray scattering experiments in solution.

Regression based D-optimality experimental design for sparse kernel density estimation

This paper derives an efficient algorithm for constructing sparse kernel density (SKD) estimates. The algorithm first selects a very small subset of significant kernels using an orthogonal forward regression (OFR) procedure based on the D-optimality experimental design criterion. The weights of the resulting sparse kernel model are then calculated using a modified multiplicative nonnegative quadratic programming algorithm. Unlike most of the SKD estimators, the proposed D-optimality regression approach is an unsupervised construction algorithm and it does not require an empirical desired response for the kernel selection task. The strength of the D-optimality OFR is owing to the fact that the algorithm automatically selects a small subset of the most significant kernels related to the largest eigenvalues of the kernel design matrix, which counts for the most energy of the kernel training data, and this also guarantees the most accurate kernel weight estimate. The proposed method is also computationally attractive, in comparison with many existing SKD construction algorithms. Extensive numerical investigation demonstrates the ability of this regression-based approach to efficiently construct a very sparse kernel density estimate with excellent test accuracy, and our results show that the proposed method compares favourably with other existing sparse methods, in terms of test accuracy, model sparsity and complexity, for constructing kernel density estimates.

The TAMORA algorithm: satellite rainfall estimates over West Africa using multi-spectral SEVIRI data

A multi-spectral rainfall estimation algorithm has been developed for the Sahel region of West Africa with the purpose of producing accumulated rainfall estimates for drought monitoring and food security. Radar data were used to calibrate multi-channel SEVIRI data from MSG, and a probability of rainfall at several different rain-rates was established for each combination of SEVIRI radiances. Radar calibrations from both Europe (the SatPrecip algorithm) and Niger (TAMORA algorithm) were used. 10 day estimates were accumulated from SatPrecip and TAMORA and compared with kriged gauge data and TAMSAT satellite rainfall estimates over West Africa. SatPrecip was found to produce large overestimates for the region, probably because of its non-local calibration. TAMORA was negatively biased for areas of West Africa with relatively high rainfall, but its skill was comparable to TAMSAT for the low-rainfall region climatologically similar to its calibration area around Niamey. These results confirm the high importance of local calibration for satellite-derived rainfall estimates. As TAMORA shows no improvement in skill over TAMSAT for dekadal estimates, the extra cloud-microphysical information provided by multi-spectral data may not be useful in determining rainfall accumulations at a ten day timescale. Work is ongoing to determine whether it shows improved accuracy at shorter timescales.