Numerical prediction of gas concentrations and fluctuations above a triangular hill within a turbulent boundary layer

The objective of the present work is to propose a numerical and statistical approach, using computational fluid dynamics, for the study of the atmospheric pollutant dispersion. Modifications in the standard k-epsilon turbulence model and additional equations for the calculation of the variance of concentration are introduced to enhance the prediction of the flow field and scalar quantities. The flow field, the mean concentration and the variance of a flow over a two-dimensional triangular hill, with a finite-size point pollutant source, are calculated by a finite volume code and compared with published experimental results. A modified low Reynolds k-epsilon turbulence model was employed in this work, using the constant of the k-epsilon model C(mu)=0.03 to take into account the inactive atmospheric turbulence. The numerical results for the velocity profiles and the position of the reattachment point are in good agreement with the experimental results. The results for the mean and the variance of the concentration are also in good agreement with experimental results from the literature. (C) 2009 Elsevier Ltd. All rights reserved.

Pitfalls driven by the sole use of local updates in dynamical

The recent claim that the exit probability (EP) of a slightly modified version of the Sznadj model is a continuous function of the initial magnetization is questioned. This result has been obtained analytically and confirmed by Monte Carlo simulations, simultaneously and independently by two different groups (EPL, 82 (2008) 18006; 18007). It stands at odds with an earlier result which yielded a step function for the EP (Europhys. Lett., 70 (2005) 705). The dispute is investigated by proving that the continuous shape of the EP is a direct outcome of a mean-field treatment for the analytical result. As such, it is most likely to be caused by finite-size effects in the simulations. The improbable alternative would be a signature of the irrelevance of fluctuations in this system. Indeed, evidence is provided in support of the stepwise shape as going beyond the mean-field level. These findings yield new insight in the physics of one-dimensional systems with respect to the validity of a true equilibrium state when using solely local update rules. The suitability and the significance to perform numerical simulations in those cases is discussed. To conclude, a great deal of caution is required when applying updates rules to describe any system especially social systems. Copyright (C) EPLA, 2011

CRITICAL BEHAVIOR AND THRESHOLD OF COEXISTENCE OF A PREDATOR-PREY STOCHASTIC MODEL IN A 2D LATTICE

We investigate the critical behavior of a stochastic lattice model describing a predator-prey system. By means of Monte Carlo procedure we simulate the model defined on a regular square lattice and determine the threshold of species coexistence, that is, the critical phase boundaries related to the transition between an active state, where both species coexist and an absorbing state where one of the species is extinct. A finite size scaling analysis is employed to determine the order parameter, order parameter fluctuations, correlation length and the critical exponents. Our numerical results for the critical exponents agree with those of the directed percolation universality class. We also check the validity of the hyperscaling relation and present the data collapse curves.

Extending the use of canonical and microcanonical Monte Carlo algorithms to spin models

We describe the canonical and microcanonical Monte Carlo algorithms for different systems that can be described by spin models. Sites of the lattice, chosen at random, interchange their spin values, provided they are different. The canonical ensemble is generated by performing exchanges according to the Metropolis prescription whereas in the microcanonical ensemble, exchanges are performed as long as the total energy remains constant. A systematic finite size analysis of intensive quantities and a comparison with results obtained from distinct ensembles are performed and the quality of results reveal that the present approach may be an useful tool for the study of phase transitions, specially first-order transitions. (C) 2009 Elsevier B.V. All rights reserved.

Shared information in stationary states at criticality

We consider bipartitions of one-dimensional extended systems whose probability distribution functions describe stationary states of stochastic models. We define estimators of the information shared between the two subsystems. If the correlation length is finite, the estimators stay finite for large system sizes. If the correlation length diverges, so do the estimators. The definition of the estimators is inspired by information theory. We look at several models and compare the behaviors of the estimators in the finite-size scaling limit. Analytical and numerical methods as well as Monte Carlo simulations are used. We show how the finite-size scaling functions change for various phase transitions, including the case where one has conformal invariance.

A conformal invariant growth model

We present a one-parameter extension of the raise and peel one-dimensional growth model. The model is defined in the configuration space of Dyck (RSOS) paths. Tiles from a rarefied gas hit the interface and change its shape. The adsorption rates are local but the desorption rates are non-local; they depend not only on the cluster hit by the tile but also on the total number of peaks (local maxima) belonging to all the clusters of the configuration. The domain of the parameter is determined by the condition that the rates are non-negative. In the finite-size scaling limit, the model is conformal invariant in the whole open domain. The parameter appears in the sound velocity only. At the boundary of the domain, the stationary state is an adsorbing state and conformal invariance is lost. The model allows us to check the universality of non-local observables in the raise and peel model. An example is given.

Phase diagram and spectral properties of a new exactly integrable spin-1 quantum chain

