A low‐dimensional spectral approach for transient axisymmetric free‐surface flow inside thin cavities of arbitrary shape

A low‐dimensional spectral method is used to solve the transient axisymmetric free surface flow inside thin cavities of arbitrary shape. The flow field is obtained on the basis of the lubrication equations, which are expanded in terms of orthonormal functions over the cavity gap. The formulation accounts for nonlinearities stemming from inertia and front location. The work is of close relevance to the filling stage during die casting, and injection molding, or the flow inside annular (extrusion) dies. Both flows under an imposed flow rate, and an imposed pressure at the cavity entrance are examined. The influence of inertia, aspect ratio, gravity, and wall geometry on the evolution of the front, flow rate, and pressure is assessed particularly in the early stage of flow, when a temporal behavior of the ‘boundary‐layer’ type develops. The multiple‐scale method is applied to obtain an approximate solution at small Reynolds number, Re. Comparison with the exact (numerical) solution indicates a wide range of validity for the multiple‐scale approach, including the moderately small Re range. Copyright © 2002 John Wiley & Sons, Ltd.     

A hybrid spectral/finite‐difference to transient free‐surface flow inside thin symmetric cavities of longitudinally varying thickness

The lubrication theory is extended for transient free‐surface flow of a viscous fluid inside three‐dimensional symmetric thin cavities of thickness that varies in the flow direction. The problem is first formulated for a cavity of arbitrary shape. The moving domain is mapped onto a rectangular domain at each time step of the computation. The pressure, which in this case is governed by the modified Laplace's equation, is expanded in a Fourier series in the spanwise direction. The expansion coefficients are obtained using the finite‐difference method. Only a few modes are usually needed to secure convergence. The flow behaviour is strongly influenced by the cavity thickness. The flows inside a straight, contracting, expanding, and modulated cavities are examined. It is found that while the evolution of the front is always monotonic with time, that of the velocity at the front can be oscillatory if the degree of contraction of the cavity (whether modulated or not) is significant. The velocity of the contact point along the lateral walls is usually larger than that at the front, leading to the straightening of the front. However, for modulated cavities, the front may advance at a faster rate, leading to its own undulation. Copyright © 2005 John Wiley & Sons, Ltd.     

A Lagrangian boundary element approach to transient three‐dimensional free surface flow in thin cavities

The lubrication theory is extended for transient free‐surface flow of a viscous fluid inside a three‐dimensional thin cavity. The problem is closely related to the filling stage during the injection molding process. The pressure, which in this case is governed by the Laplace's equation, is determined using the boundary element method. A fully Lagrangian approach is implemented for the tracking of the evolving free surface. The domain of computation is the projection of the physical domain onto the (x, y) plane. This approach is valid for simple and complex cavities as illustrated for the cases of a flat plate and a curved plate. It is found that the flow behavior is strongly influenced by the shape of the initial fluid domain, the shape of the cavity, and inlet flow pressure. Copyright © 2001 John Wiley & Sons, Ltd.     

Three‐dimensional finite element analysis of transient fluid flow with free‐surface using marker surface method and adaptive grid refinement

For three‐dimensional finite element analysis of transient fluid flow with free‐surface, a new marker surface method is proposed, in which the fluid flow is represented by the marker surface composed of marker elements instead of marker particles used in the marker particle method. This also involves an adaptive grid that is created under a new criterion of element categorization of filling states and the locations in the total region at each time step. The marker surface is used in order to represent the free‐surface accurately, as well as to decrease the memory and computation time, and to effectively display the predicted three‐dimensional free‐surface. By using the adaptive grid in which the elements, finer than those in internal and external regions, are distributed in the surface region through refinement and coarsening procedures, the analysis of three‐dimensional transient fluid flow with free‐surface is achieved more efficiently. Through three‐dimensional analysis of two kinds of problems using several grids, the efficiency of the proposed marker surface method and the adaptive grid are verified. Copyright © 1999 John Wiley & Sons, Ltd.     

A low‐dimensional spectral approach for transient free‐surface flow inside thin cavities of symmetric shape

The lubrication theory is extended for transient free‐surface flow of a viscous fluid inside three‐dimensional cavities of general symmetric shape but of small thickness. The problem is closely related to the filling stage during the injection molding process. The moving domain is mapped onto a rectangular domain at each time step of the computation. A modified pressure is introduced, which in this case is governed by the Laplace's equation, and it is expanded in a Fourier series along the flow direction. The expansion coefficients are obtained using the finite‐difference method. This approach is valid for simple and complex cavities as illustrated for the cases of a flat plate and a curved plate. Only a few modes are needed to secure convergence in general. It is found that the flow behaviour is strongly influenced by the shape of the initial fluid domain, the shape of the cavity, cavity thickness, and the inlet flow. Copyright © 2002 John Wiley & Sons, Ltd.     

