No Arabic abstract
The simulation of large open water surface is challenging for a uniform volumetric discretization of the Navier-Stokes equation. The water splashes near moving objects, which height field methods for water waves cannot capture, necessitates high resolution simulation such as the Fluid-Implicit-Particle (FLIP) method. On the other hand, FLIP is not efficient for the long-lasting water waves that propagates to long distances, which requires sufficient depth for correct dispersion relationship. This paper presents a new method to tackle this dilemma through an efficient hybridization of volumetric and surface-based advection-projection discretizations. We design a hybrid time-stepping algorithm that combines a FLIP domain and an adaptively remeshed Boundary Element Method (BEM) domain for the incompressible Euler equations. The resulting framework captures the detailed water splashes near moving objects with FLIP, and produces convincing water waves with correct dispersion relationship at modest additional cost.
We calculate the rate of ocean waves energy dissipation due to whitecapping by numerical simulation of deterministic phase resolving model for dynamics of ocean surface. Two independent numerical experiments are performed. First, we solve the $3D$ Hamiltonian equation that includes three- and four-wave interactions. This model is valid for moderate values of surface steepness only, $mu < 0.09$. Then we solve the exact Euler equation for non-stationary potential flow of an ideal fluid with a free surface in $2D$ geometry. We use the conformal mapping of domain filled with fluid onto the lower half-plane. This model is applicable for arbitrary high levels of steepness. The results of both experiments are close. The whitecapping is the threshold process that takes place if the average steepness $mu > mu_{cr} simeq 0.055$. The rate of energy dissipation grows dramatically with increasing of steepness. Comparison of our results with dissipation functions used in the operational models of wave forecasting shows that these models overestimate the rate of wave dissipation by order of magnitude for typical values of steepness.
A formulation is developed to assimilate ocean-wave data into the Numerical Flow Analysis (NFA) code. NFA is a Cartesian-based implicit Large-Eddy Simulation (LES) code with Volume of Fluid (VOF) interface capturing. The sequential assimilation of data into NFA permits detailed analysis of ocean-wave physics with higher bandwidths than is possible using either other formulations, such as High-Order Spectral (HOS) methods, or field measurements. A framework is provided for assimilating the wavy and vortical portions of the flow. Nudging is used to assimilate wave data at low wavenumbers, and the wave data at high wavenumbers form naturally through nonlinear interactions, wave breaking, and wind forcing. Similarly, the vertical profiles of the mean vortical flow in the wind and the wind drift are nudged, and the turbulent fluctuations are allowed to form naturally. As a demonstration, the results of a HOS of a JONSWAP wave spectrum are assimilated to study short-crested seas in equilibrium with the wind. Log profiles are assimilated for the mean wind and the mean wind drift. The results of the data assimilations are (1) Windrows form under the action of breaking waves and the formation of swirling jets; (2) The crosswind and cross drift meander; (3) Swirling jets are organized into Langmuir cells in the upper oceanic boundary layer; (4) Swirling jets are organized into wind streaks in the lower atmospheric boundary layer; (5) The length and time scales of the Langmuir cells and the wind streaks increase away from the free surface; (6) Wave growth is very dynamic especially for breaking waves; (7) The effects of the turbulent fluctuations in the upper ocean on wave growth need to be considered together with the turbulent fluctuations in the lower atmosphere; and (8) Extreme events are most likely when waves are not in equilibrium.
The influence of forward speed on stochastic free-surface crossing, in a Gaussian wave field, is investigated. The case of a material point moving with a constant forward speed is considered; the wave field is assumed stationary in time, and homogeneous in space. The focus is on up-crossing events, which are defined as the material point crossing the free surface, into the water domain. The effect of the Doppler shift (induced by the forward speed) on the up-crossing frequency, and the related conditional joint distribution of wave kinematic variables is analytically investigated. Some general trends are illustrated through different examples, where three kinds of wave direction distribution are considered: unidirectional, short-crested anisotropic, and isotropic. The way the developed approach may be used in the context of slamming on marine structures is briefly discussed.
We investigate a new two-dimensional compressible Navier-Stokes hydrodynamic model design to explain and study large scale ice swirls formation at the surface of the ocean. The linearized model generates a basis of Bessel solutions from where various types of spiral patterns can be generated and their evolution and stability in time analyzed. By restricting the nonlinear system of equations to its quadratic terms we obtain swirl solutions emphasizing logarithmic spiral geometry. The resulting solutions are analyzed and validated using three mathematical approaches: one predicting the formation of patterns as Townes solitary modes, another approach mapping the nonlinear system into a sine-Gordon equation, and a third approach uses a series expansion. Pure radial, azimuthal and spiral modes are obtained from the fully nonlinear equations. Combinations of multiple-spiral solutions are also obtained, matching the experimental observations. The nonlinear stability of the spiral patterns is analyzed by Arnolds convexity method, and the Hamiltonian of the solutions is plotted versus some order parameters showing the existence of geometric phase transitions.
The multifractal theory of turbulence is used to investigate the energy cascade in the Northwestern Atlantic ocean. The statistics of singularity exponents of velocity gradients computed from in situ measurements are used to show that the anomalous scaling of the velocity structure functions at depths between 50 ad 500 m has a linear dependence on the exponent characterizing the strongest velocity gradient, with a slope that decreases with depth. Since the distribution of exponents is asymmetric about the mode at all depths, we use an infinitely divisible asymmetric model of the energy cascade, the log-Poisson model, to derive the functional dependence of the anomalous scaling with dissipation. Using this model we can interpret the vertical change of the linear slope as a change in the energy cascade.