Mechanical Balance Laws for Boussinesq Models of Surface Water Waves

Depth-integrated long-wave models, such as the shallow-water and Boussinesq equations, are standard fare in the study of small amplitude surface waves in shallow water. While the shallow-water theory features conservation of mass, momentum and energy for smooth solutions, mechanical balance equations are not widely used in Boussinesq scaling, and it appears that the expressions for many of these quantities are not known. This work presents a systematic derivation of mass, momentum and energy densities and fluxes associated with a general family of Boussinesq systems. The derivation is based on a reconstruction of the velocity field and the pressure in the fluid column below the free surface, and the derivation of differential balance equations which are of the same asymptotic validity as the evolution equations. It is shown that all these mechanical quantities can be expressed in terms of the principal dependent variables of the Boussinesq system: the surface excursion ?? and the horizontal velocity w at a given level in the fluid.     

两层流体界面上的孤立波

本文讨论两水平固壁间两层不可压无粘流体界面上的孤立波,计及界面上的表面张力效应.首先建立了适用于这种模型的基本方程组,并在弱色散近似下应用约化摄动法,导得了一阶界面升高所满足的Korteweg-de Vries方程,指出了按该方程系数α和μ的符号的异同,KdV孤立波可能凸向上或凸向下.然后详细讨论了原有近似下非线性效应与色散效应不能平衡的两种临界情形.在采用了适当的近似之后,对第一种临界情形(α=0)得到了修正的KdV方程,并指出,在所考虑的情形中,当μ>0时孤立波不存在,当μ<0时,孤立波仍可能存在,其形式与KdV孤立波不同;对第二种临界情形(μ=0),导得了推广的KdV方程,这时存在振荡型孤立波.文中还对近临界情形作了讨论.本文结果与一些经典结果完全一致,并把它们作了拓广.     

Numerical stability of the KdV equation with a negative forcing

The waves at the free surface waves of an incompressible and inviscid fluid in a two dimensional domain with horizontal rigid flat bottom with a small obstruction are considered. A time dependent KdV equation with a negative forcing is derived and studied both theoretically and numerically. The existence of a negative solitary-wave-like solution of the equation near the Froude number is proved and the numerical stability of the solution is also studied. The numerical stability of the positive both symmetric and unsymmetric solitary-wave-like solutions are also studied. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rigorous estimates on mechanical balance laws in the Boussinesq–Peregrine equations

It is shown that the Boussinesq–Peregrine system, which describes long waves of small amplitude at the surface of an inviscid fluid with variable depth, admits a number of approximate conservation equations. Notably, this paper provides accurate estimations for the approximate conservation of the mechanical balance laws associated with mass, momentum, and energy. These precise estimates offer valuable insights into the behavior and dynamics of the system, shedding light on the conservation principles governing the wave motion.     

Solitary wave solutions of the high-order KdV models for bi-directional water waves

An approach, which allows us to construct specific closed-form solitary wave solutions for the KdV-like water-wave models obtained through the Boussinesq perturbation expansion for the two-dimensional water wave problem in the limit of long wavelength/small amplitude waves, is developed. The models are relevant to the case of the bi-directional waves with the amplitude of the left-moving wave of O(ϵ) (ϵ is the amplitude parameter) as compared with that of the right-moving wave. We show that, in such a case, the Boussinesq system can be decomposed into a system of coupled equations for the right- and left-moving waves in which, to any order of the expansion, one of the equations is dependent only on the (main) right-wave elevation and takes the form of the high-order KdV equation with arbitrary coefficients whereas the second equation includes both elevations. Then the explicit solitary wave solutions constructed via our approach may be treated as the exact solutions of the infinite-order perturbed KdV equations for the right-moving wave with the properly specified high-order coefficients. Such solutions include, in a sense, contributions of all orders of the asymptotic expansion and therefore may be considered to a certain degree as modelling the solutions of the original water wave problem under proper initial conditions. Those solitary waves, although stemming from the KdV solitary waves, possess features found neither in the KdV solitons nor in the solutions of the first order perturbed KdV equations.     

On linear and nonlinear Rossby waves in an ocean

