A numerical method for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field

A finite difference technique, based on a projection method, is developed for solving the dynamic three-dimensional Ericksen-Leslie equations for nematic liquid crystals subject to a strong magnetic field. The governing equations in this situation are derived using primitive variables and are solved using the ideas behind the GENSMAC methodology (Tome and McKee [32]; Tome et al. [34]). The resulting numerical technique is then validated by comparing the numerical solution against an analytic solution for steady three-dimensional flow between two-parallel plates subject to a strong magnetic field. The validated code is then employed to solve channel flow for which there is no analytic solution. (C) 2009 Elsevier B.V. All rights reserved.

A Floating-point Extended Kalman Filter Implementation for Autonomous Mobile Robots

Localization and Mapping are two of the most important capabilities for autonomous mobile robots and have been receiving considerable attention from the scientific computing community over the last 10 years. One of the most efficient methods to address these problems is based on the use of the Extended Kalman Filter (EKF). The EKF simultaneously estimates a model of the environment (map) and the position of the robot based on odometric and exteroceptive sensor information. As this algorithm demands a considerable amount of computation, it is usually executed on high end PCs coupled to the robot. In this work we present an FPGA-based architecture for the EKF algorithm that is capable of processing two-dimensional maps containing up to 1.8 k features at real time (14 Hz), a three-fold improvement over a Pentium M 1.6 GHz, and a 13-fold improvement over an ARM920T 200 MHz. The proposed architecture also consumes only 1.3% of the Pentium and 12.3% of the ARM energy per feature.

Numerical solution of the PTT constitutive equation for unsteady three-dimensional free surface flows

This work deals with the development of a numerical technique for simulating three-dimensional viscoelastic free surface flows using the PTT (Phan-Thien-Tanner) nonlinear constitutive equation. In particular, we are interested in flows possessing moving free surfaces. The equations describing the numerical technique are solved by the finite difference method on a staggered grid. The fluid is modelled by a Marker-and-Cell type method and an accurate representation of the fluid surface is employed. The full free surface stress conditions are considered. The PTT equation is solved by a high order method, which requires the calculation of the extra-stress tensor on the mesh contours. To validate the numerical technique developed in this work flow predictions for fully developed pipe flow are compared with an analytic solution from the literature. Then, results of complex free surface flows using the FIT equation such as the transient extrudate swell problem and a jet flowing onto a rigid plate are presented. An investigation of the effects of the parameters epsilon and xi on the extrudate swell and jet buckling problems is reported. (C) 2010 Elsevier B.V. All rights reserved.

Gamma Knife(A (R)) 3-D Dose Distribution Mapping by Magnetic Resonance Imaging

In medical processes where ionizing radiation is used, dose planning and dose delivery are the key elements to patient safety and treatment success, particularly, when the delivered dose in a single session of treatment can be an order of magnitude higher than the regular doses of radiotherapy. Therefore, the radiation dose should be well defined and precisely delivered to the target while minimizing radiation exposure to surrounding normal tissues [1]. Several methods have been proposed to obtain three-dimensional (3-D) dose distribution [2, 3]. In this paper, we propose an alternative method, which can be easily implemented in any stereotactic radiosurgery center with a magnetic resonance imaging (MRI) facility. A phantom with or without scattering centers filled with Fricke gel solution is irradiated with Gamma Knife(A (R)) system at a chosen spot. The phantom can be a replica of a human organ such as head, breast or any other organ. It can even be constructed from a real 3-D MR image of an organ of a patient using a computer-aided construction and irradiated at a specific region corresponding to the tumor position determined by MRI. The spin-lattice relaxation time T (1) of different parts of the irradiated phantom is determined by localized spectroscopy. The T (1)-weighted phantom images are used to correlate the image pixels intensity to the absorbed dose and consequently a 3-D dose distribution with a high resolution is obtained.

One-dimensional trapped atoms: critical coupling parameter and critical number of particles

