Global well‐posedness and finite‐dimensional global attractor for a 3‐D planetary geostrophic viscous model

In this paper we consider a three‐dimensional planetary geostrophic viscous model of the gyre‐scale mid‐latitude ocean. We show the global existence and uniqueness of the weak and strong solutions to this model. Moreover, we establish the existence of a finite‐dimensional global attractor to this dissipative evolution system. © 2002 Wiley Periodicals, Inc.     

Surface waves on steady perfect‐fluid flows with vorticity

This is a theory of two‐dimensional steady periodic surface waves on flows under gravity in which the given data are three quantities that are independent of time in the corresponding evolution problem: the volume of fluid per period, the circulation per period on the free stream line, and the rearrangement class (equivalently, the distribution function) of the vorticity field. A minimizer of the total energy per period among flows satisfying these three constraints is shown to be a weak solution of the surface wave problem for which the vorticity is a decreasing function of the stream function. This decreasing function can be thought of as an infinite‐dimensional Lagrange multiplier corresponding to the vorticity rearrangement class being specified in the minimization problem. (Note that functional dependence of vorticity on the stream function was not specified a priori but is part of the solution to the problem and ensures the flow is steady.) To illustrate the idea with a minimum of technical difficulties, the existence of nontrivial waves on the surface of a fluid flowing with a prescribed distribution of vorticity and confined beneath an elastic sheet is proved. The theory applies equally to irrotational flows and to flows with locally square‐integrable vorticity. © 2011 Wiley Periodicals, Inc.     

Nonlinear Stability of Baroclinic Fronts in a Channel with Variable Topography

The governing equations describing baroclinic bottom-trapped fronts in a channel with variable bottom topography are shown to be a noncanonical Hamiltonian system. The Hamiltonian formalism is exploited to derive a variational principle for arbitrary steady solutions based on an appropriately constrained energy functional. The variational principle is exploited to obtain formal and nonlinear stability conditions. In the infinitesimal amplitude limit, these stability conditions reduce to previously obtained normal mode results for the transverse gradient of the mean frontal potential vorticity.     

Rotational waves generated by current‐topography interaction

We study nonlinear free‐surface rotational waves generated through the interaction of a vertically sheared current with a topography. Equivalently, the waves may be generated by a pressure distribution along the free surface. A forced Korteweg–de Vries equation (fKdV) is deduced incorporating these features. The weakly nonlinear, weakly dispersive reduced model is valid for small amplitude topographies. To study the effect of gradually increasing the topography amplitude, the free surface Euler equations are formulated in the presence of a variable depth and a sheared current of constant vorticity. Under constant vorticity, the harmonic velocity component is formulated in a simplified canonical domain, through the use of a conformal mapping which flattens both the free surface as well as the bottom topography. Critical, supercritical, and subcritical Froude number regimes are considered, while the bottom amplitude is gradually increased in both the irrotational and rotational wave regimes. Solutions to the fKdV model are compared to those from the Euler equations. We show that for rotational waves the critical Froude number is shifted away from 1. New stationary solutions are found and their stability tested numerically.     

Asymptotic and numerical solutions of three‐dimensional boundary‐layer flow past a moving wedge

We consider a laminar boundary‐layer flow of a viscous and incompressible fluid past a moving wedge in which the wedge is moving either in the direction of the mainstream flow or opposite to it. The mainstream flows outside the boundary layer are approximated by a power of the distance from the leading boundary layer. The variable pressure gradient is imposed on the boundary layer so that the system admits similarity solutions. The model is described using 3‐dimensional boundary‐layer equations that contains 2 physical parameters: pressure gradient (β) and shear‐to‐strain‐rate ratio parameter (α). Two methods are used: a linear asymptotic analysis in the neighborhood of the edge of the boundary layer and the Keller‐box numerical method for the full nonlinear system. The results show that the flow field is divided into near‐field region (mainly dominated by viscous forces) and far‐field region (mainstream flows); the velocity profiles form through an interaction between 2 regions. Also, all simulations show that the subsequent dynamics involving overshoot and undershoot of the solutions for varying parameter characterizing 3‐dimensional flows. The pressure gradient (favorable) has a tendency of decreasing the boundary‐layer thickness in which the velocity profiles are benign. The wall shear stresses increase unboundedly for increasing α when the wedge is moving in the x‐direction, while the case is different when it is moving in the y‐direction. Further, both analysis show that 3‐dimensional boundary‐layer solutions exist in the range −1<α<∞. These are some interesting results linked to an important class of boundary‐layer flows.     

