COMPARISON RESULTS FOR CONJUGATE AND FOCAL POINTS IN SEMI-RIEMANNIAN GEOMETRY VIA MASLOV INDEX

We prove an estimate on the difference of Maslov indices relative to the choice of two distinct reference Lagrangians of a continuous path in the Lagrangian Grassmannian of a symplectic space. We discuss some applications to the study of conjugate and focal points along a geodesic in a semi-Riemannian manifold.

A Lagrangian identification of major sources of moisture over Central Brazil and La Plata Basin

[1] This work examines the main sources of moisture over Central Brazil and La Plata Basin during the year through a new Lagrangian diagnosis method which identifies the humidity contributions to the moisture budget over a region. This methodology computes budgets of evaporation minus precipitation by calculating changes in the specific humidity along back-trajectories for the previous 10 d. The origin of all air masses residing over each region was tracked during a period of 5 years (2000-2004). These regions were selected because they coincide with two centers of action of a known dipole precipitation variability mode observed in different temporal scales (from intra seasonal up to inter decadal timescales) and are related to the climatic variability of the South American Monsoon System. The results suggested the importance of the tropical south Atlantic as a moisture source for Central Brazil, and of recycling for La Plata basin. It seems that the Tropical South Atlantic plays an important role as a moisture source for Central Brazil and La Plata basin along the year, particularly during the austral summer. The north Atlantic is also an additional source for both regions during the austral summer.

Dynamics of two planets in co-orbital motion

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities; we also assume that both planets share the same orbital plane. Initially, we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analysed in more detail using a semi-analytical model. Apart from the well-known quasi-satellite orbits and the classical equilibrium Lagrangian points L(4) and L(5), we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at (delta lambda, delta pi) = (+/- 60 degrees, -/+ 120 degrees), where delta lambda is the difference in mean longitudes and delta pi is the difference in longitudes of pericentre. The position of these anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities and are found for eccentricities as high as similar to 0.7. Finally, we also applied a slow mass variation to one of the planets and analysed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.

Dietary copper absorption and excretion in three semi-terrestrial grapsoid crabs with different levels of terrestrial adaptation