In this paper we consider the case of a Bose gas in low dimension in order to illustrate the applicability of a method that allows us to construct analytical relations, valid for a broad range of coupling parameters, for a function which asymptotic expansions are known. The method is well suitable to investigate the problem of stability of a collection of Bose particles trapped in one- dimensional configuration for the case where the scattering length presents a negative value. The eigenvalues for this interacting quantum one-dimensional many particle system become negative when the interactions overcome the trapping energy and, in this case, the system becomes unstable. Here we calculate the critical coupling parameter and apply for the case of Lithium atoms obtaining the critical number of particles for the limit of stability.

Influence from the sulfonate group (RSO(3)(-)) on the conversion process and emission efficiency of poly(p-phenylene vinylene)

The discovery of an alternative route to convert poly(xylyliden tetrahydrothiophenium chloride) (PTHT) into poly(p-phenylene vinylene) (PPV) using dodecylbenzenesulfonate (DBS) has allowed the formation of ultrathin films with unprecedented control of architecture and emission properties. In this work, we show that this route may be performed with several sufonated compounds where RSO(3)(-) replaces the counter-ion (Cl(-)) of PTHT, some of which are even more efficient than DBS. Spin-coating films were produced from PTHT and azo-dye molecules, an azo-polymer and organic salts as counter-ions of PTHT. The effects of the thermal annealing step of PTHT/RSO(3)(-) films at 110 and 230 degrees C were monitored by measuring the absorption and emission spectra. The results indicate that the exchange of the counterion Cl(-) of PTHT by a linear long chain with RSO(3)(-) group is a general procedure to obtain PPV polymer at lower conversion temperature (ca. 110 degrees C) with significant increase in the emission efficiency, regardless of the chemical position and the number of sulfonate groups. With the enhanced emission caused by Congo Red and Tinopal as counter-ions, it is demonstrated that the new synthetic route is entirely generic, which may allow accurate control of conversion and emission properties. (C) 2010 Elsevier B.V. All rights reserved.

Spectroscopic Investigation on the Formation of a Dinuclear Me(2)SO-Ru(2)-Mercaptopyridine Bridged Complex: [Ru(2)Cl(3)(mu-pyS)(mu-dmso)(dmso)(4)]center dot 2H(2)O

A dinuclear ruthenium(II) complex double-bridged by an N-aromatic ligand 2-mercaptopyridine (2-pyridinethiol or 2-pyridyl mercaptan) and a methyl sulfoxide (dmso) have been characterized by X-ray crystallography. The reported compound with formula [Ru(2)Cl(3) (mu-pyS)(mu-dmso)(dmso)(4)] center dot 2H(2)O, [C(15)H(36)Cl(3)NO(7)S(6)Ru(2)] (P2/c, a = 13.8175(2) angstrom, b = 10.5608(2) angstrom, c = 21.3544 (3) angstrom, beta = 106.090(1)degrees, V = 2,994.05(8) angstrom(3), Z = 4) represents a seven-membered ring system with both rutheniums in an octahedral geometry. All the hydrogen bonds (C-H-Cl) and the van der Waals contacts give rise to three-dimensional network in the structure and add stability to the dinuclear compound. To our knowledge, this is the first time that the formation of a dinuclear ruthenium(II) complex double-bridged by an N-aromatic ligand 2-mercaptopyridine and dmso have been reported. The study also provided valuable insight into bioinorganic chemistry as continuing efforts are being made to develop metal-based cancer chemotherapeutics. A major feature of this paper is the resolution of a double bridged ruthenium structure which contributes to a better understanding of ruthenium reactivity.

N-H...S=C hydrogen bond in O-alkyl N-methoxycarbonyl thiocarbamates, ROC(S)N(H)C(O)O)CH(3) (R = CH(3)(-), CH(3)CH(2)-)

