An official website of the United States government
We consider the Euler equations for the potential flow of an ideal incompressible fluid of infinite depth with a free surface in two-dimensional geometry. Both gravity and surface tension forces are taken into account. A time-dependent conformal mapping is used which maps the lower complex half-plane of the auxiliary complex variable $$w$$ into the fluid’s area, with the real line of $$w$$ mapped into the free fluid’s surface. We reformulate the exact Eulerian dynamics through a non-canonical non-local Hamiltonian structure for a pair of the Hamiltonian variables. These two variables are the imaginary part of the conformal map and the fluid’s velocity potential, both evaluated at the fluid’s free surface. The corresponding Poisson bracket is non-degenerate, i.e. it does not have any Casimir invariant. Any two functionals of the conformal mapping commute with respect to the Poisson bracket. The new Hamiltonian structure is a generalization of the canonical Hamiltonian structure of Zakharov ( J. Appl. Mech. Tech. Phys. , vol. 9(2), 1968, pp. 190–194) which is valid only for solutions for which the natural surface parametrization is single-valued, i.e. each value of the horizontal coordinate corresponds only to a single point on the free surface. In contrast, the new non-canonical Hamiltonian equations are valid for arbitrary nonlinear solutions (including multiple-valued natural surface parametrization) and are equivalent to the Euler equations. We also consider a generalized hydrodynamics with the additional physical terms in the Hamiltonian beyond the Euler equations. In that case we identify powerful reductions that allow one to find general classes of particular solutions.
more »
« less
Poisson structure of the three-dimensional Euler equations in Fourier space
We derive a simple Poisson structure in the space of Fourier modes for the vorticity formulation of the Euler equations on a three-dimensional periodic domain. This allows us to analyse the structure of the Euler equations using a Hamiltonian framework. The Poisson structure is valid on the divergence free subspace only, and we show that using a projection operator it can be extended to be valid in the full space. We then restrict the simple Poisson structure to the divergence-free subspace on which the dynamics of the Euler equations take place, reducing the size of the system of ODEs by a third. The projected and the restricted Poisson structures are shown to have the helicity as a Casimir invariant. We conclude by showing that periodic shear flows in three dimensions are equilibria that correspond to singular points of the projected Poisson structure, and hence that the usual approach to study their nonlinear stability through the Energy-Casimir method fails.
more »
« less
- Award ID(s):
- 1812481
- PAR ID:
- 10106808
- Date Published:
- Journal Name:
- Journal of Physics A: Mathematical and Theoretical
- ISSN:
- 1751-8113
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We propose Coadjoint Orbit FLIP (CO-FLIP), a high order accurate, structure preserving fluid simulation method in the hybrid Eulerian-Lagrangian framework. We start with a Hamiltonian formulation of the incompressible Euler Equations, and then, using a local, explicit, and high order divergence free interpolation, construct a modified Hamiltonian system that governs our discrete Euler flow. The resulting discretization, when paired with a geometric time integration scheme, is energy and circulation preserving (formally the flow evolves on a coadjoint orbit) and is similar to the Fluid Implicit Particle (FLIP) method. CO-FLIP enjoys multiple additional properties including that the pressure projection is exact in the weak sense, and the particle-to-grid transfer is an exact inverse of the grid-to-particle interpolation. The method is demonstrated numerically with outstanding stability, energy, and Casimir preservation. We show that the method produces benchmarks and turbulent visual effects even at low grid resolutions.more » « less
-
null (Ed.)We adapt the multiple scale method introduced over 40 years ago for the homogenization of periodic structures [1], to the quasiperiodic (cut-and-projection) setting. We make use of partial differential operators (gradient, divergence and curl) acting on periodic functions of m variables in a higher-dimensional space that are projected onto operators acting on quasiperiodic functions in the n-dimensional physical space (m>n). We replace heterogeneous quasiperiodic structures, coined irrational metamaterials in [2], by homogeneous media with anisotropic permittivity and permeability tensors, obtained from the solution of annex problems of electrostatic type in a periodic cell in higher dimensional space. This approach is valid when the wavelength is much larger than the period of the higher dimensional elementary cell.more » « less
-
This paper is concerned with the critical threshold phenomenon for one-dimensional damped, pressureless Euler–Poisson equations with electric force induced by a constant background, originally studied in [S. Engelberg and H. Liu and E. Tadmor, Indiana Univ. Math. J. 50 (2001) 109–157]. A simple transformation is used to linearize the characteristic system of equations, which allows us to study the geometrical structure of critical threshold curves for three damping cases: overdamped, underdamped and borderline damped through phase plane analysis. We also derive the explicit form of these critical curves. These sharp results state that if the initial data is within the threshold region, the solution will remain smooth for all time, otherwise it will have a finite time breakdown. Finally, we apply these general results to identify critical thresholds for a non-local system subjected to initial data on the whole line.more » « less
-
This paper considered divergence-free basis methods to solve the viscous Stokes equations. A discrete divergence-free subspace was constructed to reduce the saddle point problem of the Stokes problem to a smaller-sized symmetric and positive definite system solely depending on the velocity components. Then, the system could decouple the unknowns in velocity and pressure and solve them independently. However, such a scheme may not ensure an accurate numerical solution to the velocity. In order to obtain satisfactory accuracy, we used a velocity reconstruction technique to enhance the divergence-free scheme to achieve the desired pressure and viscosity robustness. Numerical results were presented to demonstrate the robustness and accuracy of this discrete divergence-free method.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1088/1751-8121/ab3363