This paper deals with recent developments of linear and nonlinear Rossby waves in an ocean. Included are also linear Poincaré, Rossby, and Kelvin waves in an ocean. The dispersion diagrams for Poincaré, Kelvin and Rossby waves are presented. Special attention is given to the nonlinear Rossby waves on a β-plane ocean. Based on the perturbation analysis, it is shown that the nonlinear evolution equation for the wave amplitude satisfies a modified nonlinear Schrödinger equation. The solution of this equation represents solitary waves in a dispersive medium. In other words, the envelope of the amplitude of the waves has a soliton structure and these envelope solitons propagate with the group velocity of the Rossby waves. Finally, a nonlinear analytical model is presented for long Rossby waves in a meridional channel with weak shear. A new nonlinear wave equation for the amplitude of large Rossby waves is derived in a region where fluid flows over the recirculation core. It is shown that the governing amplitude equations for the inner and outer zones are both KdV type, where weak nonlinearity is balanced by weak dispersion. In the inner zone, the nonlinear amplitude equation has a new term proportional to the 3/2 power of the difference between the wave amplitude and the critical amplitude, and this term occurs to account for a nonlinearity due to the flow over the vortex core. The solution of the amplitude equations with the linear shear flow represents the solitary waves. The present study deals with the lowest mode (n=1) analysis. An extension of the higher modes (n?2) of this work will be made in a subsequent paper.     

On conservative spatial discretizations of the barotropic quasi-gasdynamic system of equations with a potential body force

A multidimensional barotropic quasi-gasdynamic system of equations in the form of mass and momentum conservation laws with a general gas equation of state p = p(ρ) with p′(ρ) > 0 and a potential body force is considered. For this system, two new symmetric spatial discretizations on nonuniform rectangular grids are constructed (in which the density and velocity are defined on the basic grid, while the components of the regularized mass flux and the viscous stress tensor are defined on staggered grids). These discretizations involve nonstandard approximations for ?p(ρ), div(ρu), and ρ. As a result, a discrete total mass conservation law and a discrete energy inequality guaranteeing that the total energy does not grow with time can be derived. Importantly, these discretizations have the additional property of being well-balanced for equilibrium solutions. Another conservative discretization is discussed in which all mass flux components and viscous stresses are defined on the same grid. For the simpler barotropic quasi-hydrodynamic system of equations, the corresponding simplifications of the constructed discretizations have similar properties.     

On the dynamics of a fluid-particle interaction model: The bubbling regime

This article deals with the issues of global-in-time existence and asymptotic analysis of a fluid-particle interaction model in the so-called bubbling regime. The mixture occupies the physical space Ω⊂R³ which may be unbounded. The system under investigation describes the evolution of particles dispersed in a viscous compressible fluid and is expressed through the conservation of fluid mass, the balance of momentum and the balance of particle density often referred as the Smoluchowski equation. The coupling between the dispersed and dense phases is obtained through the drag forces that the fluid and the particles exert mutually by the action-reaction principle. We show that solutions exist globally in time under reasonable physical assumptions on the initial data, the physical domain, and the external potential. Furthermore, we prove the large-time stabilization of the system towards a unique stationary state fully determined by the masses of the initial density of particles and fluid and the external potential.     

On turbulence in fractal porous media

We examine three fundamental equations governing turbulence of an incompressible Newtonian fluid in a fractal porous medium: continuity, linear momentum balance and energy balance. We find that the Reynolds stress is modified when a local, rather than an integral, balance law is considered. The heat flux is modified from its classical form when either the integral or local form of the energy density balance law is studied, but the energy density is always unchanged. The modifications of Reynolds stress and heat flux are expressed directly in terms of the resolution length scale, the fractal dimension of mass distribution and the fractal dimension of a fractal’s surface. When both fractal dimensions become integer (respectively 3 and 2), classical equations are recovered.      

Mathematical and physical simulation of the interaction between a gas jet and a liquid free surface

In this work a two phase 3D mathematical model was developed using the volume of fluid (VOF) algorithm, which is able to accurately describe the cavity geometry and size as well as the liquid flow patterns created when a gas jet that impinges on a liquid free surface. These phenomena are commonly found in steelmaking operations such as in the Electric Arc Furnace (EAF) and the Basic Oxygen Furnace (BOF) where oxygen jets impinge on a steel bath and they control heat, momentum and mass transfer. The model was successfully validated with measurements made on a physical model through velocity fields obtained by Particle Image Velocimetry (PIV) and high speed camera images of the cavity. Agreement between model predictions and experimental measurements is excellent in both x-velocity component of the liquid and cavity sizes. The cavity formed in the liquid by the impinging jet depends on a force balance at the free surface where the inertial force of the jet governs this phenomena, while the liquid circulation depends on also the jet inertial force of the jet, but its angle plays an important role, being the lowest angle the best choice to shear the bath and promote stronger circulation and better mixing in the liquid.     

