This paper studies outflow of a light fluid from a point source, starting from an initially spherical bubble. This region of light fluid is embedded in a heavy fluid, from which it is separated by a thin interface. A gravitational force directed radially inward toward the mass source is permitted. Because the light inner fluid is pushing the heavy outer fluid, the interface between them may be unstable to small perturbations, in the Rayleigh–Taylor sense. An inviscid model of this two-layer flow is presented, and a linearized solution is developed for early times. It is argued that the inviscid solution develops a point of infinite curvature at the interface within finite time, after which the solution fails to exist. A Boussinesq viscous model is then presented as a means of quantifying the precise effects of viscosity. The interface is represented as a narrow region of large density gradient. The viscous results agree well with the inviscid theory at early times, but the curvature singularity of the inviscid theory is instead replaced by jet formation in the viscous case. This may be of relevance to underwater explosions and stellar evolution.

We consider the linear wave equation with Dirichlet boundary conditions in a bounded interval, and with a control acting on a moving point. We give sufficient conditions on the trajectory of the control in order to have the exact controllability property.

An approximation procedure for time optimal control problems for the linear wave equation is analyzed. Its asymptotic behavior is investigated and an optimality system including the maximum principle and the transversality conditions for the regularized and unregularized problems are derived.

The paper deals with deterministic optimal control problems with state constraints and non-linear dynamics. It is known for such problems that the value function is in general discontinuous and its characterization by means of a Hamilton-Jacobi equation requires some controllability assumptions involving the dynamics and the set of state constraints. Here, we first adopt the viability point of view and look at the value function as its epigraph. Then, we prove that this epigraph can always be described by an auxiliary optimal control problem free of state constraints, and for which the value function is Lipschitz continuous and can be characterized, without any additional assumptions, as the unique viscosity solution of a Hamilton-Jacobi equation. The idea introduced in this paper bypasses the regularity issues on the value function of the constrained control problem and leads to a constructive way to compute its epigraph by a large panel of numerical schemes. Our approach can be extended to more general control problems. We study in this paper the extension to the infinite horizon problem as well as for the two-player game setting. Finally, an illustrative numerical example is given to show the relevance of the approach.

We develop a Discrete Element Method (DEM) for elastodynamics using polyhedral elements. We show that for a given choice of forces and torques, we recover the equations of linear elastodynamics in small deformations. Furthermore, the torques and forces derive from a potential energy, and thus the global equation is an Hamiltonian dynamics. The use of an explicit symplectic time integration scheme allows us to recover conservation of energy, and thus stability over long time simulations. These theoretical results are illustrated by numerical simulations of test cases involving large displacements.

We derive and analyze adaptive solvers for boundary value problems in which thedifferential operator depends affinely on a sequence of parameters. These methods convergeuniformly in the parameters and provide an upper bound for the maximal error. Numericalcomputations indicate that they are more efficient than similar methods that control theerror in a mean square sense.

This paper is concerned with the numerical approximation of mean curvature flow t → Ω(t) satisfying an additional inclusion-exclusion constraint Ω1 ⊂ Ω(t) ⊂ Ω2. Classical phase field model to approximate these evolving interfaces consists in solving the Allen-Cahn equation with Dirichlet boundary conditions. In this work, we introduce a new phase field model, which can be viewed as an Allen Cahn equation with a penalized double well potential. We first justify this method by a Γ-convergence result and then show some numerical comparisons of these two different models.

In this paper we deal with the null controllability problem for the heat equation with amemory term by means of boundary controls. For each positive final time Tand when the control is acting on the whole boundary, we prove that there exists a set ofinitial conditions such that the null controllability property fails.

Spatial data sets can be analysed by counting the number of objects in equally sized bins. The bin counts are related to the Pólya urn process, where coloured balls (for example, white or black) are removed from the urn at random. If there are insufficient white or black balls for the prescribed number of trials, the Pólya urn process becomes untenable. In this case, we modify the Pólya urn process so that it continues to describe the removal of volume within a spatial distribution of objects. We determine when the standard formula for the variance of the standard Pólya distribution gives a good approximation to the true variance. The variance quantifies an index for assessing whether a spatial point data set is at its most randomly distributed state, called the complete spatial randomness (CSR) state. If the bin size is an order of magnitude larger than the size of the objects, then the standard formula for the CSR limit is indicative of when the CSR state has been attained. For the special case when the object size divides the bin size, the standard formula is in fact exact.

In Carnot groups of step  ≤ 3, all subriemannian geodesics are proved to be normal. Theproof is based on a reduction argument and the Goh condition for minimality of singularcurves. The Goh condition is deduced from a reformulation and a calculus of the end-pointmapping which boils down to the graded structures of Carnot groups.