Direct Computation of Operational Matrices for Polynomial Bases

Several numerical methods for boundary value problems use integral and differential operational matrices, expressed in polynomial bases in a Hilbert space of functions. This work presents a sequence of matrix operations allowing a direct computation of operational matrices for polynomial bases, orthogonal or not, starting with any previously known reference matrix. Furthermore, it shows how to obtain the reference matrix for a chosen polynomial base. The results presented here can be applied not only for integration and differentiation, but also for any linear operation.

Antecedentes do relacionamento e da performance em empresas da cadeia de suprimentos: estruturação e aplicação de modelos de equações estruturais

The relationship between companies is an important issue in the management of supply chains. Several aspects relating to the flow and exchange of information along the chain are considered as having a decisive influence on the success of this relationship. The main objective of this work was to structured and test models that link aspects of this nature with performance and the purchaser-supplier relationship in the supply chain. Aspects relevant to communication and the use do IT in relationships between companies were investigated. The importance of performance in this relationship was also investigated. The research were based on empirical data obtained by means of structural equation modeling. The results show that some aspects contribute in a significant way to the success of this relationship while others that, a priori, are considered important make no contribution.

Existence Results for Abstract Partial Neutral Integro-differential Equation with Unbounded Delay

In this paper we study the existence and regularity of mild solutions for a class of abstract partial neutral integro-differential equations with unbounded delay.

Shallow landslide prediction in the Serra do Mar, Sao Paulo, Brazil

Various methods are currently used in order to predict shallow landslides within the catchment scale. Among them, physically based models present advantages associated with the physical description of processes by means of mathematical equations. The main objective of this research is the prediction of shallow landslides using TRIGRS model, in a pilot catchment located at Serra do Mar mountain range, Sao Paulo State, southeastern Brazil. Susceptibility scenarios have been simulated taking into account different mechanical and hydrological values. These scenarios were analysed based on a landslide scars map from the January 1985 event, upon which two indexes were applied: Scars Concentration (SC - ratio between the number of cells with scars, in each class, and the total number of cells with scars within the catchment) and Landslide Potential (LP - ratio between the number of cells with scars, in each class, and the total number of cells in that same class). The results showed a significant agreement between the simulated scenarios and the scar's map. In unstable areas (SF <= 1), the SC values exceeded 50% in all scenarios. Based on the results, the use of this model should be considered an important tool for shallow landslide prediction, especially in areas where mechanical and hydrological properties of the materials are not well known.

Estimating cardiovascular risk in patients with type 2 diabetes: a national multicenter study in Brazil

According to Brazilian National Data Survey diabetes is the fifth cause for hospitalization and is one of the ten major causes of mortality in this country. Aims to stratify the estimated cardiovascular risk (eCVR) in a population of type 2 diabetics (T2DM) according to the Framingham prediction equations as well as to determine the association between eCVR with metabolic and clinical control of the disease. Methods From 2000 to 2001 a cross-sectional multicenter study was conducted in 13 public out-patients diabetes/endocrinology clinics from 8 Brazilian cities. The 10-year risk of developing coronary heart disease (CHD) was estimated by the prediction equations described by Wilson et al (Circulation 1998). LDL equations were preferably used; when patients missed LDL data we used total cholesterol equations instead. Results Data from 1382 patients (59.0% female) were analyzed. Median and inter-quartile range (IQ) of age and duration of diabetes were 57.4 (51-65) and 8.8 (3-13) years, respectively without differences according to the gender. Forty-two percent of these patients were overweight and 35.4% were obese (the prevalence of higher BMI and obesity in this T2DM group was significantly higher in women than in men; p < 0.001). The overall estimated eCVR in T2DM patients was 21.4 (13.5-31.3). The eCVR was high (> 20%) in 738 (53.4%), intermediate in 202 (14.6%) and low in 442 (32%) patients. Men [25.1(15.4-37.3)] showed a higher eCVR than women [18.8 (12.4-27.9) p < 0.001]. The most common risk factor was high LDL-cholesterol (80.8%), most frequently found in women than in men (p = 0.01). The median of risk factors present was three (2-4) without gender differences. Overall we observed that 60 (4.3%) of our patients had none, 154(11.1%) one, 310 (22.4%) two, 385 (27.9%) three, 300 (21.7%) four, 149 (10.5%) five and six, (2%) six risk factors. A higher eCVR was noted in overweight or obese patients (p = 0.01 for both groups). No association was found between eCVR with age or a specific type of diabetes treatment. A correlation was found between eCVR and duration of diabetes (p < 0.001), BMI (p < 0.001), creatinine (p < 0.001) and triglycerides levels (p < 0.001) but it was not found with HbA1c, fasting blood glucose and postprandial glucose. A higher eCVR was observed in patients with retinopathy (p < 0.001) and a tendency in patients with microalbuminuria (p = 0.06). Conclusion: our study showed that in this group of Brazilian T2DM the eCVR was correlated with the lipid profile and it was higher in patients with microvascular chronic complications. No correlation was found with glycemic control parameters. These data could explain the failure of intensive glycemic control programs aiming to reduce cardiovascular events observed in some studies.