Three‐dimensional modeling of non‐hydrostatic free‐surface flows on unstructured grids

A three‐dimensional numerical model is presented for the simulation of unsteady non‐hydrostatic shallow water flows on unstructured grids using the finite volume method. The free surface variations are modeled by a characteristics‐based scheme, which simulates sub‐critical and super‐critical flows. Three‐dimensional velocity components are considered in a collocated arrangement with a σ‐coordinate system. A special treatment of the pressure term is developed to avoid the water surface oscillations. Convective and diffusive terms are approximated explicitly, and an implicit discretization is used for the pressure term to ensure exact mass conservation. The unstructured grid in the horizontal direction and the σ coordinate in the vertical direction facilitate the use of the model in complicated geometries. Solution of the non‐hydrostatic equations enables the model to simulate short‐period waves and vertically circulating flows. Copyright © 2015 John Wiley & Sons, Ltd.     

The solution of two‐dimensional free‐surface problems using automatic mesh generation

A new method is described for the iterative solution of two‐dimensional free‐surface problems, with arbitrary initial geometries, in which the interior of the domain is represented by an unstructured, triangular Eulerian mesh and the free surface is represented directly by the piecewise‐quadratic edges of the isoparametric quadratic‐velocity, linear‐pressure Taylor–Hood elements. At each time step, the motion of the free surface is computed explicitly using the current velocity field and, once the new free‐surface location has been found, the interior nodes of the mesh are repositioned using a continuous deformation model that preserves the original connectivity. In the event that the interior of the domain must be completely remeshed, a standard Delaunay triangulation algorithm is used, which leaves the initial boundary discretisation unchanged. The algorithm is validated via the benchmark viscous flow problem of the coalescence of two infinite cylinders of equal radius, in which the motion is due entirely to the action of capillary forces on the free surface. This problem has been selected for a variety of reasons: the initial and final (steady state) geometries differ considerably; in the passage from the former to the latter, large free‐surface curvatures—requiring accurate modelling—are encountered; an analytical solution is known for the location of the free surface; there exists a large body of literature on alternative numerical simulations. A novel feature of the present work is its geometric generality and robustness; it does not require a priori knowledge of either the evolving domain geometry or the solution contained therein. Copyright © 1999 John Wiley & Sons, Ltd.     

Transient free‐surface flow of a viscoelastic fluid in a narrow channel

The interplay between inertia and elasticity is examined for transient free‐surface flow inside a narrow channel. The lubrication theory is extended for the flow of viscoelastic fluids of the Oldroyd‐B type (consisting of a Newtonian solvent and a polymeric solute). While the general formulation accounts for non‐linearities stemming from inertia effects in the momentum conservation equation, and the upper‐convected terms in the constitutive equation, only the front movement contributes to non‐linear coupling for a flow inside a straight channel. In this case, it is possible to implement a spectral representation in the depthwise direction for the velocity and stress. The evolution of the flow field is obtained locally, but the front movement is captured only in the mean sense. The influence of inertia, elasticity and viscosity ratio is examined for pressure‐induced flow. The front appears to progress monotonically with time. However, the velocity and stress exhibit typically a strong overshoot upon inception, accompanied by a plug‐flow behaviour in the channel core. The flow intensity eventually diminishes with time, tending asymptotically to Poiseuille conditions. For highly elastic liquids the front movement becomes oscillatory, experiencing strong deceleration periodically. A multiple‐scale solution is obtained for fluids with no inertia and small elasticity. Comparison with the exact (numerical) solution indicates a wide range of validity for the analytical result. Copyright © 2004 John Wiley & Sons, Ltd.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

An adaptive Lagrangian boundary element approach for three‐dimensional transient free‐surface Stokes flow as applied to extrusion,thermoforming, and rheometry

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional simulation of confined free‐surface Stokes flow. The method is stable as it includes remeshing capabilities of the deforming free surface and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. Several flow problems are presented to illustrate the utility of the approach, particularly as encountered in polymer processing and rheology. These problems illustrate the transient nature of the flow during the processes of extrusion and thermoforming, the elongation of a fluid sample in an extensional rheometer, and the coating of a sphere. Surface tension effects are also explored. Copyright © 2001 John Wiley & Sons, Ltd.     

A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 John Wiley & Sons, Ltd.     

Shape optimization approaches to free‐surface problems

Three different reformulations of a free‐surface problem as shape optimization problems are considered. These give rise to three different cost functionals that apparently have not been exploited in literature. The shape derivatives of the cost functionals are explicitly determined. The gradient information is combined with the boundary variation method in a preconditioned steepest descent algorithm to solve the shape optimization problems. Numerical results that compare the performance of the proposed cost functionals are presented. Copyright © 2014 John Wiley & Sons, Ltd.     

