Finite time blow-up and breaking of solitary wind waves

The evolution of surface water waves in finite depth under wind forcing is reduced to an antidissipative Korteweg-de Vries-Burgers equation. We exhibit its solitary wave solution. Antidissipation accelerates and increases the amplitude of the solitary wave and leads to blow-up and breaking. Blow-up occurs in finite time for infinitely large asymptotic space so it is a nonlinear, dispersive, and antidissipative equivalent of the linear instability which occurs for infinite time. Due to antidissipation two given arbitrary and adjacent planes of constant phases of the solitary wave acquire different velocities and accelerations inducing breaking. Soliton breaking occurs in finite space in a time prior to the blow-up. We show that the theoretical growth in amplitude and the time of breaking are both testable in an existing experimental facility.

A family of dissipative two-dimensional mappings: Chaotic, regular and steady state dynamics investigation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A dissipative Fermi-Ulam model under two different kinds of dissipation

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

About the minimum time transference in a fractional differential equation

The use of fractional calculus when modeling phenomena allows new queries concerning the deepest parts of the physical laws involved in. Here we will be dealing with an apparent paradox in which the time of transference from zero in a system with fractional derivatives can be strictly shortened relatively to the minimal time transference done in an equivalent system in the frame of the entire derivatives.

Ordinary fractional differential equations are in fact usual entire ordinary differential equations on time scales

Despite the huge number of works considering fractional derivatives or derivatives on time scales some basic facts remain to be evaluated. Here we will be showing that the fractional derivative of monomials is in fact an entire derivative considered on an appropriate time scale.

Space-time description of dengue epidemic in Cruzeiro, SP, in 2006 and 2011

Objective: to identify patterns in the spatial and temporal distribution of cases of dengue fever occurring in the city of Cruzeiro, state of Sao Paulo (SP).Methods: an ecological and exploratory study was undertaken using spatial analysis tools and data from dengue cases obtained on the SinanNet. The analysis was carried out by area, using the IBGE census sector as a unit. The months of March to June 2006 and 2011 were assessed, revealing progress of the disease. TerraView 3.3.1 was used to calculate the Global Moran's I, month to month, and the Kernel estimator.Results: in the year 2006, 691 cases of dengue fever (rate of 864.2 cases/100,000 inhabitants) were georeferenced; and the Moran's I and p-values were significant in the months of April and May (TM = 0.28; p = 0.01; I-M = 0.20; p = 0.01) with higher densities in the central, north, northeast and south regions. In the year 2011, 654 cases of dengue fever (rate of 886.8 cases/100,000 inhabitants) were georeferenced; and the Moran's I and p-values were significant in the months of April and May (I, = 0.28; p = 0.01; I-M = 0.16; p = 0.05) with densities in the same regions as 2006. The Global Moran's I is a global measure of spatial autocorrelation, which indicates the degree of spatial association in the set of information from the product in relation to the average. The I varies between -1 and +1 and can be attributed to a level of significance (p-value). The positive value points to a positive or direct spatial autocorrelation.Conclusion: we were able to identify patterns in the spatial and temporal distribution of dengue cases occurring in the city of Cruzeiro, SP, and locate the census sectors where the outbreak began and how it evolved.

A two-layer approximation to the Brazil Current-Intermediate Western Boundary Current System between 20 degrees S and 28 degrees S

It has been shown that the vertical structure of the Brazil Current (BC)-Intermediate Western Boundary Current (IWBC) System is dominated by the first baroclinic mode at 22 degrees S-23 degrees S. In this work, we employed the Miami Isopycnic Coordinate Ocean Model to investigate whether the rich mesoscale activity of this current system, between 20 degrees S and 28 degrees S, is reproduced by a two-layer approximation of its vertical structure. The model results showed cyclonic and anticyclonic meanders propagating southwestward along the current axis, resembling the dynamical pattern of Rossby waves superposed on a mean flow. Analysis of the upper layer zonal velocity component, using a space-time diagram, revealed a dominant wavelength of about 450 km and phase velocity of about 0.20 ms(-1) southwestward. The results also showed that the eddy-like structures slowly grew in amplitude as they moved downstream. Despite the simplified design of the numerical experiments conducted here, these results compared favorably with observations and seem to indicate that weakly unstable long baroclinic waves are responsible for most of the variability observed in the BC-IWBC system. (C) 2009 Elsevier Ltd. All rights reserved.

Color-space distortions following long-term occupational exposure to mercury vapor