Condylar displacement between centric relation and maximum intercuspation in symptomatic and asymptomatic individuals

Objective: To measure condylar displacement between centric relation (CR) and maximum intercuspation (MIC) in symptomatic and asymptomatic subjects. Materials and Methods: The sample comprised 70 non-deprogrammed individuals, divided equally into two groups, one symptomatic and the other asymptomatic, grouped according to the research diagnostic criteria for temporomandibular disorders (RDC/TMD). Condylar displacement was measured in three dimensions with the condylar position indicator (CPI) device. Dahlberg's index, intraclass correlation coefficient, repeated measures analysis of variance, analysis of variance, and generalized estimating equations were used for statistical analysis. Results: A greater magnitude of difference was observed on the vertical plane on the left side in both symptomatic and asymptomatic individuals (P = .033). The symptomatic group presented higher measurements on the transverse plane (P = .015). The percentage of displacement in the mesial direction was significantly higher in the asymptomatic group than in the symptomatic one (P = .049). Both groups presented a significantly higher percentage of mesial direction on the right side than on the left (P = .036). The presence of bilateral condylar displacement (left and right sides) in an inferior and distal direction was significantly greater in symptomatic individuals (P = .012). However, no statistical difference was noted between genders. Conclusion: Statistically significant differences between CR and MIC were quantifiable at the condylar level in asymptomatic and symptomatic individuals. (Angle Orthod. 2010;80:835-842.)

Newtonian perturbations on models with matter creation

Creation of cold dark matter (CCDM) can macroscopically be described by a negative pressure, and, therefore, the mechanism is capable to accelerate the Universe, without the need of an additional dark energy component. In this framework, we discuss the evolution of perturbations by considering a Neo-Newtonian approach where, unlike in the standard Newtonian cosmology, the fluid pressure is taken into account even in the homogeneous and isotropic background equations (Lima, Zanchin, and Brandenberger, MNRAS 291, L1, 1997). The evolution of the density contrast is calculated in the linear approximation and compared to the one predicted by the Lambda CDM model. The difference between the CCDM and Lambda CDM predictions at the perturbative level is quantified by using three different statistical methods, namely: a simple chi(2)-analysis in the relevant space parameter, a Bayesian statistical inference, and, finally, a Kolmogorov-Smirnov test. We find that under certain circumstances, the CCDM scenario analyzed here predicts an overall dynamics (including Hubble flow and matter fluctuation field) which fully recovers that of the traditional cosmic concordance model. Our basic conclusion is that such a reduction of the dark sector provides a viable alternative description to the accelerating Lambda CDM cosmology.

Equation of state for the magnetic-color-flavor-locked phase and its implications for compact star models

Using the solutions of the gap equations of the magnetic-color-flavor-locked (MCFL) phase of paired quark matter in a magnetic field, and taking into consideration the separation between the longitudinal and transverse pressures due to the field-induced breaking of the spatial rotational symmetry, the equation of state of the MCFL phase is self-consistently determined. This result is then used to investigate the possibility of absolute stability, which turns out to require a field-dependent ""bag constant"" to hold. That is, only if the bag constant varies with the magnetic field, there exists a window in the magnetic field vs bag constant plane for absolute stability of strange matter. Implications for stellar models of magnetized (self-bound) strange stars and hybrid (MCFL core) stars are calculated and discussed.

Fatores gravimétricos no Brasil para um modelo anelástico de Terra

Geodetic observations are affected by the disturbing potential of the luni-solar tide. Among those observations, the value of g obtained from gravimetric survey needs correction by the gravimetric factor. This correction is derived from the Numbers of Love, which depend on the adopted model of Earth. Because of this, it is necessary to update the correction since the gravimetric factor widely used in Brazil as delta = 1.20 does not consider local rheological variations and they are latitude dependent. A discrepancy of about 1% between the observed tidal gravimetric factors d of the ""Trans World Tidal Gravity Profiles"" (TWTGP), related to Brussels fundamental station, and those obtained by recent observations reported by Freitas and Ducarme ( 1991). Experiments based on inertial force effects also reveal a variation of about 0.5% in the observed d. A same order of magnitude difference is obtained for an anelastic Earth model when compared with a viscous-elastic model and even when different frequencies of tidal perturbations are considered. In this paper regression models are presented for gravimetric factors for the lunar components O(1) and M(2) in Brazil. These models were obtained from observations performed at stations belonging to the Brazilian segment of the TWTGP.

LIFE CYCLE AND POPULATION STRUCTURE OF AEGLA PAULENSIS (DECAPODA: ANOMURA: AEGLIDAE)

