Geometric characteristics of conics in Bézier form

Derive coordinate-free expressions for geometric characteristics of conics written in Bézier form in terms of their control points and weights.

Convective transport and stability in films of binary mixtures

Thin polymer films are increasingly used in advanced technological applications. The use of these films as coatings is often limited by their lack of stability due to their wettability properties on the substrates

A second order in time modified Lagrange-Galerkin finite element method for the incompressible Navier-Stokes equations

We introduce a second order in time modified Lagrange--Galerkin (MLG) method for the time dependent incompressible Navier--Stokes equations. The main ingredient of the new method is the scheme proposed to calculate in a more efficient manner the Galerkin projection of the functions transported along the characteristic curves of the transport operator. We present error estimates for velocity and pressure in the framework of mixed finite elements when either the mini-element or the $P2/P1$ Taylor--Hood element are used.

Lagrangian skeleton on the Gulf of Mexico oil spill

The understanding of the circulation of ocean currents, the exchange of CO2 between atmosphere and oceans, and the influence of the oceans on the distribution of heat on a global scale is key to our ability to predict and assess the future evolution of climate [1, 2]. Global climate change is affecting sea breathing through mechanisms not yet understood.

A Semi-Lagrangian Particle Level Set Finite Element Method for Interface Problems

We present a quasi-monotone semi-Lagrangian particle level set (QMSL-PLS) method for moving interfaces. The QMSL method is a blend of first order monotone and second order semi-Lagrangian methods. The QMSL-PLS method is easy to implement, efficient, and well adapted for unstructured, either simplicial or hexahedral, meshes. We prove that it is unconditionally stable in the maximum discrete norm, � · �h,∞, and the error analysis shows that when the level set solution u(t) is in the Sobolev space Wr+1,∞(D), r ≥ 0, the convergence in the maximum norm is of the form (KT/Δt)min(1,Δt � v �h,∞ /h)((1 − α)hp + hq), p = min(2, r + 1), and q = min(3, r + 1),where v is a velocity. This means that at high CFL numbers, that is, when Δt > h, the error is O( (1−α)hp+hq) Δt ), whereas at CFL numbers less than 1, the error is O((1 − α)hp−1 + hq−1)). We have tested our method with satisfactory results in benchmark problems such as the Zalesak’s slotted disk, the single vortex flow, and the rising bubble.

On a two species chemotaxis model with slow chemical diffucion.

In this paper we consider a system of three parabolic equations modeling the behavior of two biological species moving attracted by a chemical factor. The chemical substance verifies a parabolic equation with slow diffusion. The system contains second order terms in the first two equations modeling the chemotactic effects. We apply an iterative method to obtain the global existence of solutions using that the total mass of the biological species is conserved. The stability of the homogeneous steady states is studied by using an energy method. A final example is presented to illustrate the theoretical results.

Lagrangian tools and the assessment of the predictive capacity of geophysical data sets

We examine, with recently developed Lagrangian tools, altimeter data and numerical simulations obtained from the HYCOM model in the Gulf of Mexico. Our data correspond to the months just after the Deepwater Horizon oil spill in the year 2010. Our Lagrangian analysis provides a skeleton that allows the interpretation of transport routes over the ocean surface. The transport routes are further verified by the simultaneous study of the evolution of several drifters launched during those months in the Gulf of Mexico. We find that there exist Lagrangian structures that justify the dynamics of the drifters, although the agreement depends on the quality of the data. We discuss the impact of the Lagrangian tools on the assessment of the predictive capacity of these data sets.

Asymptotic analysis of vertical geothermal boreholes in the limit of slowly varyng heat injection rates

Theoretical models for the thermal response of vertical geothermal boreholes often assume that the characteristic time of variation of the heat injection rate is much larger than the characteristic diffusion time across the borehole. In this case, heat transfer inside the borehole and in its immediate surroundings is quasi-steady in the first approximation, while unsteady effects enter only in the far field. Previous studies have exploited this disparity of time scales, incorporating approximate matching conditions to couple the near-borehole region with the outer unsteady temperatura field. In the present work matched asymptotic expansion techniques are used to analyze the heat transfer problem, delivering a rigorous derivation of the true matching condition between the two regions and of the correct definition of the network of thermal resistances that represents the quasi-steady solution near the borehole. Additionally, an apparent temperature due to the unsteady far field is identified that needs to be taken into account by the near-borehole region for the correct computation of the heat injection rate. This temperature differs from the usual mean borehole temperature employed in the literatura.

A fast, high resolution, second-order central scheme for incompressible flows

A high resolution, second-order central difference method for incompressible flows is presented. The method is based on a recent second-order extension of the classic Lax–Friedrichs scheme introduced for hyperbolic conservation laws (Nessyahu H. & Tadmor E. (1990) J. Comp. Physics. 87, 408-463; Jiang G.-S. & Tadmor E. (1996) UCLA CAM Report 96-36, SIAM J. Sci. Comput., in press) and augmented by a new discrete Hodge projection. The projection is exact, yet the discrete Laplacian operator retains a compact stencil. The scheme is fast, easy to implement, and readily generalizable. Its performance was tested on the standard periodic double shear-layer problem; no spurious vorticity patterns appear when the flow is underresolved. A short discussion of numerical boundary conditions is also given, along with a numerical example.