Mixed formulation of the two‐layer quasi‐geostrophic equations of the ocean

In this article, we propose a mixed variational formulation for the streamfunction vorticity potential form for the two‐layer quasi‐geostrophic model of the ocean. We prove the existence and uniqueness of solutions of the mixed variational problem. © 1999 John Wiley & Sons, Inc. Numer Methods Partial Differential Eq 15: 489–502, 1999     

Rate‐type Model for electro‐rheological Material Behaviour consistent with Extended Thermodynamics: Application to a steady viscometric flow

Based on extended thermodynamics, a class of rate‐type constitutive equations for electrorheological materials is presented. It is shown unlike many recent works using standard methods of continuum mechanics, that this approach, which consists in postulating expressions which relate the time derivative of the Cauchy stress to the motion, takes well account of the viscoelastic‐plastic effects as well as the coupling between the motion and the electric field. By considering special cases it turns out that the derived nonlinear algebraic model embeds a variety of existing differential models in the literature. The accuracy of the model is illustrated by considering a steady shear flow in which the electric field is applied perpendicularly to the flow direction. Analytical solutions with various parameter variations are obtained and discussed. The model predicts well the behaviour associated with material instabilities observed in viscometric flows for which the existence of the yield stress is required as confirmed by many experiments.     

On the structure and growth rate of unstable modes to the Rossby‐Haurwitz wave

The normal mode instability study of a steady Rossby‐Haurwitz wave is considered both theoretically and numerically. This wave is exact solution of the nonlinear barotropic vorticity equation describing the dynamics of an ideal fluid on a rotating sphere, as well as the large‐scale barotropic dynamics of the atmosphere. In this connection, the stability of the Rossby‐Haurwitz wave is of considerable mathematical and meteorological interest. The structure of the spectrum of the linearized operator in case of an ideal fluid is studied. A conservation law for perturbations to the Rossby‐Haurwitz wave is obtained and used to get a necessary condition for its exponential instability. The maximum growth rate of unstable modes is estimated. The orthogonality of the amplitude of a non‐neutral or non‐stationary mode to the Rossby‐Haurwitz wave is shown in two different inner products. The analytical results obtained are used to test and discuss the accuracy of a numerical spectral method used for the normal mode stability study of arbitrary flow on a sphere. The comparison of the numerical and theoretical results shows that the numerical instability study method works well in case of such smooth solutions as the zonal flows and Rossby‐Haurwitz waves. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Eulerian Simulation of Gas‐Solid Particle Flow by BDIM

The numerical scheme based on the boundary domain integral method (BDIM) for the numerical simulation of twophase two‐component flows is presented. A program is being developed to model the hydrodynamics of fluidized bed systems by using the Eulerian approach in terms of velocity‐vorticity variables formulation. With the vorticity vector both phases motion computation scheme is partitioned into its kinematic and kinetic aspect. Influence of the drag coefficient on the two‐phase two‐component flow field is studied on the two‐phase gas‐solid particles vertical channel flow.     

The Bernoulli Integral for a Certain Class of Non‐Stationary Viscous Vortical Flows of Incompressible Fluid

It has been shown in our paper [1] that there is a wide class of 3D motions of incompressible viscous fluid which can be described by one scalar function dabbed the quasi‐potential. This class of fluid flows is characterized by three‐component velocity field having two‐component vorticity field; both these fields can depend of all three spatial variables and time, in general. Governing equations for the quasi‐potential have been derived and simple illustrative example of 3D flow has been presented. Here, we derive the Bernoulli integral for that class of flows and compare it against the known Bernoulli integrals for the potential flows or 2D stationary vortical flows of inviscid fluid. We show that the Bernoulli integral for this class of fluid motion possesses unusual features: it is valid for the vortical nonstationary motions of a viscous incompressible fluid. We present a new very nontrivial analytical example of 3D flow with two‐component vorticity which hardly can be obtained by any of known methods. In the last section, we suggest a generalization of the developed concept which allows one to describe a certain class of 3D flows with the 3D vorticity.     

On mean flow generation due to oblique reflection of internal waves at a slope

