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]

Homogeneous structures and rigidity of isoparametric submanifolds in Hilbert space

We study isoparametric submanifolds of rank at least two in a separable Hilbert space, which are known to be homogeneous by the main result in [E. Heintze and X. Liu, Ann. of Math. (2), 149 (1999), 149-181], and with such a submanifold M and a point x in M we associate a canonical homogeneous structure I" (x) (a certain bilinear map defined on a subspace of T (x) M x T (x) M). We prove that I" (x) , together with the second fundamental form alpha (x) , encodes all the information about M, and we deduce from this the rigidity result that M is completely determined by alpha (x) and (Delta alpha) (x) , thereby making such submanifolds accessible to classification. As an essential step, we show that the one-parameter groups of isometries constructed in [E. Heintze and X. Liu, Ann. of Math. (2), 149 (1999), 149-181] to prove their homogeneity induce smooth and hence everywhere defined Killing fields, implying the continuity of I" (this result also seems to close a gap in [U. Christ, J. Differential Geom., 62 (2002), 1-15]). Here an important tool is the introduction of affine root systems of isoparametric submanifolds.

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.

Polar foliations and isoparametric maps

A singular Riemannian foliation F on a complete Riemannian manifold M is called a polar foliation if, for each regular point p, there is an immersed submanifold Sigma, called section, that passes through p and that meets all the leaves and always perpendicularly. A typical example of a polar foliation is the partition of M into the orbits of a polar action, i.e., an isometric action with sections. In this article we prove that the leaves of H : M -> Sigma, coincide with the level sets of a smooth map H: M -> Sigma, if M is simply connected. In particular, the orbits of a polar action on a simply connected space are level sets of an isoparametric map. This result extends previous results due to the author and Gorodski, Heintze, Liu and Olmos, Carter and West, and Terng.

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

On the lifecycle of aerosol particles : Sources and dispersion over Scandinavia

Aerosol particles are likely important contributors to our future climate. Further, during recent years, effects on human health arising from emissions of particulate material have gained increasing attention. In order to quantify the effect of aerosols on both climate and human health we need to better quantify the interplay between sources and sinks of aerosol particle number and mass on large spatial scales. So far long-term, regional observations of aerosol properties have been scarce, but argued necessary in order to bring the knowledge of regional and global distribution of aerosols further. In this context, regional studies of aerosol properties and aerosol dynamics are truly important areas of investigation. This thesis is devoted to investigations of aerosol number size distribution observations performed through the course of one year encompassing observational data from five stations covering an area from southern parts of Sweden up to northern parts of Finland. This thesis tries to give a description of aerosol size distribution dynamics from both a quantitative and qualitative point of view. The thesis focuses on properties and changes in aerosol size distribution as a function of location, season, source area, transport pathways and links to various meteorological conditions. The investigations performed in this thesis show that although the basic behaviour of the aerosol number size distribution in terms of seasonal and diurnal characteristics is similar at all stations in the measurement network, the aerosol over the Nordic countries is characterised by a typically sharp gradient in aerosol number and mass. This gradient is argued to derive from geographical locations of the stations in relation to the dominant sources and transport pathways. It is clear that the source area significantly determine the aerosol size distribution properties, but it is obvious that transport condition in terms of frequency of precipitation and cloudiness in some cases even more strongly control the evolution of the number size distribution. Aerosol dynamic processes under clear sky transport are however likewise argued to be highly important. Southerly transport of marine air and northerly transport of air from continental sources is studied in detail under clear sky conditions by performing a pseudo-Lagrangian box model evaluation of the two type cases. Results from both modelling and observations suggest that nucleation events contribute to integral number increase during southerly transport of comparably clean marine air, while number depletion dominates the evolution of the size distribution during northerly transport. This difference is largely explained by different concentration of pre-existing aerosol surface associated with the two type cases. Mass is found to be accumulated in many of the individual transport cases studied. This mass increase was argued to be controlled by emission of organic compounds from the boreal forest. This puts the boreal forest in a central position for estimates of aerosol forcing on a regional scale.