Three species of phylogenetically related semi-terrestrial crabs (Superfamily Grapsoidea - Sesarma rectum, Goniopsis cruentata and Neohelice granulata (formerly: Chasmagnathus granulatus) with different degrees of terrestriality were studied to quantify the accumulation of copper (Cu) in hemolymph, gills, hepatopancreas and antennal gland, and its excretion through the faeces. These crabs were fed for 15 days practical diets containing 0 (A), 0.5 (B), 1.0 (C), and 1.5% (D) of added CuCl2 (corresponding to 0, 0.2, 0.5 and 0.7% of Cu2+, respectively). The amount of food ingested was directly proportional to the degree of terrestriality: S. rectum, the most terrestrial species, ate around 2-3 times more than the other crabs, whereas G. cruentata ate 1.5-2 times more than N. granulata, the least terrestrial. The amount of Cu excreted in the feces was proportional to Cu ingestion, and was 76.8% and 64.2% higher for Sesarma fed diet D compared to G. cruentata and N. granulata, respectively. Sesarma also displayed higher Cu concentration in the haemolymph, gills and antennal glands, but not in the hepatopancreas. A detoxifying mechanism followed by elimination was probably present at this last organ, preventing Cu accumulation. More terrestrial crabs, such as Sesarma, may accumulate more Cu in hemolymph and tissues, showing a correlation between metal accumulation and increased terrestriality. In this aspect, contaminated feed sources with Cu may have more impact in conservation of terrestrial crabs. (C) 2008 Elsevier Inc. All rights reserved.

The importance of dietary calcium consumption in two species of semi-terrestrial grapsoid crabs

Calcium (Ca) is essential for crustaceans, due to calcium carbonate (CaCO(3)) deposition in the new exoskeleton to harden it. The purpose of this work was to study short term Ca balance in terms of dietary Ca ingestion in two phylogenetically related crabs (Superfamily Grapsoidea) showing different degrees of terrestrial adaptations: Sesarma rectum Randall, 1840 and Neohelice granulata (Dana, 1851). Dietary Ca ingestion was studied using purified diets with different Ca concentrations (0, 2.2 and 6.66% Ca), together with measurements of Ca excretion and Ca hemolymph levels. The results showed that both crabs had the same response to foods containing different levels of Ca, with both species eating more of the high Ca diet. However, S. rectum consumed more per mg body mass at all Ca concentrations (6 mg.g(-1) for S. rectum against 3 mg.g(-1) for N. granulata). Both species excreted/egested Ca differently: S. rectum excreted Ca proportionally to ingestion, whereas N. granulata maintained constant faecal Ca output at all dietary Ca levels. Moreover, Ca hemolymph levels for crabs fed the different diets were independent of dietary Ca. In conclusion, both S. rectum and N. granulata seem to regulate the consumption of diets containing more Ca, which suggests a fine balance for Ca intake.

Donor bone marrow cells play a role in the prevention of accelerated graft rejection induced by semi-allogeneic spleen cells in transplantation

Spleen or spleen plus bone marrow cells from (BALB/c x C57Bl/6)F1 donors were transferred into BALB/c recipients 21 days before skin or cardiac transplantation. Prolonged graft survival was observed on recipients treated with the mixture of donor-derived cells as compared to those treated with spleen cells alone. We evaluated the expression of CD45RB and CD44 by splenic CD4(+) and CD8(+) T cells 7 and 21 days after donor cell transfer. The populations of CD8(+)CD45RB(low) and CD8(+)CD44(high) cells were significantly decreased in mice pre-treated with donor spleen and bone marrow cells as compared to animals treated with spleen cells only, although these cells expanded in both groups when compared to an earlier time-point. No differences were observed regarding CD4+ T cell population when recipients of donor-derived cells were compared. An enhanced production of IL-10 was observed seven days after transplantation in the supernatants of spleen cell cultures of mice treated with spleen and bone marrow cells. Taken together these data suggest that donor-derived bone marrow cells modulate the sensitization of the recipient by semi-allogeneic spleen cells in part by delaying the generation of activated/memory CD8(+) T cells leading to enhanced graft survival. (c) 2007 Elsevier B.V. All rights reserved.

Reduced gradient method combined with augmented Lagrangian and barrier for the optimal power flow problem

A new approach for solving the optimal power flow (OPF) problem is established by combining the reduced gradient method and the augmented Lagrangian method with barriers and exploring specific characteristics of the relations between the variables of the OPF problem. Computer simulations on IEEE 14-bus and IEEE 30-bus test systems illustrate the method. (c) 2007 Elsevier Inc. All rights reserved.

A Lagrangian relaxation approach to a coupled lot-sizing and cutting stock problem

Industrial production processes involving both lot-sizing and cutting stock problems are common in many industrial settings. However, they are usually treated in a separate way, which could lead to costly production plans. In this paper, a coupled mathematical model is formulated and a heuristic method based on Lagrangian relaxation is proposed. Computational results prove its effectiveness. (C) 2009 Elsevier B.V. All rights reserved.

Aerosol inorganic composition at a tropical site: Discrepancies between filter-based sampling and a semi-continuous method

The concentrations of the water-soluble inorganic aerosol species, ammonium (NH4+), nitrate (NO3-), chloride (Cl-), and sulfate (SO42-), were measured from September to November 2002 at a pasture site in the Amazon Basin (Rondnia, Brazil) (LBA-SMOCC). Measurements were conducted using a semi-continuous technique (Wet-annular denuder/Steam-Jet Aerosol Collector: WAD/SJAC) and three integrating filter-based methods, namely (1) a denuder-filter pack (DFP: Teflon and impregnated Whatman filters), (2) a stacked-filter unit (SFU: polycarbonate filters), and (3) a High Volume dichotomous sampler (HiVol: quartz fiber filters). Measurements covered the late dry season (biomass burning), a transition period, and the onset of the wet season (clean conditions). Analyses of the particles collected on filters were performed using ion chromatography (IC) and Particle-Induced X-ray Emission spectrometry (PIXE). Season-dependent discrepancies were observed between the WAD/SJAC system and the filter-based samplers. During the dry season, when PM2.5 (D-p <= 2.5 mu m) concentrations were similar to 100 mu g m(-3), aerosol NH4+ and SO42- measured by the filter-based samplers were on average two times higher than those determined by the WAD/SJAC. Concentrations of aerosol NO3- and Cl- measured with the HiVol during daytime, and with the DFP during day- and nighttime also exceeded those of the WAD/SJAC by a factor of two. In contrast, aerosol NO3- and Cl- measured with the SFU during the dry season were nearly two times lower than those measured by the WAD/SJAC. These differences declined markedly during the transition period and towards the cleaner conditions during the onset of the wet season (PM2.5 similar to 5 mu g m(-3)); when filter-based samplers measured on average 40-90% less than the WAD/SJAC. The differences were not due to consistent systematic biases of the analytical techniques, but were apparently a result of prevailing environmental conditions and different sampling procedures. For the transition period and wet season, the significance of our results is reduced by a low number of data points. We argue that the observed differences are mainly attributable to (a) positive and negative filter sampling artifacts, (b) presence of organic compounds and organosulfates on filter substrates, and (c) a SJAC sampling efficiency of less than 100%.

Global minimization using an Augmented Lagrangian method with variable lower-level constraints

A novel global optimization method based on an Augmented Lagrangian framework is introduced for continuous constrained nonlinear optimization problems. At each outer iteration k the method requires the epsilon(k)-global minimization of the Augmented Lagrangian with simple constraints, where epsilon(k) -> epsilon. Global convergence to an epsilon-global minimizer of the original problem is proved. The subproblems are solved using the alpha BB method. Numerical experiments are presented.

Structured minimal-memory inexact quasi-Newton method and secant preconditioners for augmented Lagrangian optimization

Augmented Lagrangian methods for large-scale optimization usually require efficient algorithms for minimization with box constraints. On the other hand, active-set box-constraint methods employ unconstrained optimization algorithms for minimization inside the faces of the box. Several approaches may be employed for computing internal search directions in the large-scale case. In this paper a minimal-memory quasi-Newton approach with secant preconditioners is proposed, taking into account the structure of Augmented Lagrangians that come from the popular Powell-Hestenes-Rockafellar scheme. A combined algorithm, that uses the quasi-Newton formula or a truncated-Newton procedure, depending on the presence of active constraints in the penalty-Lagrangian function, is also suggested. Numerical experiments using the Cute collection are presented.

Improving ultimate convergence of an augmented Lagrangian method

Optimization methods that employ the classical Powell-Hestenes-Rockafellar augmented Lagrangian are useful tools for solving nonlinear programming problems. Their reputation decreased in the last 10 years due to the comparative success of interior-point Newtonian algorithms, which are asymptotically faster. In this research, a combination of both approaches is evaluated. The idea is to produce a competitive method, being more robust and efficient than its `pure` counterparts for critical problems. Moreover, an additional hybrid algorithm is defined, in which the interior-point method is replaced by the Newtonian resolution of a Karush-Kuhn-Tucker (KKT) system identified by the augmented Lagrangian algorithm. The software used in this work is freely available through the Tango Project web page:http://www.ime.usp.br/similar to egbirgin/tango/.

Augmented Lagrangian methods under the constant positive linear dependence constraint qualification

Two Augmented Lagrangian algorithms for solving KKT systems are introduced. The algorithms differ in the way in which penalty parameters are updated. Possibly infeasible accumulation points are characterized. It is proved that feasible limit points that satisfy the Constant Positive Linear Dependence constraint qualification are KKT solutions. Boundedness of the penalty parameters is proved under suitable assumptions. Numerical experiments are presented.

A semi-Markov multistate model for estimation of the mean quality-adjusted survival for non-progressive processes

We discuss the estimation of the expected value of the quality-adjusted survival, based on multistate models. We generalize an earlier work, considering the sojourn times in health states are not identically distributed, for a given vector of covariates. Approaches based on semiparametric and parametric (exponential and Weibull distributions) methodologies are considered. A simulation study is conducted to evaluate the performance of the proposed estimator and the jackknife resampling method is used to estimate the variance of such estimator. An application to a real data set is also included.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.