On the skew-normal calibration model

In this article, we present the EM-algorithm for performing maximum likelihood estimation of an asymmetric linear calibration model with the assumption of skew-normally distributed error. A simulation study is conducted for evaluating the performance of the calibration estimator with interpolation and extrapolation situations. As one application in a real data set, we fitted the model studied in a dimensional measurement method used for calculating the testicular volume through a caliper and its calibration by using ultrasonography as the standard method. By applying this methodology, we do not need to transform the variables to have symmetrical errors. Another interesting aspect of the approach is that the developed transformation to make the information matrix nonsingular, when the skewness parameter is near zero, leaves the parameter of interest unchanged. Model fitting is implemented and the best choice between the usual calibration model and the model proposed in this article was evaluated by developing the Akaike information criterion, Schwarz`s Bayesian information criterion and Hannan-Quinn criterion.

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.

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.

Two- and three-dimensional shape fabric analysis by the intercept method in grey levels

The count intercept is a robust method for the numerical analysis of fabrics Launeau and Robin (1996). It counts the number of intersections between a set of parallel scan lines and a mineral phase, which must be identified on a digital image. However, the method is only sensitive to boundaries and therefore supposes the user has some knowledge about their significance. The aim of this paper is to show that a proper grey level detection of boundaries along scan lines is sufficient to calculate the two-dimensional anisotropy of grain or crystal distributions without any particular image processing. Populations of grains and crystals usually display elliptical anisotropies in rocks. When confirmed by the intercept analysis, a combination of a minimum of 3 mean length intercept roses, taken on 3 more or less perpendicular sections, allows the calculation of 3-dimensional ellipsoids and the determination of their standard deviation with direction and intensity in 3 dimensions as well. The feasibility of this quick method is attested by numerous examples on theoretical objects deformed by active and passive deformation, on BSE images of synthetic magma flow, on drawing or direct analysis of thin section pictures of sandstones and on digital images of granites directly taken and measured in the field. (C) 2010 Elsevier B.V. All rights reserved.

ENO adaptive method for solving one-dimensional conservation laws

In this work an efficient third order non-linear finite difference scheme for solving adaptively hyperbolic systems of one-dimensional conservation laws is developed. The method is based oil applying to the solution of the differential equation an interpolating wavelet transform at each time step, generating a multilevel representation for the solution, which is thresholded and a sparse point representation is generated. The numerical fluxes obtained by a Lax-Friedrichs flux splitting are evaluated oil the sparse grid by an essentially non-oscillatory (ENO) approximation, which chooses the locally smoothest stencil among all the possibilities for each point of the sparse grid. The time evolution of the differential operator is done on this sparse representation by a total variation diminishing (TVD) Runge-Kutta method. Four classical examples of initial value problems for the Euler equations of gas dynamics are accurately solved and their sparse solutions are analyzed with respect to the threshold parameters, confirming the efficiency of the wavelet transform as an adaptive grid generation technique. (C) 2008 IMACS. Published by Elsevier B.V. All rights reserved.

Use of the DNA Checkerboard hybridization method for detection and quantitation of Candida species in oral microbiota

The DNA Checkerboard method enables the simultaneous identification of distinct microorganisms in a large number of samples and employs up to 45 whole genomic DNA probes to gram-negative and gram-positive bacterial species present in subgingival biofilms. Collectively, they account for 55%-60% of the bacteria in subgingival biofilms. In this study, we present the DNA Checkerboard hybridization as an alternative method for the detection and quantitation of Candida species in oral cavities. Our results reveal that DNA Checkerboard is sensitive enough and constitutes a powerful and appropriate method for detecting and quantifying Candida species found in the oral cavity.

Validation of an ELISA method for the serological diagnosis of canine brucellosis due to Brucella canis

In the present study, the validation of an enzyme-linked immunosorbent assay (ELISA) for serodiagnosis of canine brucellosis is described. Two different antigenic extracts, obtained by heat or ultrasonic homogenization of microbial antigens from a wild isolate of Brucella canis bacteria, were compared by ELISA and Western blot (WB). A total of 145 canine sera were used to define sensitivity, specificity and accuracy of the ELISA as follows: (1) sera from 34 animals with natural B. canis infection, confirmed by blood culture and PCR, as well as 51 sera samples from healthy dogs with negative results by the agar-gel immunodiffusion (ACID) test for canine brucellosis, were used as the control panel for B. cants infection; and (2) to scrutinize the possibility of cross reactions with other common dog infections in the same geographical area in Brazil, 60 sera samples from dogs harboring known infections by Leptospira sp., Ehrlichia canis, canine distemper virus (CDV), Neospora caninum, Babesia canis and Leishmania chagasi (10 in each group) were included in the study. The ELISA using heat soluble bacterial extract (HE-antigen) as antigen showed the best values of sensitivity (91.18%), specificity (100%) and accuracy (96.47%). In the WB analyses, the HE-antigen showed no cross-reactivity with sera from dogs with different infections, while the B. canis sonicate had various protein bands identified by those sera. The performance of the ELISA standardized with the heat soluble B. canis antigen indicates that this assay can be used as a reliable and practical method to confirm infection by this microorganism, as well as a tool for seroepidemiological studies. (C) 2010 Elsevier Ltd. All rights reserved.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

SEM and Neurohistological Observations of Nerve Endings in the Middle Region of the Tongue of the Collared Peccary (Tayassu tajacu): A Silver Impregnation Method

The presence of lingual papillae and the nerve endings in the middle region of the tongue mucosa of collared peccary (Tayassu tajacu) were studied using scanning electron microscopy and light microscopy, based upon the silver impregnation method. The middle region of tongue mucosa revealed numerous filiform and fungiform papillae. The thick epithelial layer showed epithelial cells and a dense connective tissue layer containing nerve fibre bundles and capillaries. The sensory nerve endings, intensely stained by silver impregnation, were usually non-encapsulated and extended into the connective tissue of the filiform and fungiform papillae very close to the epithelial cells. In some regions, the sensory nerves fibres formed a dense and complex network of fine fibrils. The presence of these nerve fibrils may characterize the mechanisms of transmission of sensitive impulses to the tongue mucosa.

The geodiamolide H, derived from Brazilian sponge Geodia corticostylifera, regulates actin cytoskeleton, migration and invasion of breast cancer cells cultured in three-dimensional environment

We are investigating effects of the depsipeptide geodiamolide H, isolated from the Brazilian sponge Geodia corticostylifera, on cancer cell lines grown in 3D environment. As shown previously geodiamolide H disrupts actin cytoskeleton in both sea urchin eggs and breast cancer cell monolayers. We used a normal mammary epithelial cell line MCF 10A that in 3D assay results formation of polarized spheroids. We also used cell lines derived from breast tumors with different degrees of differentiation: MCF7 positive for estrogen receptor and the Hs578T, negative for hormone receptors. Cells were placed on top of Matrigel. Spheroids obtained from these cultures were treated with geodiamolide H. Control and treated samples were analyzed by light and confocal microscopy. Geodiamolide H dramatically affected the poorly differentiated and aggressive Hs578T cell line. The peptide reverted HsS78T malignant phenotype to polarized spheroid-like structures. MCF7 cells treated by geodiamolide H exhibited polarization compared to controls. Geodiamolide H induced striking phenotypic modifications in Hs578T cell line and disruption of actin cytoskeleton. We investigated effects of geodiamolide H on migration and invasion of Hs578T cells. Time-lapse microscopy showed that the peptide inhibited migration of these cells in a dose-dependent manner. Furthermore invasion assays revealed that geodiamolide H induced a 30% decrease on invasive behavior of Hs578T cells. Our results suggest that geodiamolide H inhibits migration and invasion of Hs578T cells probably through modifications in actin cytoskeleton. The fact that normal cell lines were not affected by treatment with geodiamolide H stimulates new studies towards therapeutic use for this peptide.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.

A discontinuous-Galerkin-based immersed boundary method with non-homogeneous boundary conditions and its application to elasticity

We propose a discontinuous-Galerkin-based immersed boundary method for elasticity problems. The resulting numerical scheme does not require boundary fitting meshes and avoids boundary locking by switching the elements intersected by the boundary to a discontinuous Galerkin approximation. Special emphasis is placed on the construction of a method that retains an optimal convergence rate in the presence of non-homogeneous essential and natural boundary conditions. The role of each one of the approximations introduced is illustrated by analyzing an analog problem in one spatial dimension. Finally, extensive two- and three-dimensional numerical experiments on linear and nonlinear elasticity problems verify that the proposed method leads to optimal convergence rates under combinations of essential and natural boundary conditions. (C) 2009 Elsevier B.V. All rights reserved.

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.

Numerical Solution of the Upper-Convected Maxwell Model for Three-Dimensional Free Surface Flows

This work presents a finite difference technique for simulating three-dimensional free surface flows governed by the Upper-Convected Maxwell (UCM) constitutive equation. A Marker-and-Cell approach is employed to represent the fluid free surface and formulations for calculating the non-Newtonian stress tensor on solid boundaries are developed. The complete free surface stress conditions are employed. The momentum equation is solved by an implicit technique while the UCM constitutive equation is integrated by the explicit Euler method. The resulting equations are solved by the finite difference method on a 3D-staggered grid. By using an exact solution for fully developed flow inside a pipe, validation and convergence results are provided. Numerical results include the simulation of the transient extrudate swell and the comparison between jet buckling of UCM and Newtonian fluids.