Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

Kinetic energy sum spectra in nonmesonic weak decay of hypernuclei

We evaluate the coincidence spectra in the nonmesonic weak decay (NMWD) Lambda N -> nN of Lambda hypernuclei (4)(Lambda)He, (5)(Lambda)He, (12)(Lambda)C, (16)(Lambda)O, and (28)(Lambda)Si, as a function of the sum of kinetic energies E(nN)=E(n)+E(N) for N=n,p. The strangeness-changing transition potential is described by the one-meson-exchange model, with commonly used parametrization. Two versions of the independent-particle shell model (IPSM) are employed to account for the nuclear structure of the final residual nuclei. They are as follows: (a) IPSM-a, where no correlation, except for the Pauli principle, is taken into account and (b) IPSM-b, where the highly excited hole states are considered to be quasistationary and are described by Breit-Wigner distributions, whose widths are estimated from the experimental data. All np and nn spectra exhibit a series of peaks in the energy interval 110 MeV < E(nN)< 170 MeV, one for each occupied shell-model state. Within the IPSM-a, and because of the recoil effect, each peak covers an energy interval proportional to A(-1) , going from congruent to 4 MeV for (28)(Lambda)Si to congruent to 40 MeV for (4)(Lambda)He. Such a description could be pretty fair for the light (4)(Lambda)He and (5)(Lambda)He hypernuclei. For the remaining, heavier, hypernuclei it is very important, however, to consider as well the spreading in strength of the deep-hole states and bring into play the IPSM-b approach. Notwithstanding the nuclear model that is employed the results depend only very weakly on the details of the dynamics involved in the decay process proper. We propose that the IPSM is the appropriate lowest-order approximation for the theoretical calculations of the of kinetic energy sum spectra in the NMWD. It is in comparison to this picture that one should appraise the effects of the final-state interactions and of the two-nucleon-induced decay mode.

Behavior of Socket Base Connections Emphasizing Pedestal Walls

To evaluate the main design models for socket base connections of precast concrete structures, an experimental investigation was carried out on specimens of this connection with smooth and rough interfaces in contact with cast-in-place concrete. The specimens consisted of pedestal walls and were submitted to loads with large eccentricities. Based on the experimental results, two rational design models are proposed for this connection. One of these models accounts for the friction and is applied to socket bases with smooth interfaces. The main behavior model was verified for sockets with this type of interface and the design of the longitudinal walls as corbels is also suggested in this case. Because the behavior of the rough interface specimens was very close to a monolithic connection, the other proposed model is an adaptation of the bending theory to calculate the vertical reinforcement of socket bases with rough interfaces.

Robust integration of real time optimization with linear model predictive control

Here, we study the stable integration of real time optimization (RTO) with model predictive control (MPC) in a three layer structure. The intermediate layer is a quadratic programming whose objective is to compute reachable targets to the MPC layer that lie at the minimum distance to the optimum set points that are produced by the RTO layer. The lower layer is an infinite horizon MPC with guaranteed stability with additional constraints that force the feasibility and convergence of the target calculation layer. It is also considered the case in which there is polytopic uncertainty in the steady state model considered in the target calculation. The dynamic part of the MPC model is also considered unknown but it is assumed to be represented by one of the models of a discrete set of models. The efficiency of the methods presented here is illustrated with the simulation of a low order system. (C) 2010 Elsevier Ltd. All rights reserved.

Pion mass effects on axion emission from neutron stars through NN bremsstrahlung processes

The rates of axion emission by nucleon-nucleon bremsstrahlung are calculated with the inclusion of the full momentum contribution from a nuclear one pion exchange (OPE) potential. The contributions of the neutron-neutron (nn), proton-proton (pp) and neutron-proton (np) processes in both the non-degenerate and degenerate limits are explicitly given. We find that the finite-momentum corrections to the emissivities are quantitatively significant for the non-degenerate regime and temperature-dependent, and should affect the existing axion mass hounds. The trend of these nuclear effects is to diminish the emissivities. (C) 2009 Elsevier B.V. All rights reserved.

