A boundary preserving numerical algorithm for the Wright-Fisher model with mutation

The Wright-Fisher model is an Itô stochastic differential equation that was originally introduced to model genetic drift within finite populations and has recently been used as an approximation to ion channel dynamics within cardiac and neuronal cells. While analytic solutions to this equation remain within the interval [0,1], current numerical methods are unable to preserve such boundaries in the approximation. We present a new numerical method that guarantees approximations to a form of Wright-Fisher model, which includes mutation, remain within [0,1] for all time with probability one. Strong convergence of the method is proved and numerical experiments suggest that this new scheme converges with strong order 1/2. Extending this method to a multidimensional case, numerical tests suggest that the algorithm still converges strongly with order 1/2. Finally, numerical solutions obtained using this new method are compared to those obtained using the Euler-Maruyama method where the Wiener increment is resampled to ensure solutions remain within [0,1].

Anodization of Aluminium using a fast two-step process

Ultra-fast two-step anodization method is developed for obtaining ordered nano-pores on aluminium (Al) foil. First anodization was carried out for 10 min, followed by 3 min of second anodization at high voltage (150 V) compared to previous reports of anodization times of 12 h (40-60 V). The pore dimensions on anodized alumina are 180 nm for pore diameter and 130 nm for inter-pore distance. It was evident that by increasing the anodization voltage to 150 V, the diameter of the pores formed was above 150 nm. The electrolyte and its temperature affect the shape and size of the pore formation. At lower anodization temperature, controlled pore formation was observed. The anodized samples were characterized using the field emission scanning electron microscope (FE-SEM) to determine the pore diameter and inter-pore distance. Using UV-Visible spectroscopy, the reflectance spectra of anodized samples were measured. The alumina (Al2O3) peaks were identified by x-ray diffraction (XRD) technique. The x-ray photo electron spectroscopy (XPS) analysis confirmed the Al 2p peak at 73.1 eV along with the oxygen O 1s at 530.9 eV and carbon traces C 1s at 283.6 eV.

较高阶非线性抽运-探测过程中的空间效应

在飞秒抽运一探测光谱技术中，空间分辨的探测光信号反映了在不同空间位置的材料的非线性效应。当抽运光强度增大时，探测光信号中会出现明显的高阶特别是五阶非线性效应。利用劈裂算子方法直接解决了一维非线性传播方程的问题。在数值模拟中，研究了在不同抽运强度和位置下的抽运一探测过程中的五阶非线性效应。在足够高的抽运场下，探测信号出现清晰的振荡，显示了三阶和五阶非线性效应之间的干涉。当空间位置离抽运场中心足够远时，五阶比三阶非线性效应的衰减快得多，对其物理机制和趋势进行了定性的讨论。

An Arbitrary Lagrangian Eulerian Method for Three-Phase Flows with Triple Junction Points

We present a moving mesh method suitable for solving two-dimensional and axisymmetric three-liquid flows with triple junction points. This method employs a body-fitted unstructured mesh where the interfaces between liquids are lines of the mesh system, and the triple junction points (if exist) are mesh nodes. To enhance the accuracy and the efficiency of the method, the mesh is constantly adapted to the evolution of the interfaces by refining and coarsening the mesh locally; dynamic boundary conditions on interfaces, in particular the triple points, are therefore incorporated naturally and accurately in a Finite- Element formulation. In order to allow pressure discontinuity across interfaces, double-values of pressure are necessary for interface nodes and triple-values of pressure on triple junction points. The resulting non-linear system of mass and momentum conservation is then solved by an Uzawa method, with the zero resultant condition on triple points reinforced at each time step. The method is used to investigate the rising of a liquid drop with an attached bubble in a lighter liquid.

Facile synthesis of highly uniform octahedral LuVO4 microcrystals by a facile chemical conversion method

Uniform octahedral LuVO4 microcrystals have been successfully prepared through a designed two-step hydrothermal method. One-dimensional lutetium precursor was first prepared through a simple hydrothermal route. Subsequently, a well-shaped octahedral LuVO4 sample was synthesized at the expense of the wirelike precursors during the hydrothermal process. The whole process in this method was carried out in aqueous conditions without the use of any organic solvents, surfactant, or catalyst. The conversion process from nanowire precursor to octahedral product has been investigated in detail. The LuVO4 : Ln(3+) (Ln Eu, Dy, Sm, and Er) phosphors show strong light emissions with different colors coming from different activator ions under ultraviolet light excitation or low-voltage electron beam excitation. Furthermore, this general and facile method may be of much significance in the synthesis of many other lanthanide compounds with polyhedral morphology.

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

Predictor-Corrector Primal-Dual Interior Point Method for Solving Economic Dispatch Problems: A Postoptimization Analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Bose-Einstein condensation dynamics from the numerical solution of the Gross-Pitaevskii equation

We study certain stationary and time-evolution problems of trapped Bose-Einstein condensates using the numerical solution of the Gross-Pitaevskii (GP) equation with both spherical and axial symmetries. We consider time-evolution problems initiated by suddenly changing the interatomic scattering length or harmonic trapping potential in a stationary condensate. These changes introduce oscillations in the condensate which are studied in detail. We use a time iterative split-step method for the solution of the time-dependent GP equation, where all nonlinear and linear non-derivative terms are treated separately from the time propagation with the kinetic energy terms. Even for an arbitrarily strong nonlinear term this leads to extremely accurate and stable results after millions of time iterations of the original equation.

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

