A systematic method of parameterisation estimation using data assimilation

In numerical weather prediction, parameterisations are used to simulate missing physics in the model. These can be due to a lack of scientific understanding or a lack of computing power available to address all the known physical processes. Parameterisations are sources of large uncertainty in a model as parameter values used in these parameterisations cannot be measured directly and hence are often not well known; and the parameterisations themselves are also approximations of the processes present in the true atmosphere. Whilst there are many efficient and effective methods for combined state/parameter estimation in data assimilation (DA), such as state augmentation, these are not effective at estimating the structure of parameterisations. A new method of parameterisation estimation is proposed that uses sequential DA methods to estimate errors in the numerical models at each space-time point for each model equation. These errors are then fitted to pre-determined functional forms of missing physics or parameterisations that are based upon prior information. We applied the method to a one-dimensional advection model with additive model error, and it is shown that the method can accurately estimate parameterisations, with consistent error estimates. Furthermore, it is shown how the method depends on the quality of the DA results. The results indicate that this new method is a powerful tool in systematic model improvement.

Digital image correlation analysis of the load transfer by implant-supported restorations

This study compared splinted and non-splinted implant-supported prosthesis with and without a distal proximal contact using a digital image correlation method. An epoxy resin model was made with acrylic resin replicas of a mandibular first premolar and second molar and with threaded implants replacing the second premolar and first molar. Splinted and non-splinted metal-ceramic screw-retained crowns were fabricated and loaded with and without the presence of the second molar. A single-camera measuring system was used to record the in-plane deformation on the model surface at a frequency of 1.0 Hz under a load from 0 to 250 N. The images were then analyzed with specialist software to determine the direct (horizontal) and shear strains along the model. Not splinting the crowns resulted in higher stress transfer to the supporting implants when the second molar replica was absent. The presence of a second molar and an effective interproximal contact contributed to lower stress transfer to the supporting structures even for non-splinted restorations. Shear strains were higher in the region between the molars when the second molar was absent, regardless of splinting. The opposite was found for the region between the implants, which had higher shear strain values when the second molar was present. When an effective distal contact is absent, non-splinted implant-supported restorations introduce higher direct strains to the supporting structures under loading. Shear strains appear to be dependent also on the region within the model, with different regions showing different trends in strain changes in the absence of an effective distal contact. (C) 2011 Elsevier Ltd. All rights reserved.

Reappraisal of the effective elastic thickness for the sub-Andes using 3-D finite element flexural modelling, gravity and geological constraints

P>Estimates of effective elastic thickness (T(e)) for the western portion of the South American Plate using, independently, forward flexural modelling and coherence analysis, suggest different thermomechanical properties for the same continental lithosphere. We present a review of these T(e) estimates and carry out a critical reappraisal using a common methodology of 3-D finite element method to solve a differential equation for the bending of a thin elastic plate. The finite element flexural model incorporates lateral variations of T(e) and the Andes topography as the load. Three T(e) maps for the entire Andes were analysed: Stewart & Watts (1997), Tassara et al. (2007) and Perez-Gussinye et al. (2007). The predicted flexural deformation obtained for each T(e) map was compared with the depth to the base of the foreland basin sequence. Likewise, the gravity effect of flexurally induced crust-mantle deformation was compared with the observed Bouguer gravity. T(e) estimates using forward flexural modelling by Stewart & Watts (1997) better predict the geological and gravity data for most of the Andean system, particularly in the Central Andes, where T(e) ranges from greater than 70 km in the sub-Andes to less than 15 km under the Andes Cordillera. The misfit between the calculated and observed foreland basin subsidence and the gravity anomaly for the Maranon basin in Peru and the Bermejo basin in Argentina, regardless of the assumed T(e) map, may be due to a dynamic topography component associated with the shallow subduction of the Nazca Plate beneath the Andes at these latitudes.

Finite temperature effective actions

We present, from first principles, a direct method for evaluating the exact fermion propagator in the presence of a general background held at finite temperature, which can be used to determine the finite temperature effective action for the system. As applications, we determine the complete one loop finite temperature effective actions for (0 + 1)-dimensional QED as well as the Schwinger model. These effective actions, which are derived in the real time (closed time path) formalism, generate systematically all the Feynman amplitudes calculated in thermal perturbation theory and also show that the retarded (advanced) amplitudes vanish in these theories. (C) 2009 Elsevier B.V. All rights reserved.

Internal Residual Stress Measurements in a Bioactive Glass-Ceramic Using Vickers Indentation