Weak nonmesonic decay spectra of hypernuclei

We compute one- and two-nucleon kinetic-energy spectra and opening-angle distributions for the nonmesonic weak decay of several hypernuclei, and compare our results with some recent data. The decay dynamics is described by transition potentials of the one-meson-exchange type, and the nuclear structure aspects by two versions of the independent-particle shell model (IPSM). In version IPSM-a, the bole states are treated as stationary, while in version IPSM-b the deep-hole ones are considered to be quasi-stationary and are described by Breit-Wigner distributions.

Nonmesonic weak decay spectra of (4)(Lambda)He

To comprehend the recent Brookhaven National Laboratory experiment E788 on (4)(Lambda)He, we have outlined a simple theoretical framework. based on the independent-particle shell model, for the one-nucleon-induced nonmesonic weak decay spectra. Basically, the shapes of all the spectra are tailored by the kinematics of the corresponding phase space, depending very weakly on the dynamics, which is gauged here by the one-meson-exchange potential. In spite of the straightforwardness of the approach a good agreement with data is achieved. This might be an indication that the final-state-interactions and the two-nucleon induced processes are not very important in the decay of this hypernucleus. We have also found that the pi + K exchange potential with soft vertex-form-factor cutoffs (Lambda(pi) approximate to 0.7 GeV, Lambda(K) approximate to 0.9 GeV), is able to account simultaneously for the available experimental data related to Gamma(p) and Gamma(n) for (4)(Lambda)H, (4)(Lambda)H, and (5)(Lambda)H. (C) 2009 Elsevier B.V. All rights reserved.

Effects of Micelles and Vesicles on the Oximolysis of p-Nitrophenyl Diphenyl Phosphate: A Model System for Surfactant-Based Skin-defensive Formulations against Organophosphates

The rates of oximolysis of p-nitrophenyl diphenyl phosphate (PNPDPP) by Acetophenoxime; 10-phenyl-10-hydi-oxyiminodecanoic acid; 4-(9-carboxynonanyl)-1-(9-carboxy-1-hydroyiminononanyl) benzene; 1-dodecyl-2-[(hydroxyimino)methyl]-pyridinium chloride (IV) and N-methylpyridinium-2-aldoxime chloride were determined in micelles of N-hexadecyl-N,N,N-trimethylammonium chloride (CTAC), N-hexadecyl-N,N-dimethylammonium propanesulfonate and dioctadecyldimethylammonium chloride (DODAC) vesicles. The effects of CTAC micelles and DODAC vesicles on the rates of oxymolysis of O,O-Diethyl O-(4-nitrophenyl) phosphate (paraoxon) by oxime IV were also determined. Analysis of micellar and vesicular effects on oximolysis of PNPDPP, using pseudophase or pseudophase with explicit consideration of ion exchange models, required the determination of the aggregate`s effects on the pK(a), of oximes and on the rates of PNPDPP hydrolysis. All aggregates increased the rate of oximolysis of PNPDPP and the results were analyzed quantitatively. In particular, DODAC vesicles catalyzed the reaction and increased the rate of oximolysis of PNPDPP by IV several million fold at pH`s compatible with pharmaceutical formulations. The rate increase produced by DODAC vesicles on the rate of oximolysis paraoxon by IV demonstrates the pharmaceutical potential of this system, since the substrate is used as an agricultural defensive agent and the surfactant is extensively employed in cosmetic formulations. (C) 2008 Wiley-Liss, Inc. and the American Pharmacists Association J Pharm Sci 98:1040-1052, 2009

Merging Resource Availability with Isotope Mixing Models: The Role of Neutral Interaction Assumptions

