A relaxed constant positive linear dependence constraint qualification and applications

In this work we introduce a relaxed version of the constant positive linear dependence constraint qualification (CPLD) that we call RCPLD. This development is inspired by a recent generalization of the constant rank constraint qualification by Minchenko and Stakhovski that was called RCRCQ. We show that RCPLD is enough to ensure the convergence of an augmented Lagrangian algorithm and that it asserts the validity of an error bound. We also provide proofs and counter-examples that show the relations of RCRCQ and RCPLD with other known constraint qualifications. In particular, RCPLD is strictly weaker than CPLD and RCRCQ, while still stronger than Abadie's constraint qualification. We also verify that the second order necessary optimality condition holds under RCRCQ.

Shearless transport barriers in magnetically confined plasmas

Shearless transport barriers appear in confined plasmas due to non-monotonic radial profiles and cause localized reduction of transport even after they have been broken. In this paper we summarize our recent theoretical and experimental research on shearless transport barriers in plasmas confined in toroidal devices. In particular, we discuss shearless barriers in Lagrangian magnetic field line transport caused by non-monotonic safety factor profiles. We also discuss evidence of particle transport barriers found in the TCABR Tokamak (University of Sao Paulo) and the Texas Helimak (University of Texas at Austin) in biased discharges with non-monotonic plasma flows.

Dirac fermions in strong electric field and quantum transport in graphene

Our previous results on the nonperturbative calculations of the mean current and of the energy-momentum tensor in QED with the T-constant electric field are generalized to arbitrary dimensions. The renormalized mean values are found, and the vacuum polarization contributions and particle creation contributions to these mean values are isolated in the large T limit; we also relate the vacuum polarization contributions to the one-loop effective Euler-Heisenberg Lagrangian. Peculiarities in odd dimensions are considered in detail. We adapt general results obtained in 2 + 1 dimensions to the conditions which are realized in the Dirac model for graphene. We study the quantum electronic and energy transport in the graphene at low carrier density and low temperatures when quantum interference effects are important. Our description of the quantum transport in the graphene is based on the so-called generalized Furry picture in QED where the strong external field is taken into account nonperturbatively; this approach is not restricted to a semiclassical approximation for carriers and does not use any statistical assumptions inherent in the Boltzmann transport theory. In addition, we consider the evolution of the mean electromagnetic field in the graphene, taking into account the backreaction of the matter field to the applied external field. We find solutions of the corresponding Dirac-Maxwell set of equations and with their help we calculate the effective mean electromagnetic field and effective mean values of the current and the energy-momentum tensor. The nonlinear and linear I-V characteristics experimentally observed in both low-and high-mobility graphene samples are quite well explained in the framework of the proposed approach, their peculiarities being essentially due to the carrier creation from the vacuum by the applied electric field. DOI: 10.1103/PhysRevD.86.125022

On Fermat's principle for causal curves in time oriented Finsler spacetimes

In this work, a version of Fermat's principle for causal curves with the same energy in time orientable Finsler spacetimes is proved. We calculate the second variation of the time arrival functional along a geodesic in terms of the index form associated with the Finsler spacetime Lagrangian. Then the character of the critical points of the time arrival functional is investigated and a Morse index theorem in the context of Finsler spacetime is presented. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4765066]

Evaluating bound-constrained minimization software

Bound-constrained minimization is a subject of active research. To assess the performance of existent solvers, numerical evaluations and comparisons are carried on. Arbitrary decisions that may have a crucial effect on the conclusions of numerical experiments are highlighted in the present work. As a result, a detailed evaluation based on performance profiles is applied to the comparison of bound-constrained minimization solvers. Extensive numerical results are presented and analyzed.

Nontwist symplectic maps in tokamaks

We review symplectic nontwist maps that we have introduced to describe Lagrangian transport properties in magnetically confined plasmas in tokamaks. These nontwist maps are suitable to describe the formation and destruction of transport barriers in the shearless region (i.e., near the curve where the twist condition does not hold). The maps can be used to investigate two kinds of problems in plasmas with non-monotonic field profiles: the first is the chaotic magnetic field line transport in plasmas with external resonant perturbations. The second problem is the chaotic particle drift motion caused by electrostatic drift waves. The presented analytical maps, derived from plasma models with equilibrium field profiles and control parameters that are commonly measured in plasma discharges, can be used to investigate long-term transport properties. (C) 2011 Elsevier B.V. All rights reserved.

A marker-and-cell approach to free surface 2-D multiphase flows

