Search

The aim of this work is to analyse a shape optimization problem in a mechanical friction context. Precisely we perform a shape sensitivity analysis of a Tresca friction problem, that is, a boundary value problem involving the usual linear elasticity equations together with the (nonsmooth) Tresca friction law on a part of the boundary. We prove that the solution to the Tresca friction problem admits a directional shape derivative which moreover coincides with the solution to a boundary value problem involving tangential Signorini's unilateral conditions. Then an explicit expression of the shape gradient of the Tresca energy functional is provided (which allows us to provide numerical simulations illustrating our theoretical results). Our methodology is not based on any regularization procedure, but rather on the twice epi-differentiability of the (nonsmooth) Tresca friction functional which is analyzed thanks to a change of variables which is well-suited in the two-dimensional case. The obstruction in the higher-dimensional case is discussed., Comment: 30 pages. arXiv admin note: text overlap with arXiv:2410.11750, arXiv:2410.12315

This paper investigates, without any regularization procedure, the sensitivity analysis of a mechanical friction problem involving the (nonsmooth) Tresca friction law in the linear elastic model. To this aim a recent methodology based on advanced tools from convex and variational analyses is used. Precisely we express the solution to the so-called Tresca friction problem thanks to the proximal operator associated with the corresponding Tresca friction functional. Then, using an extended version of twice epi-differentiability, we prove the differentiability of the solution to the parameterized Tresca friction problem, characterizing its derivative as the solution to a boundary value problem involving tangential Signorini's unilateral conditions. Finally our result is used to investigate and numerically solve an optimal control problem associated with the Tresca friction model., Comment: 37 pages. arXiv admin note: text overlap with arXiv:2410.11760, arXiv:2410.12315

This paper investigates the sensitivity analysis of a scalar mechanical contact problem described by a boundary value problem involving the Tresca's friction law. The sensitivity analysis is performed with respect to right-hand source and boundary terms perturbations. In particular the friction threshold involved in the Tresca's friction law is perturbed, which constitutes the main novelty of the present work with respect to the existing literature. Hence we introduce a parameterized Tresca friction problem and its solution is characterized by using the proximal operator associated with the corresponding perturbed nonsmooth convex Tresca friction functional. Then, by invoking the extended notion of twice epi-differentiability depending on a parameter, we prove the differentiability of the solution to the parameterized Tresca friction problem, characterizing its derivative as the solution to a boundary value problem involving Signorini unilateral conditions. Finally numerical simulations are provided in order to illustrate our main result., Comment: 30 pages

This paper investigates, without any regularization or penalization procedure, a shape optimization problem involving a simplified friction phenomena modeled by a scalar Tresca friction law. Precisely, using tools from convex and variational analysis such as proximal operators and the notion of twice epi-differentiability, we prove that the solution to a scalar Tresca friction problem admits a directional derivative with respect to the shape which moreover coincides with the solution to a boundary value problem involving Signorini-type unilateral conditions. Then we explicitly characterize the shape gradient of the corresponding energy functional and we exhibit a descent direction. Finally numerical simulations are performed to solve the corresponding energy minimization problem under a volume constraint which shows the applicability., Comment: 30 pages

In this work is considered a spectral problem, involving a second order term on the domain boundary: the Laplace-Beltrami operator. A variational formulation is presented, leading to a finite element discretization. For the Laplace-Beltrami operator to make sense on the boundary, the domain is smooth: consequently the computational domain (classically a polygonal domain) will not match the physical one. Thus, the physical domain is discretized using high order curved meshes so as to reduce the \textit{geometric error}. The \textit{lift operator}, which is aimed to transform a function defined on the mesh domain into a function defined on the physical one, is recalled. This \textit{lift} is a key ingredient in estimating errors on eigenvalues and eigenfunctions. A bootstrap method is used to prove the error estimates, which are expressed both in terms of \textit{finite element approximation error} and of \textit{geometric error}, respectively associated to the finite element degree $k\ge 1$ and to the mesh order~$r\ge 1$. Numerical experiments are led on various smooth domains in 2D and 3D, which allow us to validate the presented theoretical results.