Background: Bayesian mixing models have allowed for the inclusion of uncertainty and prior information in the analysis of trophic interactions using stable isotopes. Formulating prior distributions is relatively straightforward when incorporating dietary data. However, the use of data that are related, but not directly proportional, to diet (such as prey availability data) is often problematic because such information is not necessarily predictive of diet, and the information required to build a reliable prior distribution for all prey species is often unavailable. Omitting prey availability data impacts the estimation of a predator's diet and introduces the strong assumption of consumer ultrageneralism (where all prey are consumed in equal proportions), particularly when multiple prey have similar isotope values. Methodology: We develop a procedure to incorporate prey availability data into Bayesian mixing models conditional on the similarity of isotope values between two prey. If a pair of prey have similar isotope values (resulting in highly uncertain mixing model results), our model increases the weight of availability data in estimating the contribution of prey to a predator's diet. We test the utility of this method in an intertidal community against independently measured feeding rates. Conclusions: Our results indicate that our weighting procedure increases the accuracy by which consumer diets can be inferred in situations where multiple prey have similar isotope values. This suggests that the exchange of formalism for predictive power is merited, particularly when the relationship between prey availability and a predator's diet cannot be assumed for all species in a system.

Black-box decomposition approach for computational hemodynamics: One-dimensional models

Increasing efforts exist in integrating different levels of detail in models of the cardiovascular system. For instance, one-dimensional representations are employed to model the systemic circulation. In this context, effective and black-box-type decomposition strategies for one-dimensional networks are needed, so as to: (i) employ domain decomposition strategies for large systemic models (1D-1D coupling) and (ii) provide the conceptual basis for dimensionally-heterogeneous representations (1D-3D coupling, among various possibilities). The strategy proposed in this article works for both of these two scenarios, though the several applications shown to illustrate its performance focus on the 1D-1D coupling case. A one-dimensional network is decomposed in such a way that each coupling point connects two (and not more) of the sub-networks. At each of the M connection points two unknowns are defined: the flow rate and pressure. These 2M unknowns are determined by 2M equations, since each sub-network provides one (non-linear) equation per coupling point. It is shown how to build the 2M x 2M non-linear system with arbitrary and independent choice of boundary conditions for each of the sub-networks. The idea is then to solve this non-linear system until convergence, which guarantees strong coupling of the complete network. In other words, if the non-linear solver converges at each time step, the solution coincides with what would be obtained by monolithically modeling the whole network. The decomposition thus imposes no stability restriction on the choice of the time step size. Effective iterative strategies for the non-linear system that preserve the black-box character of the decomposition are then explored. Several variants of matrix-free Broyden`s and Newton-GMRES algorithms are assessed as numerical solvers by comparing their performance on sub-critical wave propagation problems which range from academic test cases to realistic cardiovascular applications. A specific variant of Broyden`s algorithm is identified and recommended on the basis of its computer cost and reliability. (C) 2010 Elsevier B.V. All rights reserved.

Oriented Percolation in One-dimensional 1/vertical bar x-y vertical bar(2) Percolation Models

We consider independent edge percolation models on Z, with edge occupation probabilities. We prove that oriented percolation occurs when beta > 1 provided p is chosen sufficiently close to 1, answering a question posed in Newman and Schulman (Commun. Math. Phys. 104: 547, 1986). The proof is based on multi-scale analysis.

Soluble one-dimensional particle conservation models with infinitely many absorbing states

Particle conservation lattice-gas models with infinitely many absorbing states are studied on a one-dimensional lattice. As one increases the particle density, they exhibit a phase transition from an absorbing to an active phase. The models are solved exactly by the use of the transfer matrix technique from which the critical behavior was obtained. We have found that the exponent related to the order parameter, the density of active sites, is 1 for all studied models except one of them with exponent 2.

Comparative study of dental arch width in plaster models, photocopies and digitized images