Color vision was examined in subjects with long-term occupational exposure to mercury (Hg) vapor. The color vision impairment was assessed by employing a quantitative measure of distortion of individual and group perceptual color spaces. Hg subjects (n = 18; 42.1 +/- 6.5 years old; exposure time = 10.4 +/- 5.0 years; time away from the exposure source = 6.8 +/- 4.6 years) and controls (n = 18; 46.1 +/- 8.4 years old) were examined using two arrangement tests, D-15 and D-15d, in the traditional way, and also in a triadic procedure. From each subject`s `odd-one-out` choices, matrices of inter-cap subjective dissimilarities were derived and processed by non-metric multidimensional scaling (MDS). D-15d results differed significantly between the Hg-group and the control group (p < 0.05), with the impairment predominantly along the tritan axis. 2D perceptual color spaces, individual and group, were reconstructed, with the dimensions interpreted as the red-green (RG) and the blue-yellow (BY) systems. When color configurations from the Hg-group were compared to those of the controls, they presented more fluctuations along both chromatic dimensions, indicating a statistically significant difference along the BY axis. In conclusion, the present findings confirm that color vision impairments persist in subjects that have received long-term occupational exposure to Hg-vapor although, at the time of testing, they were presenting mean urinary concentration within the normal range for non-exposed individuals. Considering the advantages of the triadic procedure in clinical evaluation of acquired color vision deficiencies, further studies should attempt to verify and/or improve its efficacy.

Using occupancy models to investigate space partitioning between two sympatric large predators, the jaguar and puma in central Brazil

Coexistence of sympatric species is mediated by resource partitioning. Pumas occur sympatrically with jaguars throughout most of the jaguar's range but few studies have investigated space partitioning between both species. Here, camera trapping and occupancy models accounting for imperfect detection were employed in a Bayesian framework to investigate space partitioning between the jaguar and puma in Emas National Park (ENP), central Brazil. Jaguars were estimated to occupy 54.1% and pumas 39.3% of the sample sites. Jaguar occupancy was negatively correlated with distance to water and positively correlated with the amount of dense habitat surrounding the camera trap. Puma occupancy only showed a weak negative correlation with distance to water and with jaguar presence. Both species were less often present at the same site than expected under independent distributions. Jaguars had a significantly higher detection probability at cameras on roads than at off-road locations. For pumas, detection was similar on and off-road. Results indicate that both differences in habitat use and active avoidance shape space partitioning between jaguars and pumas in ENP. Considering its size, the jaguar is likely the competitively dominant of the two species. Owing to its habitat preferences, suitable jaguar habitat outside the park is probably sparse. Consequently, the jaguar population is likely largely confined to the park, while the puma population is known to extend into ENP's surroundings. (C) 2011 Deutsche Gesellschaft fur Saugetierkunde. Published by Elsevier GmbH. All rights reserved.

Identifying regulational alterations in gene regulatory networks by state space representation of vector autoregressive models and variational annealing

Background: In the analysis of effects by cell treatment such as drug dosing, identifying changes on gene network structures between normal and treated cells is a key task. A possible way for identifying the changes is to compare structures of networks estimated from data on normal and treated cells separately. However, this approach usually fails to estimate accurate gene networks due to the limited length of time series data and measurement noise. Thus, approaches that identify changes on regulations by using time series data on both conditions in an efficient manner are demanded. Methods: We propose a new statistical approach that is based on the state space representation of the vector autoregressive model and estimates gene networks on two different conditions in order to identify changes on regulations between the conditions. In the mathematical model of our approach, hidden binary variables are newly introduced to indicate the presence of regulations on each condition. The use of the hidden binary variables enables an efficient data usage; data on both conditions are used for commonly existing regulations, while for condition specific regulations corresponding data are only applied. Also, the similarity of networks on two conditions is automatically considered from the design of the potential function for the hidden binary variables. For the estimation of the hidden binary variables, we derive a new variational annealing method that searches the configuration of the binary variables maximizing the marginal likelihood. Results: For the performance evaluation, we use time series data from two topologically similar synthetic networks, and confirm that our proposed approach estimates commonly existing regulations as well as changes on regulations with higher coverage and precision than other existing approaches in almost all the experimental settings. For a real data application, our proposed approach is applied to time series data from normal Human lung cells and Human lung cells treated by stimulating EGF-receptors and dosing an anticancer drug termed Gefitinib. In the treated lung cells, a cancer cell condition is simulated by the stimulation of EGF-receptors, but the effect would be counteracted due to the selective inhibition of EGF-receptors by Gefitinib. However, gene expression profiles are actually different between the conditions, and the genes related to the identified changes are considered as possible off-targets of Gefitinib. Conclusions: From the synthetically generated time series data, our proposed approach can identify changes on regulations more accurately than existing methods. By applying the proposed approach to the time series data on normal and treated Human lung cells, candidates of off-target genes of Gefitinib are found. According to the published clinical information, one of the genes can be related to a factor of interstitial pneumonia, which is known as a side effect of Gefitinib.

Three time-based scale formulations for the two-stage lot sizing and scheduling in process industries

In this paper, we propose three novel mathematical models for the two-stage lot-sizing and scheduling problems present in many process industries. The problem shares a continuous or quasi-continuous production feature upstream and a discrete manufacturing feature downstream, which must be synchronized. Different time-based scale representations are discussed. The first formulation encompasses a discrete-time representation. The second one is a hybrid continuous-discrete model. The last formulation is based on a continuous-time model representation. Computational tests with state-of-the-art MIP solver show that the discrete-time representation provides better feasible solutions in short running time. On the other hand, the hybrid model achieves better solutions for longer computational times and was able to prove optimality more often. The continuous-type model is the most flexible of the three for incorporating additional operational requirements, at a cost of having the worst computational performance. Journal of the Operational Research Society (2012) 63, 1613-1630. doi:10.1057/jors.2011.159 published online 7 March 2012

Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

Closed-loop stable model predictive control of integrating systems with dead time

Model predictive control (MPC) applications in the process industry usually deal with process systems that show time delays (dead times) between the system inputs and outputs. Also, in many industrial applications of MPC, integrating outputs resulting from liquid level control or recycle streams need to be considered as controlled outputs. Conventional MPC packages can be applied to time-delay systems but stability of the closed loop system will depend on the tuning parameters of the controller and cannot be guaranteed even in the nominal case. In this work, a state space model based on the analytical step response model is extended to the case of integrating time systems with time delays. This model is applied to the development of two versions of a nominally stable MPC, which is designed to the practical scenario in which one has targets for some of the inputs and/or outputs that may be unreachable and zone control (or interval tracking) for the remaining outputs. The controller is tested through simulation of a multivariable industrial reactor system. (C) 2012 Elsevier Ltd. All rights reserved.

Modeling the gluon propagator in Landau gauge: Lattice estimates of pole masses and dimension-two condensates

We present an analytic description of numerical results for the Landau-gauge SU(2) gluon propagator D(p(2)), obtained from lattice simulations (in the scaling region) for the largest lattice sizes to date, in d = 2, 3 and 4 space-time dimensions. Fits to the gluon data in 3d and in 4d show very good agreement with the tree-level prediction of the refined Gribov-Zwanziger (RGZ) framework, supporting a massive behavior for D(p(2)) in the infrared limit. In particular, we investigate the propagator's pole structure and provide estimates of the dynamical mass scales that can be associated with dimension-two condensates in the theory. In the 2d case, fitting the data requires a noninteger power of the momentum p in the numerator of the expression for D(p(2)). In this case, an infinite-volume-limit extrapolation gives D(0) = 0. Our analysis suggests that this result is related to a particular symmetry in the complex-pole structure of the propagator and not to purely imaginary poles, as would be expected in the original Gribov-Zwanziger scenario.

SINGULARLY PERTURBED DISCOUNTED MARKOV CONTROL PROCESSES IN A GENERAL STATE SPACE

This paper studies the asymptotic optimality of discrete-time Markov decision processes (MDPs) with general state space and action space and having weak and strong interactions. By using a similar approach as developed by Liu, Zhang, and Yin [Appl. Math. Optim., 44 (2001), pp. 105-129], the idea in this paper is to consider an MDP with general state and action spaces and to reduce the dimension of the state space by considering an averaged model. This formulation is often described by introducing a small parameter epsilon > 0 in the definition of the transition kernel, leading to a singularly perturbed Markov model with two time scales. Our objective is twofold. First it is shown that the value function of the control problem for the perturbed system converges to the value function of a limit averaged control problem as epsilon goes to zero. In the second part of the paper, it is proved that a feedback control policy for the original control problem defined by using an optimal feedback policy for the limit problem is asymptotically optimal. Our work extends existing results of the literature in the following two directions: the underlying MDP is defined on general state and action spaces and we do not impose strong conditions on the recurrence structure of the MDP such as Doeblin's condition.