We describe growth, longevity, sex ratio, reproductive period, and recruitment of Aegla paulensis from Jaragua Stale Park, Sao Paulo, Brazil (23 degrees 27'27.9 '' S; 46 degrees 45'32.3 '' W). The population was sampled monthly (September 2007 through August 2009) with the aid of traps. Over five thousand individuals were captured, sexed, measured (carapace length = CL) and inspected for reproductive traits (females only), and then released back to the sampling site. The pattern of the reproductive cycle was strongly seasonal (austral mid autumn through late winter), with a single recruitment pulse per year. The obtained von Bertalanffy growth equations were CL = 21.25[1-e(-0.041(t + 1.250))] and CL = 16.52[1-e(-0.049(t + 1.823))] for males and females, respectively. Males (mean CL +/- SD = 11.86 +/- 2.79 mm) attain larger sizes than females (mean CL +/- SD = 10.84 +/- 2.36 mm). Aegla paulensis reproduces twice during an estimated life span of 40.2 months for females and 33.9 months for males. Temporal variation of sex ratio showed a distinctive pattern characterized by a sequence of three distinct periods that repeated from one year to another, and which suggested that a behavioral component influence the proportion of sex in adult specimens sampled with traps during reproductive and non-reproductive periods.

Modeling Two-Oscillator Circadian Systems Entrained by Two Environmental Cycles

Several experimental studies have altered the phase relationship between photic and non-photic environmental, 24 h cycles (zeitgebers) in order to assess their role in the synchronization of circadian rhythms. To assist in the interpretation of the complex activity patterns that emerge from these ""conflicting zeitgeber'' protocols, we present computer simulations of coupled circadian oscillators forced by two independent zeitgebers. This circadian system configuration was first employed by Pittendrigh and Bruce (1959), to model their studies of the light and temperature entrainment of the eclosion oscillator in Drosophila. Whereas most of the recent experiments have restricted conflicting zeitgeber experiments to two experimental conditions, by comparing circadian oscillator phases under two distinct phase relationships between zeitgebers (usually 0 and 12 h), Pittendrigh and Bruce compared eclosion phase under 12 distinct phase relationships, spanning the 24 h interval. Our simulations using non-linear differential equations replicated complex non-linear phenomena, such as ""phase jumps'' and sudden switches in zeitgeber preferences, which had previously been difficult to interpret. Our simulations reveal that these phenomena generally arise when inter-oscillator coupling is high in relation to the zeitgeber strength. Manipulations in the structural symmetry of the model indicated that these results can be expected to apply to a wide range of system configurations. Finally, our studies recommend the use of the complete protocol employed by Pittendrigh and Bruce, because different system configurations can generate similar results when a ""conflicting zeitgeber experiment'' incorporates only two phase relationships between zeitgebers.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

CHARACTERIZATION OF EXPONENTIAL DIVERGENCE OF THE KALMAN FILTER FOR TIME-VARYING SYSTEMS

This paper studies semistability of the recursive Kalman filter in the context of linear time-varying (LTV), possibly nondetectable systems with incorrect noise information. Semistability is a key property, as it ensures that the actual estimation error does not diverge exponentially. We explore structural properties of the filter to obtain a necessary and sufficient condition for the filter to be semistable. The condition does not involve limiting gains nor the solution of Riccati equations, as they can be difficult to obtain numerically and may not exist. We also compare semistability with the notions of stability and stability w.r.t. the initial error covariance, and we show that semistability in a sense makes no distinction between persistent and nonpersistent incorrect noise models, as opposed to stability. In the linear time invariant scenario we obtain algebraic, easy to test conditions for semistability and stability, which complement results available in the context of detectable systems. Illustrative examples are included.

Korteveg-de Vries solitons in a cold quark-gluon plasma

The relativistic heavy ion program developed at RHIC and now at LHC motivated a deeper study of the properties of the quark-gluon plasma (QGP) and, in particular, the study of perturbations in this kind of plasma. We are interested on the time evolution of perturbations in the baryon and energy densities. If a localized pulse in baryon density could propagate throughout the QGP for long distances preserving its shape and without loosing localization, this could have interesting consequences for relativistic heavy ion physics and for astrophysics. A mathematical way to prove that this can happen is to derive (under certain conditions) from the hydrodynamical equations of the QGP a Korteveg-de Vries (KdV) equation. The solution of this equation describes the propagation of a KdV soliton. The derivation of the KdV equation depends crucially on the equation of state (EOS) of the QGP. The use of the simple MIT bag model EOS does not lead to KdV solitons. Recently we have developed an EOS for the QGP which includes both perturbative and nonperturbative corrections to the MIT one and is still simple enough to allow for analytical manipulations. With this EOS we were able to derive a KdV equation for the cold QGP.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.