The aim of this study was to comparatively assess dental arch width, in the canine and molar regions, by means of direct measurements from plaster models, photocopies and digitized images of the models. The sample consisted of 130 pairs of plaster models, photocopies and digitized images of the models of white patients (n = 65), both genders, with Class I and Class II Division 1 malocclusions, treated by standard Edgewise mechanics and extraction of the four first premolars. Maxillary and mandibular intercanine and intermolar widths were measured by a calibrated examiner, prior to and after orthodontic treatment, using the three modes of reproduction of the dental arches. Dispersion of the data relative to pre- and posttreatment intra-arch linear measurements (mm) was represented as box plots. The three measuring methods were compared by one-way ANOVA for repeated measurements (α = 0.05). Initial / final mean values varied as follows: 33.94 to 34.29 mm / 34.49 to 34.66 mm (maxillary intercanine width); 26.23 to 26.26 mm / 26.77 to 26.84 mm (mandibular intercanine width); 49.55 to 49.66 mm / 47.28 to 47.45 mm (maxillary intermolar width) and 43.28 to 43.41 mm / 40.29 to 40.46 mm (mandibular intermolar width). There were no statistically significant differences between mean dental arch widths estimated by the three studied methods, prior to and after orthodontic treatment. It may be concluded that photocopies and digitized images of the plaster models provided reliable reproductions of the dental arches for obtaining transversal intra-arch measurements.

Allometric models for estimating the phytomass of a secondary Atlantic Forest area of southeastern Brazil

The purpose of this study was to develop and validate equations to estimate the aboveground phytomass of a 30 years old plot of Atlantic Forest. In two plots of 100 m², a total of 82 trees were cut down at ground level. For each tree, height and diameter were measured. Leaves and woody material were separated in order to determine their fresh weights in field conditions. Samples of each fraction were oven dried at 80 °C to constant weight to determine their dry weight. Tree data were divided into two random samples. One sample was used for the development of the regression equations, and the other for validation. The models were developed using single linear regression analysis, where the dependent variable was the dry mass, and the independent variables were height (h), diameter (d) and d²h. The validation was carried out using Pearson correlation coefficient, paired t-Student test and standard error of estimation. The best equations to estimate aboveground phytomass were: lnDW = -3.068+2.522lnd (r² = 0.91; s y/x = 0.67) and lnDW = -3.676+0.951ln d²h (r² = 0.94; s y/x = 0.56).

Ventilation behavior in trained and untrained men during incremental test: evidence of one metabolic transition point

This study aimed to describe and compare the ventilation behavior during an incremental test utilizing three mathematical models and to compare the feature of ventilation curve fitted by the best mathematical model between aerobically trained (TR) and untrained ( UT) men. Thirty five subjects underwent a treadmill test with 1 km.h(-1) increases every minute until exhaustion. Ventilation averages of 20 seconds were plotted against time and fitted by: bi-segmental regression model (2SRM); three-segmental regression model (3SRM); and growth exponential model (GEM). Residual sum of squares (RSS) and mean square error (MSE) were calculated for each model. The correlations between peak VO2 (VO2PEAK), peak speed (Speed(PEAK)), ventilatory threshold identified by the best model (VT2SRM) and the first derivative calculated for workloads below (moderate intensity) and above (heavy intensity) VT2SRM were calculated. The RSS and MSE for GEM were significantly higher (p < 0.01) than for 2SRM and 3SRM in pooled data and in UT, but no significant difference was observed among the mathematical models in TR. In the pooled data, the first derivative of moderate intensities showed significant negative correlations with VT2SRM (r = -0.58; p < 0.01) and Speed(PEAK) (r = -0.46; p < 0.05) while the first derivative of heavy intensities showed significant negative correlation with VT2SRM (r = -0.43; p < 0.05). In UT group the first derivative of moderate intensities showed significant negative correlations with VT2SRM (r = -0.65; p < 0.05) and Speed(PEAK) (r = -0.61; p < 0.05), while the first derivative of heavy intensities showed significant negative correlation with VT2SRM (r= -0.73; p < 0.01), Speed(PEAK) (r = -0.73; p < 0.01) and VO2PEAK (r = -0.61; p < 0.05) in TR group. The ventilation behavior during incremental treadmill test tends to show only one threshold. UT subjects showed a slower ventilation increase during moderate intensities while TR subjects showed a slower ventilation increase during heavy intensities.