Pure O-methyl N-methoxycarbonyl thiocarbamate CH(3)OC(S)N(H)C(O)OCH(3) (I) and O-ethyl N-methoxycarbonyl thiocarbamate, CH(3)CH(2)OC(S)N(H)C(O)OCH(3) (II), are quantitatively prepared by the addition reaction between the CH(3)OC(O)NCS and the corresponding alcohols. The compounds are characterized by multinuclear ((1)H and (13)C) and bi-dimensional ((13)C HSQC) NMR, GC-MS and FTIR spectroscopy techniques. Structural and conformational properties are analyzed using a combined approach involving crystallographic data, vibration spectra and theoretical calculations. The low-temperature (150 K) crystal structure of II was determined by X-ray diffraction methods. The substance crystallizes in the monoclinic space group P2(1)/n with a = 4.088(1)angstrom. b = 22.346(1)angstrom, c = 8.284(1)angstrom, beta = 100.687(3)degrees and Z = 4 molecules per unit cell. The conformation adopted by the thiocarbamate group -OC(S)N(H)- is syn (C=S double bond in synperiplanar orientation with respect to the N-H single bond), while the methoxycarbonyl C=O double bond is in antiperiplanar orientation with respect to the N-H bond. The non-H atoms in II are essentially coplanar and the molecules are arranged in the crystal lattice as centro-symmetric dimeric units held by N-H center dot center dot center dot S=C hydrogen bonds Id(N center dot center dot center dot S) = 3.387(1)angstrom, <(N-H center dot center dot center dot S) = 166.4(2)degrees]. Furthermore, the effect of the it electronic resonance in the structural and vibrational properties is also discussed. (C) 2009 Elsevier Ltd. All rights reserved.

Classical and hologram QSAR studies on a series of inhibitors of trypanosomatid Glyceraldehyde-3-Phosphate Dehydrogenase

Leishmaniasis and trypanosomiasis are major causes of morbidity and mortality in both tropical and subtropical regions of the world. The current available drugs are limited, ineffective, and require long treatment regimens. Due to the high dependence of trypanosomatids on glycolysis as a source of energy, some glycolytic enzymes have been identified as attractive targets for drug design. In the present work, classical Two-Dimensional Quantitative Structure -Activity Relationships (2D QSAR) and Hologram QSAR (HQSAR) studies were performed on a series of adenosine derivatives as inhibitors of Leishmania mexicana Glyceraldehyde-3-Phosphate Dehydrogenase (LmGAPDH). Significant correlation coefficients (classical QSAR, r(2)=0.83 and q(2) =0.81; HQSAR, r(2)=0.91 and q(2) =0.86) were obtained for the 56 training set compounds, indicating the potential of the models for untested compounds. The models were then externally validated using a test set of 14 structurally related compounds and the predicted values were in good agreement with the experimental results (classical QSAR, r(pred)(2) = 0.94; HQSAR, r(pred)(2) = 0.92).

Mode of operation and low-resolution structure of a multi-domain and hyperthermophilic endo-beta-1,3-glucanase from Thermotoga petrophila

1,3-beta-Glucan depolymerizing enzymes have considerable biotechnological applications including biofuel production, feedstock-chemicals and pharmaceuticals. Here we describe a comprehensive functional characterization and low-resolution structure of a hyperthermophilic laminarinase from Thermotoga petrophila (TpLam). We determine TpLam enzymatic mode of operation, which specifically cleaves internal beta-1,3-glucosidic bonds. The enzyme most frequently attacks the bond between the 3rd and 4th residue from the non-reducing end, producing glucose, laminaribiose and laminaritriose as major products. Far-UV circular dichroism demonstrates that TpLam is formed mainly by beta structural elements, and the secondary structure is maintained after incubation at 90 degrees C. The structure resolved by small angle X-ray scattering, reveals a multi-domain structural architecture of a V-shape envelope with a catalytic domain flanked by two carbohydrate-binding modules. Crown Copyright (C) 2011 Published by Elsevier Inc. All rights reserved.