This paper is devoted to the understanding of regularisation process in the shape optimization approach to the so-called Dirichlet inverse obstacle problem for elliptic operators. More precisely, we study two different regularisations of the very classical shape optimization approach consisting in minimizing a mismatched functional. The first one is an implicit regularisation when working in the class of inclusion having a uniform $\varepsilon$-cone property, a natural class in shape optimization. As this regularity is not trivial to guarantee numerically, we discuss the regularisation by perimeter penalization. We show that this second regularisation provides a stability gain in the minimization process.

In this work is considered an elliptic problem, referred to as the Ventcel problem, involvinga second order term on the domain boundary (the Laplace-Beltrami operator). A variationalformulation of the Ventcel problem is studied, leading to a finite element discretization. Thefocus is on the construction of high order curved meshes for the discretization of the physicaldomain and on the definition of the lift operator, which is aimed to transform a functiondefined on the mesh domain into a function defined on the physical one. This lift is definedin a way as to satisfy adapted properties on the boundary, relatively to the trace operator.The Ventcel problem approximation is investigated both in terms of geometrical error and offinite element approximation error. Error estimates are obtained both in terms of the meshorder r $\ge$ 1 and to the finite element degree k $\ge$ 1, whereas such estimates usually have beenconsidered in the isoparametric case so far, involving a single parameter k = r. The numericalexperiments we led, both in dimension 2 and 3, allow us to validate the results obtained andproved on the a priori error estimates depending on the two parameters k and r. A numericalcomparison is made between the errors using the former lift definition and the lift defined inthis work establishing an improvement in the convergence rate of the error in the latter case.

In this work is provided a numerical study of a diffusion problem involving a second order term on the domain boundary (the Laplace-Beltrami operator) referred to as the \textit{Ventcel problem}.A variational formulation of the Ventcel problem is studied, leading to a finite element discretization.The focus is on the resort to high order curved meshes for the discretization of the physical domain.The computational errors are investigated both in terms of geometrical error and of finite element approximation error, respectively associated to the mesh degree $r\ge 1$ and to the finite element degree $k\ge 1$. The numerical experiments we led allow us to formulate a conjecture on the \textit{a priori} error estimates depending on the two parameters $r$ and $k$. In addition, these error estimates rely on the definition of a functional \textit{lift} with adapted properties on the boundary to move numerical solutions defined on the computational domain to the physical one.