Small‐amplitude expansions are utilized to discuss the mean flow induced by the reflection of a weakly nonlinear internal gravity wave beam at a uniform rigid slope, in the case where the beam planes of constant phase meet the slope at an arbitrary direction, not necessarily parallel to the isobaths, and the flow cannot be taken as two dimensional. Along the vertical, the Eulerian mean flow, due to such an oblique reflection, is equal and opposite to the Stokes drift so the Lagrangian mean flow vanishes, similar to a two‐dimensional reflection. The horizontal Eulerian mean flow, however, is controlled by the mean potential vorticity (PV) and the corresponding Lagrangian mean flow is generally nonzero, in contrast to two‐dimensional flow where PV identically vanishes. For an oblique reflection, furthermore, viscous dissipation can trigger generation of horizontal mean flow via irreversible production of mean PV, a phenomenon akin to streaming.     

Rossby Solitary Waves in the Presence of a Critical Layer

This study considers the evolution of weakly nonlinear long Rossby waves in a horizontally sheared zonal current. We consider a stable flow so that the nonlinear time scale is long. These assumptions enable the flow to organize itself into a large‐scale coherent structure in the régime where a competition sets in between weak nonlinearity and weak dispersion. This balance is often described by a Korteweg‐de‐Vries equation. The traditional assumption of a weak amplitude breaks down when the wave speed equals the mean flow velocity at a certain latitude, due to the appearance of a singularity in the leading‐order equation, which strongly modifies the flow in a critical layer. Here, nonlinear effects are invoked to resolve this singularity, because the relevant geophysical flows have high Reynolds numbers. Viscosity is introduced in order to render the nonlinear‐critical‐layer solution unique, but the inviscid limit is eventually taken. By the method of matched asymptotic expansions, this inner flow is matched at the edges of the critical layer with the outer flow. We will show that the critical‐layer–induced flow leads to a strong rearrangement of the related streamlines and consequently of the potential‐vorticity contours, particularly in the neighborhood of the separatrices between the open and closed streamlines. The symmetry of the critical layer vis‐à‐vis the critical level is also broken. This theory is relevant for the phenomenon of Rossby wave breaking and eventual saturation into a nonlinear wave. Spatially localized solutions are described by a Korteweg‐de‐Vries equation, modified by new nonlinear terms; depending on the critical‐layer shape, this leads to depression or elevation waves. The additional terms are made necessary at a certain order of the asymptotic expansion while matching the inner flow on the dividing streamlines. The new evolution equation supports a family of solitary waves. In this paper we describe in detail the case of a depression wave, and postpone for further discussion the more complex case of an elevation wave.     

Stokes waves with constant vorticity: I. Numerical computation

Periodic traveling waves are numerically computed in a constant vorticity flow subject to the force of gravity. The Stokes wave problem is formulated via a conformal mapping as a nonlinear pseudodifferential equation, involving a periodic Hilbert transform for a strip, and solved by the Newton‐GMRES method. For strong positive vorticity, in the finite or infinite depth, overhanging profiles are found as the amplitude increases and tend to a touching wave, whose surface contacts itself at the trough line, enclosing an air bubble; numerical solutions become unphysical as the amplitude increases further and make a gap in the wave speed versus amplitude plane; another touching wave takes over and physical solutions follow along the fold in the wave speed versus amplitude plane until they ultimately tend to an extreme wave, which exhibits a corner at the crest. Touching waves connected to zero amplitude are found to approach the limiting Crapper wave as the strength of positive vorticity increases unboundedly, while touching waves connected to the extreme waves approach the rigid body rotation of a fluid disk.     

Gamma‐convergence of gradient flows with applications to Ginzburg‐Landau

We present a method to prove convergence of gradient flows of families of energies that Γ‐converge to a limiting energy. It provides lower‐bound criteria to obtain the convergence that correspond to a sort of C¹‐order Γ‐convergence of functionals. We then apply this method to establish the limiting dynamical law of a finite number of vortices for the heat flow of the Ginzburg‐Landau energy in dimension 2, retrieving in a different way the existing results for the case without magnetic field and obtaining new results for the case with magnetic field. © 2004 Wiley Periodicals, Inc.     

Conservation Laws,Hamiltonian Structure,Modulational Instability Properties and Solitary Wave Solutions for a Higher‐Order Model Describing Nonlinear Internal Waves