Algorithms for two-dimensional cutting stock and strip packing problems using dynamic programming and column generation

We investigate several two-dimensional guillotine cutting stock problems and their variants in which orthogonal rotations are allowed. We first present two dynamic programming based algorithms for the Rectangular Knapsack (RK) problem and its variants in which the patterns must be staged. The first algorithm solves the recurrence formula proposed by Beasley; the second algorithm - for staged patterns - also uses a recurrence formula. We show that if the items are not so small compared to the dimensions of the bin, then these algorithms require polynomial time. Using these algorithms we solved all instances of the RK problem found at the OR-LIBRARY, including one for which no optimal solution was known. We also consider the Two-dimensional Cutting Stock problem. We present a column generation based algorithm for this problem that uses the first algorithm above mentioned to generate the columns. We propose two strategies to tackle the residual instances. We also investigate a variant of this problem where the bins have different sizes. At last, we study the Two-dimensional Strip Packing problem. We also present a column generation based algorithm for this problem that uses the second algorithm above mentioned where staged patterns are imposed. In this case we solve instances for two-, three- and four-staged patterns. We report on some computational experiments with the various algorithms we propose in this paper. The results indicate that these algorithms seem to be suitable for solving real-world instances. We give a detailed description (a pseudo-code) of all the algorithms presented here, so that the reader may easily implement these algorithms. (c) 2007 Elsevier B.V. All rights reserved.

STABILITY ANALYSIS WITH APPLICATIONS OF A TWO-DIMENSIONAL DYNAMICAL SYSTEM ARISING FROM A STOCHASTIC MODEL FOR AN ASSET MARKET

We analyze the stability properties of equilibrium solutions and periodicity of orbits in a two-dimensional dynamical system whose orbits mimic the evolution of the price of an asset and the excess demand for that asset. The construction of the system is grounded upon a heterogeneous interacting agent model for a single risky asset market. An advantage of this construction procedure is that the resulting dynamical system becomes a macroscopic market model which mirrors the market quantities and qualities that would typically be taken into account solely at the microscopic level of modeling. The system`s parameters correspond to: (a) the proportion of speculators in a market; (b) the traders` speculative trend; (c) the degree of heterogeneity of idiosyncratic evaluations of the market agents with respect to the asset`s fundamental value; and (d) the strength of the feedback of the population excess demand on the asset price update increment. This correspondence allows us to employ our results in order to infer plausible causes for the emergence of price and demand fluctuations in a real asset market. The employment of dynamical systems for studying evolution of stochastic models of socio-economic phenomena is quite usual in the area of heterogeneous interacting agent models. However, in the vast majority of the cases present in the literature, these dynamical systems are one-dimensional. Our work is among the few in the area that construct and study analytically a two-dimensional dynamical system and apply it for explanation of socio-economic phenomena.

A DISCRETIZATION SCHEME FOR AN ONE-DIMENSIONAL REACTION-DIFFUSION EQUATION WITH DELAY AND ITS DYNAMICS

The goal of this paper is to present an approximation scheme for a reaction-diffusion equation with finite delay, which has been used as a model to study the evolution of a population with density distribution u, in such a way that the resulting finite dimensional ordinary differential system contains the same asymptotic dynamics as the reaction-diffusion equation.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

COMMUTATIVE FINITE-DIMENSIONAL ALGEBRAS SATISFYING x(x(xy))=0 ARE NILPOTENT

We investigate the structure of commutative non-associative algebras satisfying the identity x(x(xy)) = 0. Recently, Correa and Hentzel proved that every commutative algebra satisfying above identity over a field of characteristic not equal 2 is solvable. We prove that every commutative finite-dimensional algebra u over a field F of characteristic not equal 2, 3 which satisfies the identity x(x(xy)) = 0 is nilpotent. Furthermore, we obtain new identities and properties for this class of algebras.