The spectral properties and phase diagram of the exactly integrable spin-1 quantum chain introduced by Alcaraz and Bariev are presented. The model has a U(1) symmetry and its integrability is associated with an unknown R-matrix whose dependence on the spectral parameters is not of a different form. The associated Bethe ansatz equations that fix the eigenspectra are distinct from those associated with other known integrable spin models. The model has a free parameter t(p). We show that at the special point t(p) = 1, the model acquires an extra U(1) symmetry and reduces to the deformed SU(3) Perk-Schultz model at a special value of its anisotropy q = exp(i2 pi/3) and in the presence of an external magnetic field. Our analysis is carried out either by solving the associated Bethe ansatz equations or by direct diagonalization of the quantum Hamiltonian for small lattice sizes. The phase diagram is calculated by exploring the consequences of conformal invariance on the finite-size corrections of the Hamiltonian eigenspectrum. The model exhibits a critical phase ruled by the c = 1 conformal field theory separated from a massive phase by first-order phase transitions.

Spin-density functional for exchange anisotropic Heisenberg model

Ground-state energies for anti ferromagnetic Heisenberg models with exchange anisotropy are estimated by means of a local-spin approximation made in the context of the density functional theory. Correlation energy is obtained using the non-linear spin-wave theory for homogeneous systems from which the spin functional is built. Although applicable to chains of any size, the results are shown for small number of sites, to exhibit finite-size effects and allow comparison with exact-numerical data from direct diagonalization of small chains. (C) 2009 Elsevier B.V. All rights reserved.

TEM, XRD and AFM study of poly(o-ethoxyaniline) films: new evidence for the formation of conducting islands

The existence of conducting islands in polyaniline films has long been proposed in the literature, which would be consistent with conducting mechanisms based on hopping. Obtaining direct evidence of conducting islands, however, is not straightforward. In this paper, conducting islands were visualized in poly(o-ethoxyaniline) (POEA) films prepared at low pH, using Transmission Electron Microscopy (TEM) and atomic force spectroscopy (AFS). The size of the islands varied between 67 and 470 angstrom for a pH=3.0, with a larger average being obtained with AFS, probably due to the finite size effect of the atomic force microscopy tip. In AFS, the conducting islands were denoted by regions with repulsive forces due to the double-layer forces. On the basis of X-ray diffraction (XRD) patterns for POEA in the powder form, we infer that the conducting islands are crystalline, and therefore a POEA film is believed to consist of conducting islands dispersed in an insulating, amorphous matrix. From conductivity measurements we inferred the charge transport to be governed by a typical quasi-one dimensional variable range hopping (VRH) mechanism.

Density-Profile Processes Describing Biological Signaling Networks: Almost Sure Convergence to Deterministic Trajectories

We introduce jump processes in R(k), called density-profile processes, to model biological signaling networks. Our modeling setup describes the macroscopic evolution of a finite-size spin-flip model with k types of spins with arbitrary number of internal states interacting through a non-reversible stochastic dynamics. We are mostly interested on the multi-dimensional empirical-magnetization vector in the thermodynamic limit, and prove that, within arbitrary finite time-intervals, its path converges almost surely to a deterministic trajectory determined by a first-order (non-linear) differential equation with explicit bounds on the distance between the stochastic and deterministic trajectories. As parameters of the spin-flip dynamics change, the associated dynamical system may go through bifurcations, associated to phase transitions in the statistical mechanical setting. We present a simple example of spin-flip stochastic model, associated to a synthetic biology model known as repressilator, which leads to a dynamical system with Hopf and pitchfork bifurcations. Depending on the parameter values, the magnetization random path can either converge to a unique stable fixed point, converge to one of a pair of stable fixed points, or asymptotically evolve close to a deterministic orbit in Rk. We also discuss a simple signaling pathway related to cancer research, called p53 module.

On the finite element modeling of fatigue crack growth in pressurized cylindrical shells

A methodology for the computational modeling of the fatigue crack growth in pressurized shell structures, based on the finite element method and concepts of Linear Elastic Fracture Mechanics, is presented. This methodology is based on that developed by Potyondy [Potyondy D, Wawrzynek PA, Ingraffea, AR. Discrete crack growth analysis methodology for through crack in pressurized fuselage structures. Int J Numer Methods Eng 1995;38:1633-1644], which consists of using four stress intensity factors, computed from the modified crack integral method, to predict the fatigue propagation life as well as the crack trajectory, which is computed as part of the numerical simulation. Some issues not presented in the study of Potyondy are investigated herein such as the influence of the crack increment size and the number of nodes per element (4 or 9 nodes) on the simulation results by means of a fatigue crack propagation simulation of a Boeing 737 airplane fuselage. The results of this simulation are compared with experimental results and those obtained by Potyondy [1]. (C) 2008 Elsevier Ltd. All rights reserved.

Coarse-grid higher order finite-difference time-domain algorithm with low dispersion errors