We present CosmoHub (https://cosmohub.pic.es), a web application based on Hadoop to perform interactive exploration and distribution of massive cosmological datasets. Recent Cosmology seeks to unveil the nature of both dark matter and dark energy mapping the large-scale structure of the Universe, through the analysis of massive amounts of astronomical data, progressively increasing during the last (and future) decades with the digitization and automation of the experimental techniques. CosmoHub, hosted and developed at the Port d'Informaci\'o Cient\'ifica (PIC), provides support to a worldwide community of scientists, without requiring the end user to know any Structured Query Language (SQL). It is serving data of several large international collaborations such as the Euclid space mission, the Dark Energy Survey (DES), the Physics of the Accelerating Universe Survey (PAUS) and the Marenostrum Institut de Ci\`encies de l'Espai (MICE) numerical simulations. While originally developed as a PostgreSQL relational database web frontend, this work describes the current version of CosmoHub, built on top of Apache Hive, which facilitates scalable reading, writing and managing huge datasets. As CosmoHub's datasets are seldomly modified, Hive it is a better fit. Over 60 TiB of catalogued information and $50 \times 10^9$ astronomical objects can be interactively explored using an integrated visualization tool which includes 1D histogram and 2D heatmap plots. In our current implementation, online exploration of datasets of $10^9$ objects can be done in a timescale of tens of seconds. Users can also download customized subsets of data in standard formats generated in few minutes.

We are interested in geometric approximation by parameterization of two-dimensional multiple-component shapes, in particular when the number of components is a priori unknown. Starting a standard method based on successive shape deformations with a one-component initial shape in order to approximate a multiple-component target shape usually leads the deformation flow to make the boundary evolve until it surrounds all the components of the target shape. This classical phenomenon tends to create double points on the boundary of the approximated shape. In order to improve the approximation of multiple-component shapes (without any knowledge on the number of components in advance), we use in this paper a piecewise B\'ezier parameterization and we consider two procedures called intersecting control polygons detection and flip procedure. The first one allows to prevent potential collisions between two parts of the boundary of the approximated shape, and the second one permits to change its topology by dividing a one-component shape into a two-component shape. For an experimental purpose, we include these two processes in a basic geometrical shape optimization algorithm and test it on the classical inverse obstacle problem. This new approach allows to obtain a numerical approximation of the unknown inclusion, detecting both the topology (i.e. the number of connected components) and the shape of the obstacle. Several numerical simulations are performed.

The main result of the present theoretical paper is an original decomposition formula for the proximal operator of the sum of two proper, lower semicontinuous and convex functions $f$ and $g$. For this purpose, we introduce a new operator, called $f$-proximal operator of $g$ and denoted by $\mathrm{prox}^f_g$, that generalizes the classical notion. Then we prove the decomposition formula $\mathrm{prox}_{f+g} = \mathrm{prox}_f \circ \mathrm{prox}^f_g$. After collecting several properties and characterizations of $\mathrm{prox}^f_g$, we prove that it coincides with the fixed points of a generalized version of the classical Douglas-Rachford operator. This relationship is used for the construction of a weakly convergent algorithm that computes numerically this new operator $\mathrm{prox}^f_g$, and thus, from the decomposition formula, allows to compute numerically $\mathrm{prox}_{f+g}$. It turns out that this algorithm was already considered and implemented in previous works, showing that $\mathrm{prox}^f_g$ is already present (in a hidden form) and useful for numerical purposes in the existing literature. However, to the best of our knowledge, it has never been explicitly expressed in a closed formula and neither been deeply studied from a theoretical point of view. The present paper contributes to fill this gap in the literature. Finally we give an illustration of the usefulness of the decomposition formula in the context of sensitivity analysis of linear variational inequalities of second kind in a Hilbert space., Comment: To appear in Journal of Convex Analysis, 2018

In this paper, we investigate an optimal design problem motivated by some issues arising in population dynamics. In a nutshell, we aim at determining the optimal shape of a region occupied by resources for maximizing the survival ability of a species in a given box and we consider the general case of Robin boundary conditions on its boundary. Mathematically, this issue can be modeled with the help of an extremal indefinite weight linear eigenvalue problem. The optimal spatial arrangement is obtained by minimizing the positive principal eigenvalue with respect to the weight, under a L 1 constraint standing for limitation of the total amount of resources. The specificity of such a problem rests upon the presence of nonlinear functions of the weight both in the numerator and denominator of the Rayleigh quotient. By using adapted symmetrization procedures, a well-chosen change of variable, as well as necessary optimality conditions, we completely solve this optimization problem in the unidimensional case by showing first that every minimizer is unimodal and bang-bang. This leads to investigate a finite dimensional optimization problem. This allows to show in particular that every minimizer is (up to additive constants) the characteristic function of three possible domains: an interval that sticks on the boundary of the box, an interval that is symmetrically located at the middle of the box, or, for a precise value of the Robin coefficient, all intervals of a given fixed length.

In this work, we present some new Carleman inequalities for Stokes and Oseen equations with non-homogeneous boundary conditions. These estimates lead to log type stability inequalities for the problem of recovering the solution of the Stokes and Navier-Stokes equations from both boundary and distributed observations. These inequalities fit the well-known unique continuation result of Fabre and Lebeau [18]: the distributed observation only depends on interior measurement of the velocity, and the boundary observation only depends on the trace of the velocity and of the Cauchy stress tensor measurements. Finally, we present two applications for such inequalities. First, we apply these estimates to obtain stability inequalities for the inverse problem of recovering Navier or Robin boundary coefficients from boundary measurements. Next, we use these estimates to deduce the rate of convergence of two reconstruction methods of the Stokes solution from the measurement of Cauchy data: a quasi-reversibility method and a penalized Kohn-Vogelius method.

Rate constants for the C(3P) + CH3OH reaction have been measured in a continuous supersonic flow reactor over the range 50 K to 296 K. C(3P) was created by the in-situ pulsed laser photolysis of CBr4, a multiphoton process which also produced some C(1D), allowing us to investigate simultaneously the low temperature kinetics of the C(1D) + CH3OH reaction. C(1D) atoms were followed by an indirect chemiluminescent tracer method in the presence of excess CH3OH. C(3P) atoms were detected by the same chemiluminescence technique and also by direct vacuum ultra-violet laser induced fluorescence (VUV LIF). Secondary measurements of product H(2S) atom formation have been undertaken allowing absolute H atom yields to be obtained by comparison with a suitable reference reaction. In parallel, statistical calculations have been performed based on ab-initio calculations of the complexes, adducts and transition states (TSs) relevant to the title reaction. By comparison with the experimental H atom yields, the preferred reaction pathways could be determined, placing important constraints on the statistical calculations. The experimental and theoretical work are in excellent agreement, predicting a negative temperature dependence of the rate constant increasing from 2.2 x 10-11 cm3 molecule-1 s-1 at 296 K to 20.0 x 10-11 cm3 molecule-1 s-1 at 50 K. CH3 and HCO are found to be the major products under our experimental conditions. As this reaction is not considered in current astrochemical networks, its influence on interstellar methanol abundances is tested using a model of dense interstellar clouds., Comment: Accepted to RSC Advances

Water is an important reservoir species for oxygen in interstellar space and plays a key role in the physics of star formation through cooling by far-infrared emission. Whilst water vapour is present at high abundances in the outflows of protostars, its contribution to the chemical evolution of these regions is a minor one due to its limited low temperature reactivity in the gas-phase. Here, we performed kinetic experiments on the barrierless CH + H2O reaction in a supersonic flow reactor down to 50 K. The measured rate increases rapidly below room temperature, confirming and extending the predictions of earlier statistical calculations. The open product channels for this reaction suggest that this process could be an important gas-phase route for formaldehyde formation in protostellar envelopes., Comment: Accepted to Journal of Physical Chemistry Letters. 24 pages, 8 figures, 3 tables

Rate constants for the potentially important interstellar N(4S) + CH(X2{\Pi}r) reaction have been measured in a continuous supersonic flow reactor over the range 56 K < T < 296 K using the relative rate technique employing both the N(4S) + OH(X2{\Pi}i) and N(4S) + CN(X2{\Sigma}+) reactions as references. Excess concentrations of atomic nitrogen were produced by the microwave discharge method upstream of the Laval nozzle and CH and OH radicals were created by the in-situ pulsed laser photolysis of suitable precursor molecules. In parallel, quantum dynamics calculations of the title reaction have been performed based on accurate global potential energy surfaces for the 13A' and 13A" states of HCN and HNC, brought about through a hierarchical construction scheme. Both adiabatic potential energy surfaces are barrierless, each one having two deep potential wells suggesting that this reaction is dominated by a complex-forming mechanism. The experimental and theoretical work are inexcellent agreement, predicting a positive temperature dependence of the rate constant, in contrast to earlier experimental work at low temperature. The effects of the new low temperature rate constants on interstellar N2 formation are tested using a dense cloud model, yielding N2 abundances 10-20 % lower than previously predicted., Comment: Accepted for publication in Chemical Physics Physical Chemistry