“Coarse” stability and bifurcation analysis using time-steppers: A reaction-diffusion example

Evolutionary, pattern forming partial differential equations (PDEs) are often derived as limiting descriptions of microscopic, kinetic theory-based models of molecular processes (e.g., reaction and diffusion). The PDE dynamic behavior can be probed through direct simulation (time integration) or, more systematically, through stability/bifurcation calculations; time-stepper-based approaches, like the Recursive Projection Method [Shroff, G. M. & Keller, H. B. (1993) SIAM J. Numer. Anal. 30, 1099–1120] provide an attractive framework for the latter. We demonstrate an adaptation of this approach that allows for a direct, effective (“coarse”) bifurcation analysis of microscopic, kinetic-based models; this is illustrated through a comparative study of the FitzHugh-Nagumo PDE and of a corresponding Lattice–Boltzmann model.

Quantitative stability of linear infinite inequality systems under block perturbations with applications to convex systems

The original motivation for this paper was to provide an efficient quantitative analysis of convex infinite (or semi-infinite) inequality systems whose decision variables run over general infinite-dimensional (resp. finite-dimensional) Banach spaces and that are indexed by an arbitrary fixed set J. Parameter perturbations on the right-hand side of the inequalities are required to be merely bounded, and thus the natural parameter space is l ∞(J). Our basic strategy consists of linearizing the parameterized convex system via splitting convex inequalities into linear ones by using the Fenchel–Legendre conjugate. This approach yields that arbitrary bounded right-hand side perturbations of the convex system turn on constant-by-blocks perturbations in the linearized system. Based on advanced variational analysis, we derive a precise formula for computing the exact Lipschitzian bound of the feasible solution map of block-perturbed linear systems, which involves only the system’s data, and then show that this exact bound agrees with the coderivative norm of the aforementioned mapping. In this way we extend to the convex setting the results of Cánovas et al. (SIAM J. Optim. 20, 1504–1526, 2009) developed for arbitrary perturbations with no block structure in the linear framework under the boundedness assumption on the system’s coefficients. The latter boundedness assumption is removed in this paper when the decision space is reflexive. The last section provides the aimed application to the convex case.

Calmness modulus of linear semi-infinite programs

Our main goal is to compute or estimate the calmness modulus of the argmin mapping of linear semi-infinite optimization problems under canonical perturbations, i.e., perturbations of the objective function together with continuous perturbations of the right-hand side of the constraint system (with respect to an index ranging in a compact Hausdorff space). Specifically, we provide a lower bound on the calmness modulus for semi-infinite programs with unique optimal solution which turns out to be the exact modulus when the problem is finitely constrained. The relationship between the calmness of the argmin mapping and the same property for the (sub)level set mapping (with respect to the objective function), for semi-infinite programs and without requiring the uniqueness of the nominal solution, is explored, too, providing an upper bound on the calmness modulus of the argmin mapping. When confined to finitely constrained problems, we also provide a computable upper bound as it only relies on the nominal data and parameters, not involving elements in a neighborhood. Illustrative examples are provided.

From the Farkas Lemma to the Hahn–Banach Theorem

This paper provides new versions of the Farkas lemma characterizing those inequalities of the form f(x) ≥ 0 which are consequences of a composite convex inequality (S ◦ g)(x) ≤ 0 on a closed convex subset of a given locally convex topological vector space X, where f is a proper lower semicontinuous convex function defined on X, S is an extended sublinear function, and g is a vector-valued S-convex function. In parallel, associated versions of a stable Farkas lemma, considering arbitrary linear perturbations of f, are also given. These new versions of the Farkas lemma, and their corresponding stable forms, are established under the weakest constraint qualification conditions (the so-called closedness conditions), and they are actually equivalent to each other, as well as equivalent to an extended version of the so-called Hahn–Banach–Lagrange theorem, and its stable version, correspondingly. It is shown that any of them implies analytic and algebraic versions of the Hahn–Banach theorem and the Mazur–Orlicz theorem for extended sublinear functions.

Robust Solutions of MultiObjective Linear Semi-Infinite Programs under Constraint Data Uncertainty

The multiobjective optimization model studied in this paper deals with simultaneous minimization of finitely many linear functions subject to an arbitrary number of uncertain linear constraints. We first provide a radius of robust feasibility guaranteeing the feasibility of the robust counterpart under affine data parametrization. We then establish dual characterizations of robust solutions of our model that are immunized against data uncertainty by way of characterizing corresponding solutions of robust counterpart of the model. Consequently, we present robust duality theorems relating the value of the robust model with the corresponding value of its dual problem.

II. Academic notebook, [ca. 1745]

This paper notebook contains undated and unattributed handwritten summaries and notes from published texts on world history, religion, and natural history, and was presumably kept by William Sever. The first page notes the death of poet Alexander Pope in 1744, constituting the latest date referenced in the notebook. The entries are typically short and offer condensed information on subjects including the histories of England, Japan and Siam; early Christian history; and natural religion.