This work describes a methodology to simulate free surface incompressible multiphase flows. This novel methodology allows the simulation of multiphase flows with an arbitrary number of phases, each of them having different densities and viscosities. Surface and interfacial tension effects are also included. The numerical technique is based on the GENSMAC front-tracking method. The velocity field is computed using a finite-difference discretization of a modification of the NavierStokes equations. These equations together with the continuity equation are solved for the two-dimensional multiphase flows, with different densities and viscosities in the different phases. The governing equations are solved on a regular Eulerian grid, and a Lagrangian mesh is employed to track free surfaces and interfaces. The method is validated by comparing numerical with analytic results for a number of simple problems; it was also employed to simulate complex problems for which no analytic solutions are available. The method presented in this paper has been shown to be robust and computationally efficient. Copyright (c) 2012 John Wiley & Sons, Ltd.

Holographic field theory models of dark energy in interaction with dark matter

We discuss two Lagrangian interacting dark energy models in the context of the holographic principle. The potentials of the interacting fields are constructed. The models are compared with CMB distance information, baryonic acoustic oscillations, lookback time and the Constitution supernovae sample. For both models, the results are consistent with a nonvanishing interaction in the dark sector of the Universe and the sign of coupling is consistent with dark energy decaying into dark matter, alleviating the coincidence problem-with more than 3 standard deviations of confidence for one of them. However, this is because the noninteracting holographic dark energy model is a bad fit to the combination of data sets used in this work as compared to the cosmological constant with cold dark matter model, so that one needs to introduce the interaction in order to improve this model.

Impact of the Manaus urban plume on trace gas mixing ratios near the surface in the Amazon Basin: Implications for the NO-NO2-O-3 photostationary state and peroxy radical levels

We measured the mixing ratios of NO, NO2, O-3, and volatile organic carbon as well as the aerosol light-scattering coefficient on a boat platform cruising on rivers downwind of the city of Manaus (Amazonas State, Brazil) in July 2001 (Large-Scale Biosphere-Atmosphere Experiment in Amazonia-Cooperative LBA Airborne Regional Experiment-2001). The dispersion and impact of the Manaus plume was investigated by a combined analysis of ground-based (boat platform) and airborne trace gas and aerosol measurements as well as by meteorological measurements complemented by dispersion calculations (Hybrid Single-Particle Lagrangian Integrated Trajectory model). For the cases with the least anthropogenic influence (including a location in a so far unexplored region similar to 150 km west of Manaus on the Rio Manacapuru), the aerosol scattering coefficient, sigma(s), was below 11 Mm(-1), NOx mixing ratios remained below 0.6 ppb, daytime O-3 mixing ratios were mostly below 20 ppb and maximal isoprene mixing ratios were about 3 ppb in the afternoon. The photostationary state (PSS) was not established for these cases, as indicated by values of the Leighton ratio, Phi, well above unity. Due to the influence of river breeze systems and other thermally driven mesoscale circulations, a change of the synoptic wind direction from east-northeast to south-southeast in the afternoon often caused a substantial increase of ss and trace gas mixing ratios (about threefold for sigma(s), fivefold for NOx, and twofold for O-3), which was associated with the arrival of the Manaus pollution plume at the boat location. The ratio F reached unity within its uncertainty range at NOx mixing ratios of about 3 ppb, indicating "steady-state" conditions in cases when radiation variations, dry deposition, emissions, and reactions mostly involving peroxy radicals (XO2) played a minor role. The median midday/afternoon XO2 mixing ratios estimated using the PSS method range from 90 to 120 parts per trillion (ppt) for the remote cases (sigma(s) < 11 Mm(-1) and NOx < 0.6 ppb), while for the polluted cases our estimates are 15 to 60 ppt. These values are within the range of XO2 estimated by an atmospheric chemistry box model (Chemistry As A Box model Application-Module Efficiently Calculating the Chemistry of the Atmosphere (CAABA/MECCA)-3.0).

Vortices in the extended Skyrme-Faddeev model

We construct analytical and numerical vortex solutions for an extended Skyrme-Faddeev model in a (3 + 1) dimensional Minkowski space-time. The extension is obtained by adding to the Lagrangian a quartic term, which is the square of the kinetic term, and a potential which breaks the SO(3) symmetry down to SO(2). The construction makes use of an ansatz, invariant under the joint action of the internal SO(2) and three commuting U(1) subgroups of the Poincare group, and which reduces the equations of motion to an ordinary differential equation for a profile function depending on the distance to the x(3) axis. The vortices have finite energy per unit length, and have waves propagating along them with the speed of light. The analytical vortices are obtained for a special choice of potentials, and the numerical ones are constructed using the successive over relaxation method for more general potentials. The spectrum of solutions is analyzed in detail, especially its dependence upon special combinations of coupling constants.

TWO NEW WEAK CONSTRAINT QUALIFICATIONS AND APPLICATIONS

