Models of non-Boussinesq lock-exchange flow

2011 ◽  
Vol 675 ◽  
pp. 1-26 ◽  
Author(s):  
R. ROTUNNO ◽  
J. B. KLEMP ◽  
G. H. BRYAN ◽  
D. J. MURAKI

Nearly all analytical models of lock-exchange flow are based on the shallow-water approximation. Since the latter approximation fails at the leading edges of the mutually intruding fluids of lock-exchange flow, solutions to the shallow-water equations can be obtained only through the specification of front conditions. In the present paper, analytic solutions to the shallow-water equations for non-Boussinesq lock-exchange flow are given for front conditions deriving from free-boundary arguments. Analytic solutions are also derived for other proposed front conditions – conditions which appear to the shallow-water system as forced boundary conditions. Both solutions to the shallow-water equations are compared with the numerical solutions of the Navier–Stokes equations and a mixture of successes and failures is recorded. The apparent success of some aspects of the forced solutions of the shallow-water equations, together with the fact that in a real fluid the density interface is a free boundary, shows the need for an improved theory of lock-exchange flow taking into account non-hydrostatic effects for density interfaces intersecting rigid boundaries.

2010 ◽  
Vol 20 (07) ◽  
pp. 1049-1087 ◽  
Author(s):  
BORIS HASPOT

In this paper, we consider the compressible Navier–Stokes equation with density-dependent viscosity coefficients and a term of capillarity introduced formally by van der Waals in Ref. 51. This model includes at the same time the barotropic Navier–Stokes equations with variable viscosity coefficients, shallow-water system and the model introduced by Rohde in Ref. 46. We first study the well-posedness of the model in critical regularity spaces with respect to the scaling of the associated equations. In a functional setting as close as possible to the physical energy spaces, we prove global existence of solutions close to a stable equilibrium, and local in time existence of solutions with general initial data. Uniqueness is also obtained.


2012 ◽  
Vol 14 (3) ◽  
pp. 772-783 ◽  
Author(s):  
Ljubomir Budinski

In order to improve efficiency and accuracy, while maintaining an ease of modeling flows with the lattice Boltzmann approach in domains having complex geometry, a method for modeling equations of 2D flow in curvilinear coordinates has been developed. Both the transformed shallow water equations and the transformed 2D Navier-Stokes equations in the horizontal plane were synchronized with the equilibrium distribution function and the force term in the rectangular lattice. Since the solution of these equations takes place in the classical rectangular lattice environment, boundary conditions are modeled in the standard form of already existing simple methods (bounce-back), not requiring any additional functions. Owing to this and to the fact that the proposed method ensures a more accurate fitting of equations, even to domains of interest having complex geometry, the accuracy of solution is significantly increased, while the simplicity of the standard lattice Boltzmann approach is maintained. For the shallow water equations transformed in curvilinear coordinates, the proposed procedure is verified in three different hydraulic problems, all characterized by complex geometry.


Author(s):  
Hiroshi Dan ◽  
◽  
Hiroshi Kanayama

In this paper, viscous shallow-water equations are derived from three-dimensional Navier-Stokes equations under the hydrostatic assumption. The viscous shallow-water equations are approximated by the finite element method based on our numerical scheme developed in 1978. This approach is used to simulate a tsunami in Hakata Bay. Results show a reasonable estimate of the tsunami arrival time.


1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.


2021 ◽  
Vol 153 (A2) ◽  
Author(s):  
Q Yang ◽  
W Qiu

Slamming forces on 2D and 3D bodies have been computed based on a CIP method. The highly nonlinear water entry problem governed by the Navier-Stokes equations was solved by a CIP based finite difference method on a fixed Cartesian grid. In the computation, a compact upwind scheme was employed for the advection calculations and a pressure-based algorithm was applied to treat the multiple phases. The free surface and the body boundaries were captured using density functions. For the pressure calculation, a Poisson-type equation was solved at each time step by the conjugate gradient iterative method. Validation studies were carried out for 2D wedges with various deadrise angles ranging from 0 to 60 degrees at constant vertical velocity. In the cases of wedges with small deadrise angles, the compressibility of air between the bottom of the wedge and the free surface was modelled. Studies were also extended to 3D bodies, such as a sphere, a cylinder and a catamaran, entering calm water. Computed pressures, free surface elevations and hydrodynamic forces were compared with experimental data and the numerical solutions by other methods.


1965 ◽  
Vol 69 (658) ◽  
pp. 714-718 ◽  
Author(s):  
Ronald D. Mills

The Navier-Stokes equations are solved iteratively on a small digital computer for the class of flows generated within a rectangular “cavity” by a surface passing over its open end. Solutions are presented for depth/breadth ratios ƛ=0.5 (shallow), 10 (square), 20 (deep) and Reynolds number 100. Flow photographs ore obtained which largely confirm the predicted flows. The theoretical velocity profiles and pressure distributions through the centre of the vortex in the square cavity are calculated.In an appendix an improved finite difference formula is given for the vorticity generated at a moving boundary.Since Thorn began his pioneering work some thirty-five years ago the number of numerical solutions which have been obtained for the equations of incompressible viscous fluid motion remains small (see bibliographies of Thom and Apelt, Fromm). The known solutions are principally for steady streaming flows, although two methods have now been used with success for non-steady flows (Payne jets and Fromm flow past obstacles). By contrast this paper is concerned with the class of closed flows generated in a rectangular region of varying depth/breadth ratio by a surface passing over an open end. This problem has been considered for a number of reasons.


Sign in / Sign up

Export Citation Format

Share Document