Solitary waves and a stability analysis of an equation of short and long dispersive waves

A universal model for the interaction of long nonlinear waves and packets of short waves with long linear carrier waves is given by a system in which an equation of Korteweg–de Vries (KdV) type is coupled to an equation of nonlinear Schrödinger (NLS) type. The system has solutions of steady form in which one component is like a solitary-wave solution of the KdV equation and the other component is like a ground-state solution of the NLS equation. We study the stability of solitary-wave solutions to an equation of short and long waves by using variational methods based on the use of energy–momentum functionals and the techniques of convexity type. We use the concentration compactness method to prove the existence of solitary waves. We prove that the stability of solitary waves is determined by the convexity or concavity of a function of the wave speed.     

Multiphase averaging and the inverse spectral solution of the Korteweg—de Vries equation

Inverse spectral theory is used to prescribe and study equations for the slow modulations of N-phase wave trains for the Korteweg-de Vries (KdV) equation. An invariant representation of the modulational equations is deduced. This representation depends upon certain differentials on a Riemann surface. When evaluated near ∞ on the surface, the invariant representation reduces to averaged conservations laws; when evaluated near the branch points, the representation shows that the simple eigenvalues provide Riemann invariants for the modulational equations. Integrals of the invariant representation over certain cycles on the Riemann surface yield “conservation of waves.” Explicit formulas for the characteristic speeds of the modulational equations are derived. These results generalize known results for a single-phase traveling wave, and indicate that complete integrability can induce enough structure into the modulational equations to diagonalize (in the sense of Riemann invariants) their first-order terms.     

Numerical approximation for a Baer-Nunziato model of two-phase flows

We present a well-balanced numerical scheme for approximating the solution of the Baer-Nunziato model of two-phase flows by balancing the source terms and discretizing the compaction dynamics equation. First, the system is transformed into a new one of three subsystems: the first subsystem consists of the balance laws in the gas phase, the second subsystem consists of the conservation law of the mass in the solid phase and the conservation law of the momentum of the mixture, and the compaction dynamic equation is considered as the third subsystem. In the first subsystem, stationary waves are used to build up a well-balanced scheme which can capture equilibrium states. The second subsystem is of conservative form and thus can be numerically treated in a standard way. For the third subsystem, the fact that the solid velocity is constant across the solid contact suggests us to compose the technique of the Engquist-Osher scheme. We show that our scheme is capable of capturing exactly equilibrium states. Moreover, numerical tests show the convergence of approximate solutions to the exact solution.     

Systematic derivation of an asymptotic model for the dynamics of curved viscous fibers

This paper presents a slender body theory for the dynamics of a curved inertial viscous Newtonian fiber. Neglecting surface tension and temperature dependence, the fiber flow is modeled as a three‐dimensional free boundary value problem in terms of instationary incompressible Navier–Stokes equations. From regular asymptotic expansions in powers of the slenderness parameter, leading‐order balance laws for mass (cross‐section) and momentum are derived that combine the unrestricted motion of the fiber centerline with the inner viscous transport. The physically reasonable form of the one‐dimensional fiber model results thereby from the introduction of the intrinsic velocity that characterizes the convective terms. For the numerical investigation of the viscous, gravitational and rotational effects on the fiber dynamics, a finite volume approach on a staggered grid with implicit upwind flux discretization is applied. Copyright © 2007 John Wiley & Sons, Ltd.     

Global solution to the full one-dimensional Frémond model for shape memory alloys

In this paper, we prove the existence and uniqueness of the solution to the one-dimensional initial-boundary value problem resulting from the Frémond thermomechanical model of structural phase transitions in shape memory materials. In this model, the free energy is assumed to depend on temperature, macroscopic deformation and phase fractions. The resulting equilibrium equations are the balance laws of (linear) momentum and energy, coupled with an evolution variational inequality for the phase fractions. Fourth-order regularizing terms in the quasi-stationary momentum balance equation are not necessary, and, as far as we know for the first time, all the non-linear terms of the energy balance equation are taken into account.     

Similarity solution for a cylindrical shock wave in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux

The propagation of shock waves in a rotational axisymmetric dusty gas with heat conduction and radiation heat flux, which has a variable azimuthally fluid velocity together with a variable axial fluid velocity, is investigated. The dusty gas is assumed to be a mixture of non-ideal (or perfect) gas and small solid particles, in which solid particles are continuously distributed. It is assumed that the equilibrium flow-condition is maintained and variable energy input is continuously supplied by the piston (or inner expanding surface). The fluid velocities in the ambient medium are assume to be vary and obey power laws. The density of the ambient medium is assumed to be constant, the heat conduction is express in terms of Fourier’s law and the radiation is considered to be of the diffusion type for an optically thick grey gas model. The thermal conductivity K and the absorption coefficient α_R are assumed to vary with temperature and density. In order to obtain the similarity solutions the angular velocity of the ambient medium is assume to be decreasing as the distance from the axis increases. The effects of the variation of the heat transfer parameter and non-idealness of the gas in the mixture are investigated. The effects of an increase in (i) the mass concentration of solid particles in the mixture and (ii) the ratio of the density of solid particles to the initial density of the gas on the flow variables are also investigated.     

Numerical study of transient free convective mass transfer in a Walters-B viscoelastic flow with wall suction

A transient model for the free convective, nonlinear, steady, laminar flow and mass transfer in a viscoelastic fluid from a vertical porous plate is presented. The Walters-B liquid model is employed which introduces supplementary terms into the momentum conservation equation. The transformed conservation equations are solved using the finite difference method (FDM). The influence of viscoelasticity parameter (Γ), species Grashof number (Gc), Schmidt number (Sc), distance (Y) and time (t) on the velocity (U) and also concentration distribution (C) is studied graphically. Velocity is found to increase with a rise in viscoelasticity parameter (Γ) with both time and distances close to the plate surface. An increase in Schmidt number is observed to significantly decrease both velocity and concentration in time and also with separation from the plate. Increasing species Grashof number boosts the flow velocity through all time and causes a significant rise primarily near the plate surface. The study has applications in polymer materials processing.     

Hydrodynamic and thermal slip flow boundary layers over a flat plate with constant heat flux boundary condition

In this paper the boundary layer flow over a flat plat with slip flow and constant heat flux surface condition is studied. Because the plate surface temperature varies along the x direction, the momentum and energy equations are coupled due to the presence of the temperature gradient along the plate surface. This coupling, which is due to the presence of the thermal jump term in Maxwell slip condition, renders the momentum and energy equations non-similar. As a preliminary study, this paper ignores this coupling due to thermal jump condition so that the self-similar nature of the equations is preserved. Even this fundamental problem for the case of a constant heat flux boundary condition has remained unexplored in the literature. It was therefore chosen for study in this paper. For the hydrodynamic boundary layer, velocity and shear stress distributions are presented for a range of values of the parameter characterizing the slip flow. This slip parameter is a function of the local Reynolds number, the local Knudsen number, and the tangential momentum accommodation coefficient representing the fraction of the molecules reflected diffusively at the surface. As the slip parameter increases, the slip velocity increases and the wall shear stress decreases. These results confirm the conclusions reached in other recent studies. The energy equation is solved to determine the temperature distribution in the thermal boundary layer for a range of values for both the slip parameter as well as the fluid Prandtl number. The increase in Prandtl number and/or the slip parameter reduces the dimensionless surface temperature. The actual surface temperature at any location of x is a function of the local Knudsen number, the local Reynolds number, the momentum accommodation coefficient, Prandtl number, other flow properties, and the applied heat flux.     

Exact solutions of nonlinear dispersive K(m, n) model with variable coefficients

The variable-coefficient Korteweg-de Vries (KdV) equation with additional terms contributed from the inhomogeneity in the axial direction and the strong transverse confinement of the condense was presented to describe the dynamics of nonlinear excitations in trapped quasi-one-dimensional Bose-Einstein condensates with repulsive atom-atom interactions. To understand the role of nonlinear dispersion in this variable-coefficient model, we introduce and study a new variable-coefficient KdV with nonlinear dispersion (called vc-K(m, n) equation). With the aid of symbolic computation, we obtain its compacton-like solutions and solitary pattern-like solutions. Moreover, we also present some conservation laws for both vc-K⁺(n, n) equation and vc-K⁻(n, n) equation.     

Cnoidal wave trains and solitary waves in a dissipation-modified Korteweg-de Vries equation

A generalization of the Korteweg-de Vries equation incorporating an energy input-output balance, hence a dissipation-modified KdV equation is considered. The equation is relevant to describe, for instance, nonlinear Marangoni-Bénard oscillatory instability in a liquid layer heated from above. Cnoidal waves and solitary waves of this equation are obtained both asymptotically and numerically.