Aging and stationary properties of non-equilibrium symmetrical three-state models

We consider a non-equilibrium three-state model whose dynamics is Markovian and displays the same symmetry as the three-state Potts model, i.e. the transition rates are invariant under the cyclic permutation of the states. Unlike the Potts model, detailed balance is, in general, not satisfied. The aging and the stationary properties of the model defined on a square lattice are obtained by means of large-scale Monte Carlo simulations. We show that the phase diagram presents a critical line, belonging to the three-state Potts universality class, that ends at a point whose universality class is that of the Voter model. Aging is considered on the critical line, at the Voter point and in the ferromagnetic phase.

A non-parametric vessel detection method for complex vascular structures

Modern medical imaging techniques enable the acquisition of in vivo high resolution images of the vascular system. Most common methods for the detection of vessels in these images, such as multiscale Hessian-based operators and matched filters, rely on the assumption that at each voxel there is a single cylinder. Such an assumption is clearly violated at the multitude of branching points that are easily observed in all, but the Most focused vascular image studies. In this paper, we propose a novel method for detecting vessels in medical images that relaxes this single cylinder assumption. We directly exploit local neighborhood intensities and extract characteristics of the local intensity profile (in a spherical polar coordinate system) which we term as the polar neighborhood intensity profile. We present a new method to capture the common properties shared by polar neighborhood intensity profiles for all the types of vascular points belonging to the vascular system. The new method enables us to detect vessels even near complex extreme points, including branching points. Our method demonstrates improved performance over standard methods on both 2D synthetic images and 3D animal and clinical vascular images, particularly close to vessel branching regions. (C) 2008 Elsevier B.V. All rights reserved.

Computational study of the step size parameter of the subgradient optimization method

The subgradient optimization method is a simple and flexible linear programming iterative algorithm. It is much simpler than Newton's method and can be applied to a wider variety of problems. It also converges when the objective function is non-differentiable. Since an efficient algorithm will not only produce a good solution but also take less computing time, we always prefer a simpler algorithm with high quality. In this study a series of step size parameters in the subgradient equation is studied. The performance is compared for a general piecewise function and a specific p-median problem. We examine how the quality of solution changes by setting five forms of step size parameter.

A contractive method for computing the stationary solution of the euler equation

A contractive method for computing stationary solutions of intertemporal equilibrium models is provide. The method is is implemented using a contraction mapping derived from the first-order conditions. The deterministic dynamic programming problem is used to illustrate the method. Some numerical examples are performed.

A non-invasive method for silencing gene transcription in honeybees maintained under natural conditions

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

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Non-destructive genetic sampling in fish. An improved method for DNA extraction from fish fins and scales

DNA-based studies have been one of the major interests in conservation biology of endangered species and in population genetics. As species and population genetic assessment requires a source of biological material, the sampling strategy can be overcome by non-destructive procedures for DNA isolation. An improved method for obtaining DNA from fish fins and scales with the use of an extraction buffer containing urea and further DNA purification with phenol-chloroform is described. The methodology combines the benefits of a non-destructive DNA sampling and its high efficiency. In addition, comparisons with other methodologies for isolating DNA from fish demonstrated that the present procedure also becomes a very attractive alternative to obtain large amounts of high-quality DNA for use in different molecular analyses. The DNA samples, isolated from different fish species, have been successfully used on random amplified polymorphic DNA (RAPD) experiments, as well as on amplification of specific ribosomal and mitochondrial DNA sequences. The present DNA extraction procedure represents an alternative for population approaches and genetic studies on rare or endangered taxa.

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.

Study of the thermal behavior of the transition phase of Co(II)-diclofenac compound by non-isothermal method

The Co(II)-diclofenac complex was evaluated by simultaneous thermogravimetry-differential thermal analysis (TG-DTA) and differential scanning calorimetry (DSC). The DTA curve profile shows one exothermic peak because of the transition phase of the compound between 170 and 180 A degrees C, which was confirmed by X-ray powder diffractometry. The transition phase behavior was studied by DSC curves at several heating rates of a sample mass between 1 and 10 mg in nitrogen atmosphere and in a crucible with and without a lid. Thus, the kinetic parameters were evaluated using an isoconversional non-linear fitting proposed by Capela and Ribeiro. The results show that the activation energy and pre-exponential factor for the transition phase is dependant on the different experimental conditions. Nevertheless, these results indicate that the kinetic compensation effect shows a relationship between them.

ACCELERATING ITERATIVE ALGORITHMS FOR SYMMETRICAL LINEAR COMPLEMENTARITY-PROBLEMS

We propose a method for accelerating iterative algorithms for solving symmetric linear complementarity problems. The method consists in performing a one-dimensional optimization in the direction generated by a splitting method even for non-descent directions. We give strong convergence proofs and present numerical experiments that justify using this acceleration.

Non-linear analysis of reinforced concrete plates by the boundary element method and continuum damage mechanics

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

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Investigations into the use of the finite element method and artificial neural networks in the non-destructive analysis of metallic tubes

This work presents an investigation into the use of the finite element method and artificial neural networks in the identification of defects in industrial plants metallic tubes, due to the aggressive actions of the fluids contained by them, and/or atmospheric agents. The methodology used in this study consists of simulating a very large number of defects in a metallic tube, using the finite element method. Both variations in width and height of the defects are considered. Then, the obtained results are used to generate a set of vectors for the training of a perceptron multilayer artificial neural network. Finally, the obtained neural network is used to classify a group of new defects, simulated by the finite element method, but that do not belong to the original dataset. The reached results demonstrate the efficiency of the proposed approach, and encourage future works on this subject.

On non-ideal and non-linear portal frame dynamics analysis using Bogoliubov averaging method

We apply the Bogoliubov Averaging Method to the study of the vibrations of an elastic foundation, forced by a Non-ideal energy source. The considered model consists of a portal plane frame with quadratic nonlinearities, with internal resonance 1:2, supporting a direct current motor with limited power. The non-ideal excitation is in primary resonance in the order of one-half with the second mode frequency. The results of the averaging method, plotted in time evolution curve and phase diagrams are compared to those obtained by numerically integrating of the original differential equations. The presence of the saturation phenomenon is verified by analytical procedures.

Electroosmotic pumping in rectangular microchannels: A numerical treatment by the finite element method

In this work, the analysis of electroosmotic pumping mechanisms in microchannels is performed through the solution of Poisson-Boltzmann and Navier Stokes equations by the Finite Element Method. This approach is combined with a Newton-Raphson iterative scheme, allowing a full treatment of the non-linear Poisson-Boltzmann source term which is normally approximated by linearizations in other methods.