We present two new constraint qualifications (CQs) that are weaker than the recently introduced relaxed constant positive linear dependence (RCPLD) CQ. RCPLD is based on the assumption that many subsets of the gradients of the active constraints preserve positive linear dependence locally. A major open question was to identify the exact set of gradients whose properties had to be preserved locally and that would still work as a CQ. This is done in the first new CQ, which we call the constant rank of the subspace component (CRSC) CQ. This new CQ also preserves many of the good properties of RCPLD, such as local stability and the validity of an error bound. We also introduce an even weaker CQ, called the constant positive generator (CPG), which can replace RCPLD in the analysis of the global convergence of algorithms. We close this work by extending convergence results of algorithms belonging to all the main classes of nonlinear optimization methods: sequential quadratic programming, augmented Lagrangians, interior point algorithms, and inexact restoration.

New Optimization Algorithms for Structural Reliability Analysis

Solution of structural reliability problems by the First Order method require optimization algorithms to find the smallest distance between a limit state function and the origin of standard Gaussian space. The Hassofer-Lind-Rackwitz-Fiessler (HLRF) algorithm, developed specifically for this purpose, has been shown to be efficient but not robust, as it fails to converge for a significant number of problems. On the other hand, recent developments in general (augmented Lagrangian) optimization techniques have not been tested in aplication to structural reliability problems. In the present article, three new optimization algorithms for structural reliability analysis are presented. One algorithm is based on the HLRF, but uses a new differentiable merit function with Wolfe conditions to select step length in linear search. It is shown in the article that, under certain assumptions, the proposed algorithm generates a sequence that converges to the local minimizer of the problem. Two new augmented Lagrangian methods are also presented, which use quadratic penalties to solve nonlinear problems with equality constraints. Performance and robustness of the new algorithms is compared to the classic augmented Lagrangian method, to HLRF and to the improved HLRF (iHLRF) algorithms, in the solution of 25 benchmark problems from the literature. The new proposed HLRF algorithm is shown to be more robust than HLRF or iHLRF, and as efficient as the iHLRF algorithm. The two augmented Lagrangian methods proposed herein are shown to be more robust and more efficient than the classical augmented Lagrangian method.

Numerical assessment of stability of interface discontinuous finite element pressure spaces

The stability of two recently developed pressure spaces has been assessed numerically: The space proposed by Ausas et al. [R.F. Ausas, F.S. Sousa, G.C. Buscaglia, An improved finite element space for discontinuous pressures, Comput. Methods Appl. Mech. Engrg. 199 (2010) 1019-1031], which is capable of representing discontinuous pressures, and the space proposed by Coppola-Owen and Codina [A.H. Coppola-Owen, R. Codina, Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions, Int. J. Numer. Methods Fluids, 49 (2005) 1287-1304], which can represent discontinuities in pressure gradients. We assess the stability of these spaces by numerically computing the inf-sup constants of several meshes. The inf-sup constant results as the solution of a generalized eigenvalue problems. Both spaces are in this way confirmed to be stable in their original form. An application of the same numerical assessment tool to the stabilized equal-order P-1/P-1 formulation is then reported. An interesting finding is that the stabilization coefficient can be safely set to zero in an arbitrary band of elements without compromising the formulation's stability. An analogous result is also reported for the mini-element P-1(+)/P-1 when the velocity bubbles are removed in an arbitrary band of elements. (C) 2012 Elsevier B.V. All rights reserved.

Self-bound interacting QCD matter in compact stars

The quark gluon plasma (QGP) at zero temperature and high baryon number is a system that may be present inside compact stars. It is quite possible that this cold QGP shares some relevant features with the hot QGP observed in heavy ion collisions, being also a strongly interacting system. In a previous work we have derived from the QCD Lagrangian an equation of state (EOS) for the cold QGP, which can be considered an improved version of the MIT bag-model EOS. Compared to the latter, our EOS reaches higher values of the pressure at comparable baryon densities. This feature is due to perturbative corrections and also to nonperturbative effects. Here we apply this EOS to the study of neutron stars, discussing the absolute stability of quark matter and computing the mass-radius relation for self-bound (strange) stars. The maximum masses of the sequences exceed two solar masses, in agreement with the recently measured values of the mass of the pulsar PSR J1614-2230, and the corresponding radii of around 10-11 km.

Numerical simulation of an ethanol turbulent spray flame with RANS and diffusion combustion model

A detailed numerical simulation of ethanol turbulent spray combustion on a rounded jet flame is pre- sented in this article. The focus is to propose a robust mathematical model with relatively low complexity sub- models to reproduce the main characteristics of the cou- pling between both phases, such as the turbulence modulation, turbulent droplets dissipation, and evaporative cooling effect. A RANS turbulent model is implemented. Special features of the model include an Eulerian– Lagrangian procedure under a fully two-way coupling and a modified flame sheet model with a joint mixture fraction– enthalpy b -PDF. Reasonable agreement between measured and computed mean profiles of temperature of the gas phase and droplet size distributions is achieved. Deviations found between measured and predicted mean velocity profiles are attributed to the turbulent combustion modeling adopted