CuO-CeO2 catalysts synthesized in one-step: Characterization and PROX performance

PEM fuel cells seem to be the most affordable and commercially viable hydrogen-based cells, the biggest challenge being to obtain CO-free H-2 (<100 ppm) as the fuel. In this study, the use of CuO-CeO2 catalysts in preferential oxidation of CO to obtain CO-free H-2 (PROX reaction) was investigated. Ce1-xCuxO2 catalysts, with x (mol%) = 0, 0.01, 0.03, 0.05 and 0.10, were synthesized in one-step by the polymeric precursor method, to obtain a very fine dispersion and strong metal-support interaction, to favor active copper species and a preference for the PROX reaction. The results obtained from catalyzed reactions and characterization of the catalysts by XRD, Rietveld refinement, BET surface area, UV-Vis and TPR, suggest that this one-step synthesis method gives rise to catalysts with copper species selective for the PROX reaction, which reaches a maximum rate on Ce0.97Cu0.03O2 catalyst. Copyright (C) 2012, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

A new relaxation method for roll forming problems

Finite element analysis (FEA) of nonlinear problems in solid mechanics is a time consuming process, but it can deal rigorously with the problems of both geometric, contact and material nonlinearity that occur in roll forming. The simulation time limits the application of nonlinear FEA to these problems in industrial practice, so that most applications of nonlinear FEA are in theoretical studies and engineering consulting or troubleshooting. Instead, quick methods based on a global assumption of the deformed shape have been used by the roll-forming industry. These approaches are of limited accuracy. This paper proposes a new form-finding method - a relaxation method to solve the nonlinear problem of predicting the deformed shape due to plastic deformation in roll forming. This method involves applying a small perturbation to each discrete node in order to update the local displacement field, while minimizing plastic work. This is iteratively applied to update the positions of all nodes. As the method assumes a local displacement field, the strain and stress components at each node are calculated explicitly. Continued perturbation of nodes leads to optimisation of the displacement field. Another important feature of this paper is a new approach to consideration of strain history. For a stable and continuous process such as rolling and roll forming, the strain history of a point is represented spatially by the states at a row of nodes leading in the direction of rolling to the current one. Therefore the increment of the strain components and the work-increment of a point can be found without moving the object forward. Using this method we can find the solution for rolling or roll forming in just one step. This method is expected to be faster than commercial finite element packages by eliminating repeated solution of large sets of simultaneous equations and the need to update boundary conditions that represent the rolls.

Laser ablation-based one-step generation and bio-functionalization of gold nanoparticles conjugated with aptamers

Background: Bio-conjugated nanoparticles are important analytical tools with emerging biological and medical applications. In this context, in situ conjugation of nanoparticles with biomolecules via laser ablation in an aqueous media is a highly promising one-step method for the production of functional nanoparticles resulting in highly efficient conjugation. Increased yields are required, particularly considering the conjugation of cost-intensive biomolecules like RNA aptamers. Results: Using a DNA aptamer directed against streptavidin, in situ conjugation results in nanoparticles with diameters of approximately 9 nm exhibiting a high aptamer surface density (98 aptamers per nanoparticle) and a maximal conjugation efficiency of 40.3%. We have demonstrated the functionality of the aptamer-conjugated nanoparticles using three independent analytical methods, including an agglomeration-based colorimetric assay, and solid-phase assays proving high aptamer activity. To demonstrate the general applicability of the in situ conjugation of gold nanoparticles with aptamers, we have transferred the method to an RNA aptamer directed against prostate-specific membrane antigen (PSMA). Successful detection of PSMA in human prostate cancer tissue was achieved utilizing tissue microarrays. Conclusions: In comparison to the conventional generation of bio-conjugated gold nanoparticles using chemical synthesis and subsequent bio-functionalization, the laser-ablation-based in situ conjugation is a rapid, one-step production method. Due to high conjugation efficiency and productivity, in situ conjugation can be easily used for high throughput generation of gold nanoparticles conjugated with valuable biomolecules like aptamers.

Dismantling the public sector bastion : evaluation capital works

Purpose: The determination of the most appropriate procurement system for a capital works project is a challenging task for public sector clients considering the array of assessment criteria that are considered and the procurement methods that are available. This is particularly pertinent to the Western Australian public sector where they have had a propensity to use traditional lump sum as the default procurement solution despite knowing that the selection of an inappropriate procurement method may lead to cost and time overruns, claims, and disputes’ on projects. This paper presents a six step procurement method evaluation approach that requires public sector agencies to consider in detail an array of options so as to obtain value for money. Design/methodology/approach: A procurement evaluation approach is developed and is examined using a focus group of 12 participants comprising of a public sector client, project team and key stakeholders. The focus group was used to examine the developed approach in the context of a real-life capital works project. Findings: The procurement method evaluation approach was deemed to be pragmatic and enabled decision-makers to re-evaluate outcomes from previous steps in the process. All focus group participants stated the six step process enabled a recommendation that was grounded in reflection and detailed evaluation. Practical implications: The developed procurement approach has enabled the public sector client evaluate the way in which they view procurement method selection and examine how they obtain ‘value for money’. Originality/value: The six step procurement approach makes use of quantitative and qualitative techniques and is reliant on discourse and reflection in making a procurement method recommendation. Consequently, the approach enables public sector clients to account for the complexities often associated with procurement selection.