The residual stress distribution that arises in the glass matrix during cooling of a partially crystallized 17.2Na(2)O-32.1CaO-48.1SiO(2)-2.5P(2)O(5) (mol%) bioactive glass-ceramic was measured using the Vickers indentation method proposed by Zeng and Rowcliffe (ZR). The magnitude of the determined residual stress at the crystal/glass boundary was 1/4-1/3 of the values measured using X-ray diffraction (within the crystals) and calculated using Selsing`s model. A correction for the crack geometry factor, assuming a semi-elliptical shape, is proposed and then good agreement between experimental and theoretical values is found. Thus, if the actual crack geometry is taken into account, the indentation technique of ZR can be successfully used. In addition, a numerical model for the calculation of residual stresses that takes into account the hemispherical shape of the crystalline precipitates at a free surface was developed. The result is that near the sample surface, the radial component of the residual stress is increased by 70% in comparison with the residual stress calculated by Selsing`s model.

On the generalized eigenvalue method for energies and matrix elements in lattice field theory

We discuss the generalized eigenvalue problem for computing energies and matrix elements in lattice gauge theory, including effective theories such as HQET. It is analyzed how the extracted effective energies and matrix elements converge when the time separations are made large. This suggests a particularly efficient application of the method for which we can prove that corrections vanish asymptotically as exp(-(E(N+1) - E(n))t). The gap E(N+1) - E(n) can be made large by increasing the number N of interpolating fields in the correlation matrix. We also show how excited state matrix elements can be extracted such that contaminations from all other states disappear exponentially in time. As a demonstration we present numerical results for the extraction of ground state and excited B-meson masses and decay constants in static approximation and to order 1/m(b) in HQET.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

A Robust, Fully Adaptive Hybrid Level-Set/Front-Tracking Method for Two-Phase Flows with an Accurate Surface Tension Computation

We present a variable time step, fully adaptive in space, hybrid method for the accurate simulation of incompressible two-phase flows in the presence of surface tension in two dimensions. The method is based on the hybrid level set/front-tracking approach proposed in [H. D. Ceniceros and A. M. Roma, J. Comput. Phys., 205, 391400, 2005]. Geometric, interfacial quantities are computed from front-tracking via the immersed-boundary setting while the signed distance (level set) function, which is evaluated fast and to machine precision, is used as a fluid indicator. The surface tension force is obtained by employing the mixed Eulerian/Lagrangian representation introduced in [S. Shin, S. I. Abdel-Khalik, V. Daru and D. Juric, J. Comput. Phys., 203, 493-516, 2005] whose success for greatly reducing parasitic currents has been demonstrated. The use of our accurate fluid indicator together with effective Lagrangian marker control enhance this parasitic current reduction by several orders of magnitude. To resolve accurately and efficiently sharp gradients and salient flow features we employ dynamic, adaptive mesh refinements. This spatial adaption is used in concert with a dynamic control of the distribution of the Lagrangian nodes along the fluid interface and a variable time step, linearly implicit time integration scheme. We present numerical examples designed to test the capabilities and performance of the proposed approach as well as three applications: the long-time evolution of a fluid interface undergoing Rayleigh-Taylor instability, an example of bubble ascending dynamics, and a drop impacting on a free interface whose dynamics we compare with both existing numerical and experimental data.

Highly Sensitive Fluorescent Method for the Detection of Cholesterol Aldehydes Formed by Ozone and Singlet Molecular Oxygen

Cholesterol oxidation gives rise to a mixture of oxidized products. Different types of products are generated according to the reactive species being involved. Recently, attention has been focused on two cholesterol aldehydes, 3 beta-hydroxy-5 beta-hydroxy-B-norcholestane-6 beta-carboxyaldehyde (1a) and 3 beta-hydroxy-5-oxo-5,6-secocholestan-6-al (1b). These aldehydes can be generated by ozone-, as well as by singlet molecular oxygen-mediated cholesterol oxidation. It has been suggested that 1b is preferentially formed by ozone and la is preferentially formed by singlet molecular oxygen. In this study we describe the use of 1-pyrenebutyric hydrazine (PBH) as a fluorescent probe for the detection of cholesterol aldehydes. The formation of the fluorescent adduct between la with PBH was confirmed by HPLC-MS/MS. The fluorescence spectra of PBH did not change upon binding to the aldehyde. Moreover, the derivatization was also effective in the absence of an acidified medium, which is critical to avoid the formation of cholesterol aldehydes through Hock cleavage of 5 alpha-hydroperoxycholesterol. In conclusion, PBH can be used as an efficient fluorescent probe for the detection/quantification of cholesterol aldehydes in biological samples. Its analysis by HPLC coupled to a fluorescent detector provides a sensitive and specific way to quantify cholesterol aldehydes in the low femtomol range.

