Reformulações e relaxação Lagrangiana para o problema de dimensionamento de lotes com várias plantas

Pós-graduação em Matemática - IBILCE

Um panorama sobre roll waves em escoamentos laminares e turbulentos com superfície livre

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

An integrable evolution equation for surface waves in deep water

In order to describe the dynamics of monochromatic surface waves in deep water, we derive a nonlinear and dispersive system of equations for the free surface elevation and the free surface velocity from the Euler equations in infinite depth. From it, and using a multiscale perturbative method, an asymptotic model for small wave steepness ratio is derived. The model is shown to be completely integrable. The Lax pair, the first conserved quantities as well as the symmetries are exhibited. Theoretical and numerical studies reveal that it supports periodic progressive Stokes waves which peak and break in finite time. Comparison between the limiting wave solution of the asymptotic model and classical results is performed.

Otimização ecológica dos ciclos ar-padrão Otto e Diesel

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Radiação e detetores de Unruh - DeWitt

Pós-graduação em Física - IFT

Caos e termalização na teoria de Yang-Mills-Higgs em uma rede espacial

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

Effective transport barriers in nontwist systems

In fluids and plasmas with zonal flow reversed shear, a peculiar kind of transport barrier appears in the shearless region, one that is associated with a proper route of transition to chaos. These barriers have been identified in symplectic nontwist maps that model such zonal flows. We use the so-called standard nontwist map, a paradigmatic example of nontwist systems, to analyze the parameter dependence of the transport through a broken shearless barrier. On varying a proper control parameter, we identify the onset of structures with high stickiness that give rise to an effective barrier near the broken shearless curve. Moreover, we show how these stickiness structures, and the concomitant transport reduction in the shearless region, are determined by a homoclinic tangle of the remaining dominant twin island chains. We use the finite-time rotation number, a recently proposed diagnostic, to identify transport barriers that separate different regions of stickiness. The identified barriers are comparable to those obtained by using finite-time Lyapunov exponents.

Stickiness in a bouncer model: A slowing mechanism for Fermi acceleration

Some phase space transport properties for a conservative bouncer model are studied. The dynamics of the model is described by using a two-dimensional measure preserving mapping for the variables' velocity and time. The system is characterized by a control parameter epsilon and experiences a transition from integrable (epsilon = 0) to nonintegrable (epsilon not equal 0). For small values of epsilon, the phase space shows a mixed structure where periodic islands, chaotic seas, and invariant tori coexist. As the parameter epsilon increases and reaches a critical value epsilon(c), all invariant tori are destroyed and the chaotic sea spreads over the phase space, leading the particle to diffuse in velocity and experience Fermi acceleration (unlimited energy growth). During the dynamics the particle can be temporarily trapped near periodic and stable regions. We use the finite time Lyapunov exponent to visualize this effect. The survival probability was used to obtain some of the transport properties in the phase space. For large epsilon, the survival probability decays exponentially when it turns into a slower decay as the control parameter epsilon is reduced. The slower decay is related to trapping dynamics, slowing the Fermi Acceleration, i.e., unbounded growth of the velocity.

Quantum motion in superposition of Aharonov-Bohm with some additional electromagnetic fields

The structure of additional electromagnetic fields to the Aharonov-Bohm field, for which the Schrodinger, Klein-Gordon, and Dirac equations can be solved exactly are described and the corresponding exact solutions are found. It is demonstrated that aside from the known cases (a constant and uniform magnetic field that is parallel to the Aharonov-Bohm solenoid, a static spherically symmetrical electric field, and the field of a magnetic monopole), there are broad classes of additional fields. Among these new additional fields we have physically interesting electric fields acting during a finite time or localized in a restricted region of space. There are additional time-dependent uniform and isotropic electric fields that allow exact solutions of the Schrodinger equation. In the relativistic case there are additional electric fields propagating along the Aharonov-Bohm solenoid with arbitrary electric pulse shape. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4714352]

Imaging of crustal scatterers using multiple seismic arrays