An implicit two‐dimensional non‐hydrostatic model for free‐surface flows

An implicit finite volume model in sigma coordinate system is developed to simulate two‐dimensional (2D) vertical free surface flows, deploying a non‐hydrostatic pressure distribution. The algorithm is based on a projection method which solves the complete 2D Navier–Stokes equations in two steps. First the pressure term in the momentum equations is excluded and the resultant advection–diffusion equations are solved. In the second step the continuity and the momentum equation with only the pressure terms are solved to give a block tri‐diagonal system of equation with pressure as the unknown. This system can be solved by a direct matrix solver without iteration. A new implicit treatment of non‐hydrostatic pressure, similar to the lower layers is applied to the top layer which makes the model free of any hydrostatic pressure assumption all through the water column. This treatment enables the model to evaluate both free surface elevation and wave celerity more accurately. A series of numerical tests including free‐surface flows with significant vertical accelerations and nonlinear behaviour in shoaling zone are performed. Comparison between numerical results, analytical solutions and experimental data demonstrates a satisfactory performance. Copyright © 2006 John Wiley & Sons, Ltd.     

A fully hydrodynamic model for three‐dimensional,free‐surface flows

The hydrostatic pressure assumption has been widely used in studying water movements in rivers, lakes, estuaries, and oceans. While this assumption is valid in many cases and has been successfully used in numerous studies, there are many cases where this assumption is questionable. This paper presents a three‐dimensional, hydrodynamic model for free‐surface flows without using the hydrostatic pressure assumption. The model includes two predictor–corrector steps. In the first predictor–corrector step, the model uses hydrostatic pressure at the previous time step as an initial estimate of the total pressure field at the new time step. Based on the estimated pressure field, an intermediate velocity field is calculated, which is then corrected by adding the non‐hydrostatic component of the pressure to the estimated pressure field. A Poisson equation for non‐hydrostatic pressure is solved before the second intermediate velocity field is calculated. The final velocity field is found after the free surface at the new time step is computed by solving a free‐surface correction equation. The numerical method was validated with several analytical solutions and laboratory experiments. Model results agree reasonably well with analytical solutions and laboratory results. Model simulations suggest that the numerical method presented is suitable for fully hydrodynamic simulations of three‐dimensional, free‐surface flows. Copyright © 2003 John Wiley & Sons, Ltd.     

Semi‐implicit subgrid modelling of three‐dimensional free‐surface flows

In this paper, a semi‐implicit numerical model for two‐ and three‐dimensional free‐surface flows will be formulated in such a fashion as to intrinsically account for subgrid bathymetric details. It will be shown that with the proposed subgrid approach the model accuracy can be substantially improved without increasing the corresponding computational effort. Copyright © 2010 John Wiley & Sons, Ltd.     

Vortex particle‐in‐cell method for three‐dimensional viscous unbounded flow computations

A new vortex particle‐in‐cell (PIC) method is developed for the computation of three‐dimensional unsteady, incompressible viscous flow in an unbounded domain. The method combines the advantages of the Lagrangian particle methods for convection and the use of an Eulerian grid to compute the diffusion and vortex stretching. The velocity boundary conditions used in the method are of Dirichlet‐type, and can be calculated using the vorticity field on the grid by the Biot–Savart equation. The present results for the propagation speed of the single vortex ring are in good agreement with the Saffman's model. The applications of the method to the head‐on and head‐off collisions of the two vortex rings show good agreement with the experimental and numerical literature. Copyright © 2000 John Wiley & Sons, Ltd.     

An implicit three‐dimensional fully non‐hydrostatic model for free‐surface flows

An implicit method is developed for solving the complete three‐dimensional (3D) Navier–Stokes equations. The algorithm is based upon a staggered finite difference Crank‐Nicholson scheme on a Cartesian grid. A new top‐layer pressure treatment and a partial cell bottom treatment are introduced so that the 3D model is fully non‐hydrostatic and is free of any hydrostatic assumption. A domain decomposition method is used to segregate the resulting 3D matrix system into a series of two‐dimensional vertical plane problems, for each of which a block tri‐diagonal system can be directly solved for the unknown horizontal velocity. Numerical tests including linear standing waves, nonlinear sloshing motions, and progressive wave interactions with uneven bottoms are performed. It is found that the model is capable to simulate accurately a range of free‐surface flow problems using a very small number of vertical layers (e.g. two–four layers). The developed model is second‐order accuracy in time and space and is unconditionally stable; and it can be effectively used to model 3D surface wave motions. Copyright © 2004 John Wiley & Sons, Ltd.