SUMMARYThe paper describes a new approach to approximating the convection term found in typical steady-state transport equations. A polynomial-based discretization scheme is constructed around a technique called 'curvature compensation'; the resultant curvature-compensated convective transport approximation is essentially third-order accurate in regions of the solution domain where the concept of order is meaningful. In addition, in linear scalar transport problems it preserves the boundedness of solutions. Sharp changes in gradient in the dependent variable are handled particularly well. But above all, the scheme, when used in conjunction with an AD1 pentadiagonal solver, is easy to implement with relatively low computational cost, representing an effective algorithm for the simulation of multi-dimensional fluid flows. Two linear test problems, for the case of transport by pure convection, are employed in order to assess the merit of the method.

No abstract

No abstract

Article:Gaskell, P.H., Jimack, P.K., Sellier, M. et A range of two-and three-dimensional problems is explored featuring the gravitydriven flow of a continuous thin liquid film over a non-porous inclined flat surface containing well-defined topography. These are analysed principally within the framework of the lubrication approximation, where accurate numerical solution of the governing nonlinear equations is achieved using an efficient multigrid solver.Results for flow over one-dimensional steep-sided topographies are shown to be in very good agreement with previously reported data. The accuracy of the lubrication approximation in the context of such topographies is assessed and quantified by comparison with finite element solutions of the full Navier-Stokes equations, and results support the consensus that lubrication theory provides an accurate description of these flows even when its inherent assumptions are not strictly satisfied. The NavierStokes solutions also illustrate the effect of inertia on the capillary ridge/trough and the two-dimensional flow structures caused by steep topography.Solutions obtained for flow over localized topography are shown to be in excellent agreement with the recent experimental results of Decré & Baret (2003) for the motion of thin water films over finite trenches. The spread of the 'bow wave', as measured by the positions of spanwise local extrema in free-surface height, is shown to be well-represented both upstream and downstream of the topography by an inverse hyperbolic cosine function.An explanation, in terms of local flow rate, is given for the presence of the 'downstream surge' following square trenches, and its evolution as trench aspect ratio is increased is discussed. Unlike the upstream capillary ridge, this feature cannot be completely suppressed by increasing the normal component of gravity. The linearity of free-surface response to topographies is explored by superposition of the free surfaces corresponding to two 'equal-but-opposite' topographies. Results confirm the findings of Decré & Baret (2003) that, under the conditions considered, the responses behave in a near-linear fashion.

No abstract

No abstract

Three mathematical models are developed for meniscus roll coating in which there is steady flow of a Newtonian fluid in the narrow gap, or nip, between two contrarotating rolls in the absence of body forces.The zero flux model predicts a constant pressure gradient within the central core and two eddies, each with an inner structure, in qualitative agreement with observation. The small flux model takes account of a small inlet flux and employs the lubrication approximation to represent fluid velocity as a combination of Couette and Poiseuille flows. Results show that the meniscus coating regime is characterized by small flow rates (λ [Lt ] 1) and a sub-ambient pressure field generated by capillary action at the upstream meniscus. Such flows are found to exist for small modified capillary number, Ca(R/H0)1/2 [lsim ] 0.15, where Ca and R/H0 represent capillary number and the radius to semi-gap ratio, respectively.A third model incorporates the full effects of curved menisci and nonlinear free surface boundary conditions. The presence of a dynamic contact line, adjacent to the web on the upper roll, requires the imposition of an apparent contact angle and slip length. Numerical solutions for the velocity and pressure fields over the entire domain are obtained using the finite element method. Results are in accord with experimental observations that the flow domain consists of two large eddies and fluid transfer jets or ‘snakes’. Furthermore, the numerical results show that the sub-structure of each large eddy consists of a separatrix with one saddle point, two sub-eddies with centres, and an outer recirculation.Finally finite element solutions in tandem with lubrication analysis establish the existence of three critical flow rates corresponding to a transformation of the pressure field, the emergence of a ‘secondary snake’ (another fluid transfer jet) and the disappearance of a primary snake.

Published paperAbstract A range of problems is investigated, involving the gravity-driven inertial flow of a thin viscous liquid film over a planar surface containing topographical features, modelled via a depth-averaged form of the governing unsteady Navier-Stokes equations. The discrete analogue of the resulting coupled equation set, employing a staggered mesh arrangement for the dependent variables, is solved accurately using an efficient Full Approximation Storage (FAS) algorithm and Full Multigrid (FMG) technique together with adaptive timestepping and proper treatment of the nonlinear convective terms. A unique, comprehensive set of results is presented for both one-and two-dimensional topographical features, and errors quantified via detailed comparisons drawn with complementary experimental data and predictions from finite element analyses where they exist. It is found in the case of one-dimensional (spanwise) topography that for small Reynolds number and shallow/short features the depth-averaged form produces results that are in close agreement with corresponding finite element solutions of the full free-surface problem. For the case of flow over two-dimensional (localised) topography the free-surface disturbance observed is influenced significantly by the presence of inertia. It leads, as in the case of spanwise topography, to an increase in the magnitude and severity of the capillary ridge/trough patterns which form.