Array seismology is an useful tool to perform a detailed investigation of the Earth’s interior. Seismic arrays by using the coherence properties of the wavefield are able to extract directivity information and to increase the ratio of the coherent signal amplitude relative to the amplitude of incoherent noise. The Double Beam Method (DBM), developed by Krüger et al. (1993, 1996), is one of the possible applications to perform a refined seismic investigation of the crust and mantle by using seismic arrays. The DBM is based on a combination of source and receiver arrays leading to a further improvement of the signal-to-noise ratio by reducing the error in the location of coherent phases. Previous DBM works have been performed for mantle and core/mantle resolution (Krüger et al., 1993; Scherbaum et al., 1997; Krüger et al., 2001). An implementation of the DBM has been presented at 2D large-scale (Italian data-set for Mw=9.3, Sumatra earthquake) and at 3D crustal-scale as proposed by Rietbrock & Scherbaum (1999), by applying the revised version of Source Scanning Algorithm (SSA; Kao & Shan, 2004). In the 2D application, the rupture front propagation in time has been computed. In 3D application, the study area (20x20x33 km3), the data-set and the source-receiver configurations are related to the KTB-1994 seismic experiment (Jost et al., 1998). We used 60 short-period seismic stations (200-Hz sampling rate, 1-Hz sensors) arranged in 9 small arrays deployed in 2 concentric rings about 1 km (A-arrays) and 5 km (B-array) radius. The coherence values of the scattering points have been computed in the crustal volume, for a finite time-window along all array stations given the hypothesized origin time and source location. The resulting images can be seen as a (relative) joint log-likelihood of any point in the subsurface that have contributed to the full set of observed seismograms.

Dosimetric impact of geometric errors due to respiratory motion prediction on dynamic multileaf collimator-based four-dimensional radiation delivery

The synchronization of dynamic multileaf collimator (DMLC) response with respiratory motion is critical to ensure the accuracy of DMLC-based four dimensional (4D) radiation delivery. In practice, however, a finite time delay (response time) between the acquisition of tumor position and multileaf collimator response necessitates predictive models of respiratory tumor motion to synchronize radiation delivery. Predicting a complex process such as respiratory motion introduces geometric errors, which have been reported in several publications. However, the dosimetric effect of such errors on 4D radiation delivery has not yet been investigated. Thus, our aim in this work was to quantify the dosimetric effects of geometric error due to prediction under several different conditions. Conformal and intensity modulated radiation therapy (IMRT) plans for a lung patient were generated for anterior-posterior/posterior-anterior (AP/PA) beam arrangements at 6 and 18 MV energies to provide planned dose distributions. Respiratory motion data was obtained from 60 diaphragm-motion fluoroscopy recordings from five patients. A linear adaptive filter was employed to predict the tumor position. The geometric error of prediction was defined as the absolute difference between predicted and actual positions at each diaphragm position. Distributions of geometric error of prediction were obtained for all of the respiratory motion data. Planned dose distributions were then convolved with distributions for the geometric error of prediction to obtain convolved dose distributions. The dosimetric effect of such geometric errors was determined as a function of several variables: response time (0-0.6 s), beam energy (6/18 MV), treatment delivery (3D/4D), treatment type (conformal/IMRT), beam direction (AP/PA), and breathing training type (free breathing/audio instruction/visual feedback). Dose difference and distance-to-agreement analysis was employed to quantify results. Based on our data, the dosimetric impact of prediction (a) increased with response time, (b) was larger for 3D radiation therapy as compared with 4D radiation therapy, (c) was relatively insensitive to change in beam energy and beam direction, (d) was greater for IMRT distributions as compared with conformal distributions, (e) was smaller than the dosimetric impact of latency, and (f) was greatest for respiration motion with audio instructions, followed by visual feedback and free breathing. Geometric errors of prediction that occur during 4D radiation delivery introduce dosimetric errors that are dependent on several factors, such as response time, treatment-delivery type, and beam energy. Even for relatively small response times of 0.6 s into the future, dosimetric errors due to prediction could approach delivery errors when respiratory motion is not accounted for at all. To reduce the dosimetric impact, better predictive models and/or shorter response times are required.

