Steady solutions for feeder dykes in a density-stratified lithosphere

1991 ◽  
Vol 107 (2) ◽  
pp. 233-242 ◽  
Author(s):  
John R. Lister
Keyword(s):  
1997 ◽  
Vol 335 ◽  
pp. 233-259 ◽  
Author(s):  
P. W. DUCK ◽  
M. R. FOSTER ◽  
R. E. HEWITT

In this paper we consider the boundary layer that forms on the sloping walls of a rotating container (notably a conical container), filled with a stratified fluid, when flow conditions are changed abruptly from some initial (uniform) state. The structure of the solution valid away from the cone apex is derived, and it is shown that a similarity-type solution is appropriate. This system, which is inherently nonlinear in nature, is solved numerically for several flow regimes, and the results reveal a number of interesting and diverse features.In one case, a steady state is attained at large times inside the boundary layer. In a second case, a finite-time singularity occurs, which is fully analysed. A third scenario involves a double boundary-layer structure developing at large times, most significantly including an outer region that grows in thickness as the square-root of time.We also consider directly the nonlinear fully steady solutions to the problem, and map out in parameter space the likely ultimate flow behaviour. Intriguingly, we find cases where, when the rotation rate of the container is equal to that of the main body of the fluid, an alternative nonlinear state is preferred, rather than the trivial (uniform) solution.Finally, utilizing Laplace transforms, we re-investigate the linear initial-value problem for small differential spin-up studied by MacCready & Rhines (1991), recovering the growing-layer solution they found. However, in contrast to earlier work, we find a critical value of the buoyancy parameter beyond which the solution grows exponentially in time, consistent with our nonlinear results.


2005 ◽  
Vol 47 (2) ◽  
pp. 185-202 ◽  
Author(s):  
T. E. Stokes ◽  
G. C. Hocking ◽  
L. K. Forbes

AbstractThe unsteady axisymmetric withdrawal from a fluid with a free surface through a point sink is considered. Results both with and without surface tension are included and placed in context with previous work. The results indicate that there are two critical values of withdrawal rate at which the surface is drawn directly into the outlet, one after flow initiation and the other after the flow has been established. It is shown that the larger of these values corresponds to the point at which steady solutions no longer exist.


2016 ◽  
Vol 26 (3/4) ◽  
pp. 879-896 ◽  
Author(s):  
D. Andrew S. Rees

Purpose – The purpose of this paper is to determine the manner in which a yield stress fluid begins convecting when it saturates a porous medium. A sidewall-heated rectangular cavity is selected as the testbed for this pioneering work. Design/methodology/approach – Steady solutions are obtained using a second order accurate finite difference method, line relaxation based on the Gauss-Seidel smoother, a Full Approximation Scheme multigrid algorithm with V-cycling and a regularization of the Darcy-Bingham model to smooth the piecewise linear relation between the Darcy flux and the applied body forces. Findings – While Newtonian fluids always convect whenever the Darcy-Rayleigh number is nonzero, Bingham fluids are found to convect only when the Darcy-Rayleigh number exceeds a value which is linearly dependent on both the Rees-Bingham number and the overall perimeter of the rectangular cavity. Stagnation is always found in the centre of the cavity and in regions close to the four corners. Care must be taken over the selection of the regularization constant. Research limitations/implications – The Darcy-Rayleigh number is restricted to values which are at or below 200. Originality/value – This is the first investigation of the effect of yield stress on nonlinear convection in porous media.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Thida Pongsanguansin ◽  
Montri Maleewong ◽  
Khamron Mekchay

A well-balanced scheme with total variation diminishing Runge-Kutta discontinuous Galerkin (TVD-RK DG) method for solving shallow water equations is presented. Generally, the flux function at cell interface in the TVD-RK DG scheme is approximated by using the Harten-Lax-van Leer (HLL) method. Here, we apply the weighted average flux (WAF) which is higher order approximation instead of using the HLL in the TVD-RK DG method. The consistency property is shown. The modified well-balanced technique for flux gradient and source terms under the WAF approximations is developed. The accuracy of numerical solutions is demonstrated by simulating dam-break flows with the flat bottom. The steady solutions with shock can be captured correctly without spurious oscillations near the shock front. This presents the other flux approximations in the TVD-RK DG method for shallow water simulations.


2008 ◽  
Vol 617 ◽  
pp. 283-299 ◽  
Author(s):  
E. S. BENILOV ◽  
V. S. ZUBKOV

We consider an infinite plate being withdrawn (at an angle α to the horizontal, with a constant velocity U) from an infinite pool of viscous liquid. Assuming that the effects of inertia and surface tension are weak, Derjaguin (C. R. Dokl. Acad. Sci. URSS, vol. 39, 1943, p. 13.) conjectured that the ‘load’ l, i.e. the thickness of the liquid film clinging to the plate, is l=(μU/ρgsinα)1/2, where ρ and μ are the liquid's density and viscosity, and g is the acceleration due to gravity.In the present work, the above formula is derived from the Stokes equations in the limit of small slopes of the plate (without this assumption, the formula is invalid). It is shown that the problem has infinitely many steady solutions, all of which are stable – but only one of these corresponds to Derjaguin's formula. This particular steady solution can only be singled out by matching it to a self-similar solution describing the non-steady part of the film between the pool and the film's ‘tip’.Even though the near-pool region where the steady state has been established expands with time, the upper, non-steady part of the film (with its thickness decreasing towards the tip) expands faster and, thus, occupies a larger portion of the plate. As a result, the mean thickness of the film is 1.5 times smaller than the load.


1998 ◽  
Vol 35 (6) ◽  
pp. 2329-2349 ◽  
Author(s):  
Gunilla Kreiss ◽  
Heinz-Otto Kreiss
Keyword(s):  

2011 ◽  
Vol 113 (3) ◽  
pp. 289-303 ◽  
Author(s):  
L. Margheriti ◽  
M. P. Speciale

2010 ◽  
Vol 7 (2) ◽  
pp. 2347-2371 ◽  
Author(s):  
◽  
◽  
◽  

Abstract. The stream depletion rate (SDR) associated with pumping from vertical wells located in an aquifer is commonly estimated, where a large drawdown near the well may, however, be produced. In this paper, the solution is first developed for describing the groundwater flow associated with a point source in a confined aquifer near a stream. Based on the principle of superposition, analytical solutions for horizontal and slanted wells are then developed by integrating the point source solution along the well axis. The solutions can be simplified to quasi-steady solutions by neglecting the exponential terms to describe the late-time drawdown, which can provide useful information in designing horizontal well location and length. The direction of the well axis can be determined from the best SDR subject to the drawdown constraint. It is found that hydraulic conductivity in the direction perpendicular to the stream plays a crucial role in influencing the time required for reaching quasi-steady SDR. In addition, the effects of the well length as well as the distance between the well and stream on the SDR are also examined.


Sign in / Sign up

Export Citation Format

Share Document