🔢 mathematics

On an Initial Value Problem Describing the Small Oscillations of a Floating Cylinder

This paper establishes the well-posedness of an initial value problem for a coupled PDE-ODE system modeling the small oscillations of an infinite circular cylinder interacting with water waves, utilizing an abstract evolution equation framework and a novel analysis of a partial Dirichlet-to-Neumann map on an unbounded domain with a non-smooth boundary.

Original authors: Vicente Ocqueteau (IMB), Marius Tucsnak (IMB)

Published 2026-08-03

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Vicente Ocqueteau (IMB), Marius Tucsnak (IMB)

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Dance of the Floating Giant

Imagine the ocean as a giant, invisible trampoline. When you drop a stone, ripples spread out; when a ship bobs up and down, it creates waves that push back against the ship. This constant tango between a floating object and the water around it is the heart of fluid dynamics, a branch of physics that tries to predict how liquids move. For engineers designing floating wind turbines or massive oil rigs, understanding this dance is crucial. If they get the math wrong, the structure could shake itself apart or fail to generate power.

The tricky part is that water doesn't just sit still; it flows, and when a solid object (like a cylinder) sits in it, the water has to squeeze around the object's corners. In the real world, these corners create "kinks" in the flow that are mathematically messy. Most previous studies looked at smooth, perfect shapes or simplified, shallow water. But what happens when you have a giant, infinite cylinder bobbing in deep water, with sharp edges where the water meets the metal? That's the specific, knotty problem this paper tackles. It asks: Can we write a perfect set of rules (equations) that predict exactly how this cylinder and the water will move together over time, without the math breaking down?

The Paper's Big Discovery

This paper, written by Vicente Ocqueteau and Marius Tucsnak, is like a master locksmith finally finding the right key for a very stubborn, complex lock. The authors study a system where an infinite circular cylinder floats in deep water, with exactly half of it submerged. They are interested in "small oscillations," which is a fancy way of saying the cylinder is bobbing up and down gently, not crashing through the waves.

The main challenge they faced was the "corner" where the cylinder meets the water surface. In mathematics, corners are notorious for making equations behave badly, especially when the water goes on forever (an unbounded domain). Previous attempts to solve similar problems often hit a wall here. The authors' breakthrough is proving that, despite these corners and the infinite depth, the system is "well-posed."

In plain English, "well-posed" means three things:

Existence: A solution actually exists. The cylinder and water will move in a predictable way; the math doesn't just say "impossible."
Uniqueness: There is only one correct way for them to move. If you start with the same push, you get the same result every time.
Stability: Small changes in the starting push lead to small changes in the result, not chaotic explosions.

To prove this, the authors didn't just guess; they built a rigorous mathematical bridge. They treated the water and the cylinder as a single, coupled machine. They used a clever trick involving "potentials" (imaginary fields that describe the water's speed) to turn the messy water equations into a cleaner format. A key part of their work involved a tool called a "Dirichlet-to-Neumann map." Think of this as a translator that takes information about the water's surface height and instantly tells you how hard the water is pushing down on the bottom, without having to calculate every single drop of water in between.

The authors proved that this translator works perfectly even with the cylinder's sharp corners and the infinite ocean. They showed that this map acts like a reliable, positive force that keeps the system stable. By organizing all these pieces into a "group of operators" (a mathematical way of describing how a system evolves over time), they demonstrated that the initial value problem—the question of "if we start here, where do we go next?"—has a unique, stable solution.

They didn't just say "it works"; they proved it using standard mathematical spaces called Sobolev spaces, which are like specific containers for functions that behave nicely. They showed that if you start with a gentle bob and a calm water surface, the system will continue to evolve smoothly forever. They also showed that if you add a gentle, rhythmic push (like a wave or a motor), the system responds in a predictable, continuous way.

In short, this paper doesn't just simulate a floating cylinder; it provides the mathematical guarantee that the laws of physics governing this specific, tricky scenario are solid and reliable. It clears the path for engineers to trust the math when designing floating structures that interact with deep, complex waves, ensuring that the "dance" between the steel and the sea is one we can predict with confidence.

Technical Summary: On an Initial Value Problem Describing the Small Oscillations of a Floating Cylinder

Problem Formulation
This paper investigates the well-posedness of an initial value problem (IVP) modeling the coupled motion of a floating infinite circular cylinder and the free surface of an ideal fluid of infinite depth. The study operates within the linear potential flow framework, assuming irrotational flow, constant density, and small amplitude oscillations. The geometric configuration involves a cylinder with a unit radius, where the equilibrium position corresponds to half-immersion. The fluid domain $\Omega$ is the exterior of the unit disk in the upper half-plane ( $y>0$ ), and the free surface $E$ consists of the real axis excluding the interval $[-1, 1]$ .

The system couples a partial differential equation (PDE) governing the fluid velocity potential $\phi$ with an ordinary differential equation (ODE) describing the vertical displacement $h(t)$ of the cylinder's center. The governing equations include:

The Laplace equation for $\phi$ in $\Omega$ .
A linearized free-surface condition derived from Bernoulli's law on $E$ .
A kinematic boundary condition on the cylinder surface relating the fluid radial velocity to the cylinder's vertical velocity $\dot{h}(t)$ .
Newton's second law for the cylinder, balancing hydrostatic restoring forces, hydrodynamic pressure integrals, and external forcing.

The authors reformulate this coupled PDE-ODE system into an abstract evolution equation in a suitable Hilbert space. The primary unknowns are reduced to the trace of the potential on the free surface, $v(t, x) = \phi(t, x, 0)$ , and the cylinder displacement $h(t)$ .

Methodology
The core of the methodology relies on the rigorous construction and analysis of a Dirichlet-to-Neumann (DtN) map, denoted $\Lambda_\Omega$ , on the unbounded domain $\Omega$ with a non-smooth boundary (specifically, corners at the intersection of the cylinder and the free surface).

Decomposition and Explicit Potentials: The velocity potential is decomposed into two parts: one generated by the free surface oscillations and another by the cylinder's motion. The authors derive explicit formulas for these potentials. The cylinder-induced potential $\phi_1$ is given explicitly, while the free-surface potential is constructed via a Dirichlet operator $D_\Omega$ .
Construction of $D_\Omega$ : To handle the mixed boundary conditions (Dirichlet on the free surface, Neumann on the cylinder), the authors employ a symmetry argument. They extend the problem to the full upper half-plane $H$ by reflecting the data across the unit circle. This allows them to express $D_\Omega$ using the classical Poisson kernel for the half-plane, effectively solving the mixed problem via the standard Dirichlet problem on $H$ .
Analysis of $\Lambda_\Omega$ : The DtN map $\Lambda_\Omega$ is defined as the normal derivative of the solution on the free surface. Rather than defining it purely by duality, the authors relate $\Lambda_\Omega$ to the classical DtN map $\Lambda_H$ of the upper half-plane. They derive an explicit representation of $\Lambda_\Omega$ involving the Hilbert transform and specific correction terms arising from the geometry.
Operator Semigroup Approach: The coupled system is cast as a first-order evolution equation $\dot{z}(t) = A_1 z(t) + F(t)$ in a product Hilbert space $X$ . The operator $A_1$ is constructed from a self-adjoint, positive operator $A_0$ (involving $\Lambda_\Omega$ ) and a bounded perturbation $P$ . The well-posedness is established by proving that $A_0$ generates a $C_0$ -group (via the Lumer-Phillips theorem or similar semigroup results) and that the perturbation preserves this property.

Key Contributions and Results
The paper establishes the following main results:

Well-posedness in Sobolev Spaces: The authors prove that the initial value problem is well-posed in standard Sobolev spaces. Specifically, for initial data $v_0 \in W^{1/2, 2}(E)$ , $v_1 \in L^2(E)$ , and $h_0, h_1 \in \mathbb{R}$ , there exists a unique solution $(v, h)$ such that $v \in C([0, \infty), W^{1/2, 2}(E)) \cap C^1([0, \infty), L^2(E))$ and $h \in C^1([0, \infty))$ . Higher regularity is obtained for smoother initial data and forcing terms.
Properties of the DtN Map: A significant technical contribution is the proof that the operator $\Lambda_\Omega$ is an unbounded, self-adjoint, and positive operator on $L^2(E)$ with domain $W^{1, 2}(E)$ . This property is crucial for defining the energy space and constructing the semigroup.
Explicit Formulas: The paper provides explicit integral representations for the Dirichlet operator $D_\Omega$ and the DtN map $\Lambda_\Omega$ , avoiding purely abstract definitions and allowing for concrete analysis of the boundary terms.

Significance and Claims
The authors position this work as addressing a gap in the literature regarding the time-domain analysis of fluid-structure interactions within the linear potential flow approach. While previous works have largely focused on time-harmonic (frequency domain) motions or reduced models (shallow water), this paper explicitly treats the initial value problem for the full linear system.

The significance lies in demonstrating that, despite the geometric complexities introduced by the "corners" where the cylinder meets the free surface and the unbounded nature of the domain, the system remains well-posed in standard functional analytic settings. The authors note that the specific geometric configuration (infinite cylinder, infinite depth) allows for the derivation of explicit formulas, which simplifies the technical difficulties compared to more general configurations. The proof of well-posedness relies on the fact that the associated DtN map defines a positive operator, enabling the use of semigroup theory to describe the system's dynamics as a $C_0$ group.

The paper does not claim to solve the general nonlinear problem or extend to finite depth or deformable bodies, nor does it propose new numerical schemes or experimental validations. Its contribution is strictly theoretical, providing a rigorous mathematical foundation for the linearized time-domain model of a floating cylinder.

The Dance of the Floating Giant

The Paper's Big Discovery

More like this