Source author record

Roberto Camassa

Roberto Camassa appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2023arXiv

Critical density triplets for the arrestment of a sphere falling in a sharply stratified fluid

We study the motion of a rigid sphere falling in a two-layer stratified fluid under the action of gravity in the potential flow regime. Experiments at a moderate Reynolds number of approximately 20 to 450 indicate that a sphere with the precise critical density, higher than the bottom layer density, can display behaviors such as bounce or arrestment after crossing the interface. We experimentally demonstrate that such a critical sphere density increases linearly as the bottom fluid density increases with a fixed top fluid density. Additionally, the critical density approaches the bottom layer fluid density as the thickness of density transition layer increases. We propose an estimation of the critical density based on the potential energy. With assuming the zero layer thickness, the estimation constitutes an upper bound of the critical density with less than 0.043 relative difference within the experimental density regime 0.997 $g/cm^{3}$ $\sim $ 1.11 $g/cm^{3}$ under the zero layer thickness assumption. By matching the experimental layer thickness, we obtain a critical density estimation with less than 0.01 relative difference within the same parameter regime.

preprint2022arXiv

Numerical Algorithms for Water Waves with Background Flow over Obstacles and Topography

We present two accurate and efficient algorithms for solving the incompressible, irrotational Euler equations with a free surface in two dimensions with background flow over a periodic, multiply-connected fluid domain that includes stationary obstacles and variable bottom topography. One approach is formulated in terms of the surface velocity potential while the other evolves the vortex sheet strength. Both methods employ layer potentials in the form of periodized Cauchy integrals to compute the normal velocity of the free surface, are compatible with arbitrary parameterizations of the free surface and boundaries, and allow for circulation around each obstacle, which leads to multiple-valued velocity potentials but single-valued stream functions. We prove that the resulting second-kind Fredholm integral equations are invertible, possibly after a physically motivated finite-rank correction. In an angle-arclength setting, we show how to avoid curve reconstruction errors that are incompatible with spatial periodicity. We use the proposed methods to study gravity-capillary waves generated by flow around several elliptical obstacles above a flat or variable bottom boundary. In each case, the free surface eventually self-intersects in a splash singularity or collides with a boundary. We also show how to evaluate the velocity and pressure with spectral accuracy throughout the fluid, including near the free surface and solid boundaries. To assess the accuracy of the time evolution, we monitor energy conservation and the decay of Fourier modes and compare the numerical results of the two methods to each other. We implement several solvers for the discretized linear systems and compare their performance. The fastest approach employs a graphics processing unit (GPU) to construct the matrices and carry out iterations of the generalized minimal residual method (GMRES).

preprint2015arXiv

A class of fast geodesic shooting algorithms for template matching and its applications via the $N$-particle system of the Euler-Poincaré equations

The Euler-Poincaré (EP) equations describe the geodesic motion on the diffeomorphism group. For template matching (template deformation), the Euler-Lagrangian equation, arising from minimizing an energy function, falls into the Euler-Poincaré theory and can be recast into the EP equations. By casting the EP equations in the Lagrangian (or characteristics) form, we formulate the equations as a finite dimensional particle system. The evolution of this particle system describes the geodesic motion of landmark points on a Riemann manifold. In this paper we present a class of novel algorithms that take advantage of the structure of the particle system to achieve a fast matching process between the reference and the target templates. The strong suit of the proposed algorithms includes (1) the efficient feedback control iteration, which allows one to find the initial velocity field for driving the deformation from the reference template to the target one, (2) the use of the conical kernel in the particle system, which limits the interaction between particles and thus accelerates the convergence, and (3) the availability of the implementation of fast-multipole method for solving the particle system, which could reduce the computational cost from $O(N^2)$ to $O(N\log N)$, where $N$ is the number of particles. The convergence properties of the proposed algorithms are analyzed. Finally, we present several examples for both exact and inexact matchings, and numerically analyze the iterative process to illustrate the efficiency and the robustness of the proposed algorithms.

preprint2015arXiv

Solitary waves and $N$-particle algorithms for a class of Euler-Poincaré equations

We study a class of partial differential equations (PDEs) in the family of the so-called Euler-Poincaré differential systems, with the aim of developing a foundation for numerical algorithms of their solutions. This requires particular attention to the mathematical properties of this system when the associated class of elliptic operators possesses non-smooth kernels. By casting the system in its Lagrangian (or characteristics) form, we first formulate a particles system algorithm in free space with homogeneous Dirichlet boundary conditions for the evolving fields. We next examine the deformation of the system when non-homogeneous "constant stream" boundary conditions are assumed. We show how this simple change at the boundary deeply affects the nature of the evolution, from hyperbolic-like to dispersive with a non-trivial dispersion relation, and examine the potentially regularizing properties of singular kernels offered by this deformation. From the particle algorithm viewpoint, kernel singularities affect the existence and uniqueness of solutions to the corresponding ordinary differential equations systems. We illustrate this with the case when the operator kernel assumes a conical shape over the spatial variables, and examine in detail two-particle dynamics under the resulting lack of Lipschitz-continuity. Curiously, we find that for the conically-shaped kernels the motion of the related two-dimensional waves can become completely integrable under appropriate initial data. This reduction projects the two-dimensional system to the one-dimensional completely integrable Shallow-Water equation [Camassa, R. and Holm, D. D., Phys. Rev. Lett., 71, 1961-1964, 1993], while retaining the full dependence on two spatial dimensions for the single channel solutions.

preprint2014arXiv

A connection between the shallow-water equations and the Euler-Poincaré equations

The Euler-Poincaré differential (EPDiff) equations and the shallow water (SW) equations share similar wave characteristics. Using the Hamiltonian structure of the SW equations with flat bottom topography, we establish a connection between the EPDiff equations and the SW equations in one and multi-dimensions. Additionally, we show that the EPDiff equations can be recast in a curl formulation.

preprint2012arXiv

An inertia 'paradox' for incompressible stratified Euler fluids

The interplay between incompressibility and stratification can lead to non-conservation of horizontal momentum in the dynamics of a stably stratified incompressible Euler fluid filling an infinite horizontal channel between rigid upper and lower plates. Lack of conservation occurs even though in this configuration only vertical external forces act on the system. This apparent paradox was seemingly first noticed by Benjamin (J. Fluid Mech., vol. 165, 1986, pp. 445-474) in his classification of the invariants by symmetry groups with the Hamiltonian structure of the Euler equations in two dimensional settings, but it appears to have been largely ignored since. By working directly with the motion equations, the paradox is shown here to be a consequence of the rigid lid constraint coupling through incompressibility with the infinite inertia of the far ends of the channel, assumed to be at rest in hydrostatic equilibrium. Accordingly, when inertia is removed by eliminating the stratification, or, remarkably, by using the Boussinesq approximation of uniform density for the inertia terms, horizontal momentum conservation is recovered. This interplay between constraints,action at a distance by incompressibility, and inertia is illustrated by layer-averaged exact results, two-layer long-wave models, and direct numerical simulations of the incompressible Euler equations with smooth stratification.