EVALUATION OF METAL TOXICITY BY A MODIFIED METHOD BASED ON THE FUNGUS GERRONEMA VIRIDILUCENS BIOLUMINESCENCE IN AGAR MEDIUM

Metal cation toxicity to basidiomycete fungi is poorly understood, despite its well-known importance in terrestrial ecosystems. Moreover, there is no reported methodology for the routine evaluation of metal toxicity to basidiomycetes. In the present study, we describe the development of a procedure to assess the acute toxicity of metal cations (Na(+), K(+), Li(+), Ca(2+), Mg(2+), Co(2+), Zn(2+), Ni(2+), Mn(2+), Cd(2+), and Cu(2+)) to the bioluminescent basidiomycete fungus Gerronema viridilucens. The method is based on the decrease in the intensity of bioluminescence resulting from injuries sustained by the fungus mycelium exposed to either essential or nonessential metal toxicants. The assay described herein enables LIS to propose a metal toxicity series to Gerronenia viridilucens based on data obtained from the bioluminescence intensity (median effective concentration [EC50] values) versus metal concentration: Cd(2+) > Cu(2+) > Mn(2+) approximate to Ni(2+) approximate to Co(2+) > Zn(2+) > Mg(2+) > Li(+) > K(+) approximate to Na(+) > Ca(2+), and to shed some li-ht on the mechanism of toxic action of metal cations to basidiomycete fungi. Environ. Toxicol. Chem. 2010;29:320-326. (C) 2009 SETAC

Metabolic oxidative stress elicited by the copper(II) complex [Cu(isaepy)(2)] triggers apoptosis in SH-SY5Y cells through the induction of the AMP-activated protein kinase/p38(MAPK)/p53 signalling axis: evidence for a combined use with 3-bromopyruvate in neuroblastoma treatment