Recent theoretical advances in connecting the wave‐induced mean flow with the conserved pseudomomentum per unit mass has permitted the first rational derivation of a model that describes the weakly nonlinear propagation of internal gravity plane waves in a continuously stratified fluid. Depending on the particular parameter regime examined the new model corresponds to an extended bright or dark derivative nonlinear Schrödinger equation or an extended complex‐valued modified Korteweg‐de Vries or Sasa–Satsuma equation. Mass, momentum, and energy conservation laws are derived. A noncanonical infinite‐dimensional Hamiltonian formulation of the model is introduced. The modulational stability characteristics associated with the Stokes wave solution of the model are described. The bright and dark solitary wave solutions of the model are obtained.     

Global solutions of the Navier-Stokes equations in a uniformly rotating space

We consider the Cauchy problem for the Navier-Stokes system of equations in a three-dimensional space rotating uniformly about the vertical axis with the periodicity condition with respect to the spatial variables. Studying this problem is based on expanding given and sought vector functions in Fourier series in terms of the eigenfunctions of the curl and Stokes operators. Using the Galerkin method, we reduce the problem to the Cauchy problem for the system of ordinary differential equations, which has a simple explicit form in the basis under consideration. Its linear part is diagonal, which allows writing explicit solutions of the linear Stokes-Sobolev system, to which fluid flows with a nonzero vorticity correspond. Based on the study of the nonlinear interaction of vortical flows, we find an approach that we can use to obtain families of explicit global solutions of the nonlinear problem.     

An optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows

This article is concerned about an optimization‐based domain decomposition method for numerical simulation of the incompressible Navier‐Stokes flows. Using the method, an classical domain decomposition problem is transformed into a constrained minimization problem for which the objective functional is chosen to measure the jump in the dependent variables across the common interfaces between subdomains. The Lagrange multiplier rule is used to transform the constrained optimization problem into an unconstrained one and that rule is applied to derive an optimality system from which optimal solutions may be obtained. The optimality system is also derived using “sensitivity” derivatives instead of the Lagrange multiplier rule. We consider a gradient‐type approach to the solution of domain decomposition problem. The results of some numerical experiments are presented to demonstrate the feasibility and applicability of the algorithm developed in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Averaging over fast gravity waves for geophysical flows with arbitary

Here a mathematically rigorous framework is developed for deriving new reduced simplified dynamical equations for geophysical flows with arbitrary potential vorticity interacting with fast gravity waves. The examples include the rotating Boussinesq and rotating shallow water equations in the quasigeostrophic limit with vanishing Rossby number. For the spatial periodic case the theory implies that the quasi—geostrophic equations are valid limiting equations in the weak topology for arbitrary initial data. Furthermore, simplified reduced equations are developed for the fashion in which the vortical waves influence the gravity waves through averaging over specific gravity wave/vortical resonances.     

Numerical and asymptotic study of non‐axisymmetric magnetohydrodynamic boundary layer stagnation‐point flows

Both numerical and asymptotic analyses are performed to study the similarity solutions of three‐dimensional boundary‐layer viscous stagnation point flow in the presence of a uniform magnetic field. The three‐dimensional boundary‐layer is analyzed in a non‐axisymmetric stagnation point flow, in which the flow is developed because of influence of both applied magnetic field and external mainstream flow. Two approaches for the governing equations are employed: the Keller‐box numerical simulations solving full nonlinear coupled system and a corresponding linearized system that is obtained under a far‐field behavior and in the limit of large shear‐to‐strain‐rate parameter (λ). From these two approaches, the flow phenomena reveals a rich structure of new family of solutions for various values of the magnetic number and λ. The various results for the wall stresses and the displacement thicknesses are presented along with some velocity profiles in both directions. The analysis discovered that the flow separation occurs in the secondary flow direction in the absence of magnetic field, and the flow separation disappears when the applied magnetic field is increased. The flow field is divided into a near‐field (due to viscous forces) and far‐field (due to mainstream flows), and the velocity profiles form because of an interaction between two regions. The magnetic field plays an important role in reducing the thickness of the boundary‐layer. A physical explanation for all observed phenomena is discussed. Copyright © 2017 John Wiley & Sons, Ltd.     

Analytic approximate solution of three‐dimensional Navier‐Stokes equations of flow between two stretchable disks

The similarity transform for the steady three‐dimensional Navier‐Stokes equations of flow between two stretchable disks gives a system of nonlinear ordinary differential equations which is analytically solved by applying a newly developed method, namely, the homotopy analysis method. The analytic solutions of the system of nonlinear ordinary differential equations are constructed in the series form. The convergence of the obtained series solutions is analyzed. The validity of our solutions is verified by the numerical results. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010