Importance sampling approximations to various probabilities of ruin of spectrally negative Lévy risk processes

This article provides importance sampling algorithms for computing the probabilities of various types ruin of spectrally negative Lévy risk processes, which are ruin over the infinite time horizon, ruin within a finite time horizon and ruin past a finite time horizon. For the special case of the compound Poisson process perturbed by diffusion, algorithms for computing probabilities of ruins by creeping (i.e. induced by the diffusion term) and by jumping (i.e. by a claim amount) are provided. It is shown that these algorithms have either bounded relative error or logarithmic efficiency, as t,x→∞t,x→∞, where t>0t>0 is the time horizon and x>0x>0 is the starting point of the risk process, with y=t/xy=t/x held constant and assumed either below or above a certain constant.

Evolution and breakup of viscous rotating drops

We study the evolution of a viscous fluid drop rotating about a fixed axis at constant angular velocity $Omega$ or constant angular momentum L surrounded by another viscous fluid. The problem is considered in the limit of large Ekman number and small Reynolds number. The analysis is carried out by combining asymptotic analysis and full numerical simulation by means of the boundary element method. We pay special attention to the stability/instability of equilibrium shapes and the possible formation of singularities representing a change in the topology of the fluid domain. When the evolution is at constant $Omega$, depending on its value, drops can take the form of a flat film whose thickness goes to zero in finite time or an elongated filament that extends indefinitely. When evolution takes place at constant L and axial symmetry is imposed, thin films surrounded by a toroidal rim can develop, but the film thickness does not vanish in finite time. When axial symmetry is not imposed and L is sufficiently large, drops break axial symmetry and, depending on the value of L, reach an equilibrium configuration with a 2-fold symmetry or break up into several drops with a 2- or 3-fold symmetry. The mechanism of breakup is also described

The Lagrangian description of aperiodic flows: a case study of the Kuroshio Current

This article reviews several recently developed Lagrangian tools and shows how their com- bined use succeeds in obtaining a detailed description of purely advective transport events in general aperiodic flows. In particular, because of the climate impact of ocean transport processes, we illustrate a 2D application on altimeter data sets over the area of the Kuroshio Current, although the proposed techniques are general and applicable to arbitrary time depen- dent aperiodic flows. The first challenge for describing transport in aperiodical time dependent flows is obtaining a representation of the phase portrait where the most relevant dynamical features may be identified. This representation is accomplished by using global Lagrangian descriptors that when applied for instance to the altimeter data sets retrieve over the ocean surface a phase portrait where the geometry of interconnected dynamical systems is visible. The phase portrait picture is essential because it evinces which transport routes are acting on the whole flow. Once these routes are roughly recognised it is possible to complete a detailed description by the direct computation of the finite time stable and unstable manifolds of special hyperbolic trajectories that act as organising centres of the flow.

Active Learning and Dynamic Pricing Policies

In this paper, we address the problem of dynamic pricing to optimize the revenue coming from the sales of a limited inventory in a finite time-horizon. A priori, the demand is assumed to be unknown. The seller must learn on the fly. We first deal with the simplest case, involving only one class of product for sale. Furthermore the general situation is considered with a finite number of product classes for sale. In particular, a case in point is the sale of tickets for events related to culture and leisure; in this case, typically the tickets are sold months before the event, thus, uncertainty over actual demand levels is a very a common occurrence. We propose a heuristic strategy of adaptive dynamic pricing, based on experience gained from the past, taking into account, for each time period, the available inventory, the time remaining to reach the horizon, and the profit made in previous periods. In the computational simulations performed, the demand is updated dynamically based on the prices being offered, as well as on the remaining time and inventory. The simulations show a significant profit over the fixed-price strategy, confirming the practical usefulness of the proposed strategy. We develop a tool allowing us to test different dynamic pricing strategies designed to fit market conditions and seller s objectives, which will facilitate data analysis and decision-making in the face of the problem of dynamic pricing.