We have demonstrated previously that the complex bis[(2-oxindol-3-ylimino)-2-(2-aminoethyl)pyridine-N,N`]copper(II), named [Cu(isaepy)(2)], induces AMPK (AMP-activated protein kinase)-dependent/p53-mediated apoptosis in tumour cells by targeting mitochondria. In the present study, we found that p38(MAPK) (p38 mitogen-activated protein kinase) is the molecular link in the phosphorylation cascade connecting AMPK to p53. Transfection of SH-SY5Y cells with a dominant-negative mutant of AMPK resulted in a decrease in apoptosis and a significant reduction in phospho-active p38(MAPK) and p53. Similarly, reverse genetics of p38(MAPK) yielded a reduction in p53 and a decrease in the extent of apoptosis, confirming an exclusive hierarchy of activation that proceeds via AMPK/p38(MAPK)/p53. Fuel supplies counteracted [Cu(isaepy)(2)]-induced apoptosis and AMPK/p38(MAPK)/p53 activation, with glucose being the most effective, suggesting a role for energetic imbalance in [Cu(isaepy)(2)] toxicity. Co-administration of 3BrPA (3-bromopyruvate), a well-known inhibitor of glycolysis, and succinate dehydrogenase, enhanced apoptosis and AMPK/p38(MAPK)/p53 signalling pathway activation. Under these conditions, no toxic effect was observed in SOD (superoxide dismutase)-overexpressing SH-SY5Y cells or in PCNs (primary cortical neurons), which are, conversely, sensitized to the combined treatment with [Cu(isaepy)(2)] and 3BrPA only if grown in low-glucose medium or incubated with the glucose-6-phosphate dehydrogenase inhibitor dehydroepiandrosterone. Overall, the results suggest that NADPH deriving from the pentose phosphate pathway contributes to PCN resistance to [Cu(isaepy)(2)] toxicity and propose its employment in combination with 3BrPA as possible tool for cancer treatment.

Determination of sorbate and benzoate in beverage samples by capillary electrophoresis - Optimization of the method with inspection of ionic mobilities

The aim of this study was to develop a fast capillary electrophoresis method for the determination of benzoate and sorbate ions in commercial beverages. In the method development the pH and constituents of the background electrolyte were selected using the effective mobility versus pH curves. As the high resolution obtained experimentally for sorbate and benzoate in the studies presented in the literature is not in agreement with that expected from the ionic mobility values published, a procedure to determine these values was carried out. The salicylate ion was used as the internal standard. The background electrolyte was composed of 25 mmol L(-1) tris(hydroxymethyl)aminomethane and 12.5 mmol L(-1) 2-hydroxyisobutyric acid, atpH 8.1.Separation was conducted in a fused-silica capillary(32 cm total length and 8.5 cm effective length, 50 mu m I.D.), with short-end injection configuration and direct UV detection at 200 nm for benzoate and salicylate and 254 nm for sorbate ions. The run time was only 28 s. A few figures of merit of the proposed method include: good linearity (R(2) > 0.999), limit of detection of 0.9 and 0.3 mg L(-1) for benzoate and sorbate, respectively, inter-day precision better than 2.7% (n =9) and recovery in the range 97.9-105%. Beverage samples were prepared by simple dilution with deionized water (1:11, v/v). Concentrations in the range of 197-401 mg L(-1) for benzoate and 28-144 mg L(-1) for sorbate were found in soft drinks and tea. (c) 2008 Elsevier B.V. All rights reserved.

Development of a fast capillary electrophoresis method to determine inorganic cations in biodiesel samples

The aim of this study was to develop a fast capillary electrophoresis method for the determination of inorganic cations (Na(+), K(+), Ca(2+), Mg(2+)) in biodiesel samples, using barium (Ba(2+)) as the internal standard. The running electrolyte was optimized through effective mobility curves in order to select the co-ion and Peakmaster software was used to determine electromigration dispersion and buffer capacity. The optimum background electrolyte was composed of 10 mmol L(-1) imidazole and 40 mmol L(-1) of acetic acid. Separation was conducted in a fused-silica capillary (32 cm total length and 23.5 cm effective length, 50 mu m I.D.), with indirect UV detection at 214 nm. The migration time was only 36 s. In order to obtain the optimized conditions for extraction, a fractional factorial experimental design was used. The variables investigated were biodiesel mass, pH, extractant volume, agitation and sonication time. The optimum conditions were: biodiesel mass of 200 mg, extractant volume of 200 mu L. and agitation of 20 min. The method is characterized by good linearity in the concentration range of 0.5-20 mg kg(-1) (r > 0.999), limit of detection was equal to 0.3 mg kg(-1), inter-day precision was equal to 1.88% and recovery in the range of 88.0-120%. The developed method was successfully applied to the determination of cations in biodiesel samples. (c) 2010 Elsevier B.V. All rights reserved.

Development of a fast capillary electrophoresis method for the determination of propranolol-Total analysis time reduction strategies

The aim of this study was to develop a fast capillary electrophoresis method for the determination of propranolol in pharmaceutical preparations. In the method development the pH and constituents of the background electrolyte were selected using the effective mobility versus pH curves. Benzylamine was used as the internal standard. The background electrolyte was composed of 60 mmol L(-1) tris(hydroxymethyl)aminomethane and 30 mmol L(-1) 2-hydroxyisobutyric acid,at pH 8.1. Separation was conducted in a fused-silica capillary (32 cm total length and 8.5 cm effective length, 50 mu m I.D.) with a short-end injection configuration and direct UV detection at 214 nm. The run time was only 14 s. Three different strategies were studied in order to develop a fast CE method with low total analysis time for propranolol analysis: low flush time (Lflush) 35 runs/h, without flush (Wflush) 52 runs/h, and Invert (switched polarity) 45 runs/h. Since the three strategies developed are statistically equivalent, Mush was selected due to the higher analytical frequency in comparison with the other methods. A few figures of merit of the proposed method include: good linearity (R(2) > 0.9999); limit of detection of 0.5 mg L(-1): inter-day precision better than 1.03% (n = 9) and recovery in the range of 95.1-104.5%. (C) 2009 Elsevier B.V. All rights reserved.

Rapid method for the determination. of organic acids in wine by capillary electrophoresis with indirect UV detection

A capillary electrophoresis method for organic acids in wine was developed and validated. The optimal electrolyte consisted of 10 mmol/L 3,5-dinitrobenzoic acid (DNB) at pH 3.6 containing 0.2 mmol/L cetyltrimethylammonium bromide as flow reverser. DNB was chosen because it has an effective mobility similar to the organic acids under investigation, good buffering capacity at pH 3.6, and good chromophoric characteristics for indirect UV-absorbance detection at 254 nm. Sample preparation involved dilution and filtration. The method showed good performance characteristics: Linearity at 6 to 285 mg/L (r > 0.99); detection and quantification limits of 0.64 to 1.55 and 2.12 to 5.15 mg/L, respectively; separation time of less than 5.5 min. Coefficients of variation for ten injections were less than 5% and recoveries varied from 95% to 102%. Application to 23 samples of Brazilian wine confirmed good repeatability and demonstrated wide variation in the organic acid concentrations. (C) 2008 Elsevier Ltd. All rights reserved.