Higher order (2,4) FDTD schemes used for numerical solutions of Maxwell`s equations are focused on diminishing the truncation errors caused by the Taylor series expansion of the spatial derivatives. These schemes use a larger computational stencil, which generally makes use of the two constant coefficients, C-1 and C-2, for the four-point central-difference operators. In this paper we propose a novel way to diminish these truncation errors, in order to obtain more accurate numerical solutions of Maxwell`s equations. For such purpose, we present a method to individually optimize the pair of coefficients, C-1 and C-2, based on any desired grid size resolution and size of time step. Particularly, we are interested in using coarser grid discretizations to be able to simulate electrically large domains. The results of our optimization algorithm show a significant reduction in dispersion error and numerical anisotropy for all modeled grid size resolutions. Numerical simulations of free-space propagation verifies the very promising theoretical results. The model is also shown to perform well in more complex, realistic scenarios.

Size and power properties of some tests in the Birnbaum-Saunders regression model

The Birnbaum-Saunders distribution has been used quite effectively to model times to failure for materials subject to fatigue and for modeling lifetime data. In this paper we obtain asymptotic expansions, up to order n(-1/2) and under a sequence of Pitman alternatives, for the non-null distribution functions of the likelihood ratio, Wald, score and gradient test statistics in the Birnbaum-Saunders regression model. The asymptotic distributions of all four statistics are obtained for testing a subset of regression parameters and for testing the shape parameter. Monte Carlo simulation is presented in order to compare the finite-sample performance of these tests. We also present two empirical applications. (C) 2010 Elsevier B.V. All rights reserved.

Influence of cervical preflaring on determinationof apical file size in mandibular molars: SEM analysis

This study investigated the influence of cervical preflaring with different rotary instruments on determination of the initial apical file (IAF) in mesiobuccal roots of mandibular molars. Fifty human mandibular molars whose mesial roots presented two clearly separated apical foramens (mesiobuccal and mesiolingual) were used. After standard access opening and removal of pulp tissue, the working length (WL) was determined at 1 mm short of the root apex. Five groups (n=10) were formed at random, according to the type of instrument used for cervical preflaring. In group 1, the size of the IAF was determined without preflaring of the cervical and middle root canal thirds. In groups 2 to 5, preflaring was performed with Gates-Glidden drills, ProTaper instruments, EndoFlare instruments and LA Axxes burs, respectively. Canals were sized manually with K-files, starting with size 08 K-files, inserted passively up to the WL. File sizes were increased until a binding sensation was felt at the WL and the size of the file was recorded. The instrument corresponding to the IAF was fixed into the canal at the WL with methylcyanoacrylate. The teeth were then sectioned transversally 1 mm short of the apex, with the IAF in position. Cross-sections of the WL region were examined under scanning electron microscopy and the discrepancies between canal diameter and the diameter of IAF were calculated using the tool "rule" (FEG) of the microscope's proprietary software. The measurements (µm) were analyzed statistically by Kruskal-Wallis and Dunn's tests at 5% significance level. There were statistically significant differences among the groups (p<0.05). The non-flared group had the greatest discrepancy (125.30 ± 51.54) and differed significantly from all flared groups (p<0.05). Cervical preflaring with LA Axxess burs produced the least discrepancies (55.10 ± 48.31), followed by EndoFlare instruments (68.20 ± 42.44), Gattes Glidden drills (68.90 ± 42.46) and ProTaper files (77.40 ± 73.19). However, no significant differences (p>0.05) were found among the rotary instruments. In conclusion, cervical preflaring improved IAF fitting to the canals at the WL in mesiobuccal roots of maxillary first molars. The rotary instruments evaluated in this study did not differ from each other regarding the discrepancies produced between the IAF size and canal diameter at the WL.

Repair of critical-size defects with autogenous periosteum-derived cells combined with bovine anorganic apatite/collagen: an experimental study in rat calvaria

The aim of this study was to evaluate the bone repair using autogenous periosteum-derived cells (PDC) and bovine anorganic apatite and collagen (HA-COL). PDC from Wistar rats (n=10) were seeded on HA-COL discs and subjected to osteoinduction during 6 days. Critical-size defects in rat calvarias were treated with blood clot (G1), autogenous bone (G2), HA-COL (G3) and HA-COL combined with PDC (G4) (n=40), and then analyzed 1 and 3 months after surgeries. Radiographic analysis exhibited no significant temporal change. G1 and G2 had discrete new marginal bone, but the radiopacity of graft materials in G2, G3 and G4 impaired the detection of osteogenesis. At 3 months, histopathological analysis showed the presence of ossification islets in G1, which was more evident in G2, homogeneous new bone around HA-COL in G3 and heterogeneous new bone around HA-COL in G4 in addition to moderate presence of foreign body cells in G3 and G4. Histomorphometric analysis showed no change in the volume density of xenograft (p>0.05) and bone volume density in G2 was twice greater than in G1 and G4 after 3 months (p<0.05), but similar to G3. The PDC did not increase bone formation in vivo, although the biomaterial alone showed biocompatibility and osteoconduction capacity.