Shape optimization of pressurized air bearings

Use of externally pressurized air bearings allows for the design of mechanical systems requiring extreme precision in positioning. One application is the fine control for the positioning of mirrors in large-scale optical telescopes. Other examples come from applications in robotics and computer hard-drive manufacturing. Pressurized bearings maintain a finite separation between mechanical components by virtue of the presence of a pressurized flow of air through the gap between the components. An everyday example is an air hockey table, where a puck is levitated above the table by an array of vertical jets of air. Using pressurized bearings there is no contact between “moving parts” and hence there is no friction and no wear of sensitive components. This workshop project is focused on the problem of designing optimal static air bearings subject to given engineering constraints. Recent numerical computations of this problem, done at IBM by Robert and Hendriks, suggest that near-optimal designs can have unexpected complicated and intricate structures. We will use analytical approaches to shed some light on this situation and to offer some guides for the design process. In Section 2 the design problem is stated and formulated as an optimization problem for an elliptic boundary value problem. In Section 3 the general problem is specialized to bearings with rectangular bases. Section 4 addresses the solutions of this problem that can be obtained using variational formulations of the problem. Analysis showing the sensitive dependence to perturbations (in numerical computations or manufacturing constraints) of near-optimal designs is given in Section 5. In Section 6, a restricted class of “groove network” designs motivated by the original results of Robert and Hendriks is examined. Finally, in Section 7, we consider the design problem for circular axisymmetric air bearings

Comparison of Perron and Floquet eigenvalues in age structured cell division cycle models

We study the growth rate of a cell population that follows an age-structured PDE with time-periodic coefficients. Our motivation comes from the comparison between experimental tumor growth curves in mice endowed with intact or disrupted circadian clocks, known to exert their influence on the cell division cycle. We compare the growth rate of the model controlled by a time-periodic control on its coefficients with the growth rate of stationary models of the same nature, but with averaged coefficients. We firstly derive a delay differential equation which allows us to prove several inequalities and equalities on the growth rates. We also discuss about the necessity to take into account the structure of the cell division cycle for chronotherapy modeling. Numerical simulations illustrate the results.Comment: 26 page

Noise reconstruction for the inverse heat conduction problem

AbstractA new automatic procedure to numerically recover the sample root mean square norm of the data error for the linear inverse heat conduction problem (IHCP)—when this information is not readily available—is presented. Numerical results are described which illustrate the accuracy of the algorithm

Semi-proximal Mirror-Prox for Nonsmooth Composite Minimization

We propose a new first-order optimisation algorithm to solve high-dimensional non-smooth composite minimisation problems. Typical examples of such problems have an objective that decomposes into a non-smooth empirical risk part and a non-smooth regularisation penalty. The proposed algorithm, called Semi-Proximal Mirror-Prox, leverages the Fenchel-type representation of one part of the objective while handling the other part of the objective via linear minimization over the domain. The algorithm stands in contrast with more classical proximal gradient algorithms with smoothing, which require the computation of proximal operators at each iteration and can therefore be impractical for high-dimensional problems. We establish the theoretical convergence rate of Semi-Proximal Mirror-Prox, which exhibits the optimal complexity bounds, i.e. $O(1/\epsilon^2)$, for the number of calls to linear minimization oracle. We present promising experimental results showing the interest of the approach in comparison to competing methods

Lambek vs. Lambek: Functorial Vector Space Semantics and String Diagrams for Lambek Calculus

The Distributional Compositional Categorical (DisCoCat) model is a mathematical framework that provides compositional semantics for meanings of natural language sentences. It consists of a computational procedure for constructing meanings of sentences, given their grammatical structure in terms of compositional type-logic, and given the empirically derived meanings of their words. For the particular case that the meaning of words is modelled within a distributional vector space model, its experimental predictions, derived from real large scale data, have outperformed other empirically validated methods that could build vectors for a full sentence. This success can be attributed to a conceptually motivated mathematical underpinning, by integrating qualitative compositional type-logic and quantitative modelling of meaning within a category-theoretic mathematical framework. The type-logic used in the DisCoCat model is Lambek's pregroup grammar. Pregroup types form a posetal compact closed category, which can be passed, in a functorial manner, on to the compact closed structure of vector spaces, linear maps and tensor product. The diagrammatic versions of the equational reasoning in compact closed categories can be interpreted as the flow of word meanings within sentences. Pregroups simplify Lambek's previous type-logic, the Lambek calculus, which has been extensively used to formalise and reason about various linguistic phenomena. The apparent reliance of the DisCoCat on pregroups has been seen as a shortcoming. This paper addresses this concern, by pointing out that one may as well realise a functorial passage from the original type-logic of Lambek, a monoidal bi-closed category, to vector spaces, or to any other model of meaning organised within a monoidal bi-closed category. The corresponding string diagram calculus, due to Baez and Stay, now depicts the flow of word meanings.Comment: 29 pages, pending publication in Annals of Pure and Applied Logi

Bayesian nonparametric multivariate convex regression

In many applications, such as economics, operations research and reinforcement learning, one often needs to estimate a multivariate regression function f subject to a convexity constraint. For example, in sequential decision processes the value of a state under optimal subsequent decisions may be known to be convex or concave. We propose a new Bayesian nonparametric multivariate approach based on characterizing the unknown regression function as the max of a random collection of unknown hyperplanes. This specification induces a prior with large support in a Kullback-Leibler sense on the space of convex functions, while also leading to strong posterior consistency. Although we assume that f is defined over R^p, we show that this model has a convergence rate of log(n)^{-1} n^{-1/(d+2)} under the empirical L2 norm when f actually maps a d dimensional linear subspace to R. We design an efficient reversible jump MCMC algorithm for posterior computation and demonstrate the methods through application to value function approximation

Convolutional Dictionary Learning: Acceleration and Convergence

Convolutional dictionary learning (CDL or sparsifying CDL) has many applications in image processing and computer vision. There has been growing interest in developing efficient algorithms for CDL, mostly relying on the augmented Lagrangian (AL) method or the variant alternating direction method of multipliers (ADMM). When their parameters are properly tuned, AL methods have shown fast convergence in CDL. However, the parameter tuning process is not trivial due to its data dependence and, in practice, the convergence of AL methods depends on the AL parameters for nonconvex CDL problems. To moderate these problems, this paper proposes a new practically feasible and convergent Block Proximal Gradient method using a Majorizer (BPG-M) for CDL. The BPG-M-based CDL is investigated with different block updating schemes and majorization matrix designs, and further accelerated by incorporating some momentum coefficient formulas and restarting techniques. All of the methods investigated incorporate a boundary artifacts removal (or, more generally, sampling) operator in the learning model. Numerical experiments show that, without needing any parameter tuning process, the proposed BPG-M approach converges more stably to desirable solutions of lower objective values than the existing state-of-the-art ADMM algorithm and its memory-efficient variant do. Compared to the ADMM approaches, the BPG-M method using a multi-block updating scheme is particularly useful in single-threaded CDL algorithm handling large datasets, due to its lower memory requirement and no polynomial computational complexity. Image denoising experiments show that, for relatively strong additive white Gaussian noise, the filters learned by BPG-M-based CDL outperform those trained by the ADMM approach.Comment: 21 pages, 7 figures, submitted to IEEE Transactions on Image Processin