Three-dimensional free surface modelling in an unstructured mesh environment for metal processing applications

In the casting of metals, tundish flow, welding, converters, and other metal processing applications, the behaviour of the fluid surface is important. In aluminium alloys, for example, oxides formed on the surface may be drawn into the body of the melt where they act as faults in the solidified product affecting cast quality. For this reason, accurate description of wave behaviour, air entrapment, and other effects need to be modelled, in the presence of heat transfer and possibly phase change. The authors have developed a single-phase algorithm for modelling this problem. The Scalar Equation Algorithm (SEA) (see Refs. 1 and 2), enables the transport of the property discontinuity representing the free surface through a fixed grid. An extension of this method to unstructured mesh codes is presented here, together with validation. The new method employs a TVD flux limiter in conjunction with a ray-tracing algorithm, to ensure a sharp bound interface. Applications of the method are in the filling and emptying of mould cavities, with heat transfer and phase change.

Free surface and mesh control for the FV simulation of semi-levitation melting

Abstract not available

Linear stern waves in finite depth channels

This paper formulates an analytically tractable problem for the wake generated by a long flat bottom ship by considering the steady free surface flow of an inviscid, incompressible fluid emerging from beneath a semi-infinite rigid plate. The flow is considered to be irrotational and two-dimensional so that classical potential flow methods can be exploited. In addition, it is assumed that the draft of the plate is small compared to the depth of the channel. The linearised problem is solved exactly using a Fourier transform and the Wiener-Hopf technique, and it is shown that there is a family of subcritical solutions characterised by a train of sinusoidal waves on the downstream free surface. The amplitude of these waves decreases as the Froude number increases. Supercritical solutions are also obtained, but, in general, these have infinite vertical velocities at the trailing edge of the plate. Consideration of further terms in the expansions suggests a way of canceling the singularity for certain values of the Froude number.

What is the apparent angle of a Kelvin ship wave pattern?

While the half-angle which encloses a Kelvin ship wave pattern is commonly accepted to be 19.47 degrees, recent observations and calculations for sufficiently fast-moving ships suggest that the apparent wake angle decreases with ship speed. One explanation for this decrease in angle relies on the assumption that a ship cannot generate wavelengths much greater than its hull length. An alternative interpretation is that the wave pattern that is observed in practice is defined by the location of the highest peaks; for wakes created by sufficiently fast-moving objects, these highest peaks no longer lie on the outermost divergent waves, resulting in a smaller apparent angle. In this paper, we focus on the problems of free surface flow past a single submerged point source and past a submerged source doublet. In the linear version of these problems, we measure the apparent wake angle formed by the highest peaks, and observe the following three regimes: a small Froude number pattern, in which the divergent waves are not visible; standard wave patterns for which the maximum peaks occur on the outermost divergent waves; and a third regime in which the highest peaks form a V-shape with an angle much less than the Kelvin angle. For nonlinear flows, we demonstrate that nonlinearity has the effect of increasing the apparent wake angle so that some highly nonlinear solutions have apparent wake angles that are greater than Kelvin's angle. For large Froude numbers, the effect on apparent wake angle can be more dramatic, with the possibility of strong nonlinearity shifting the wave pattern from the third regime to the second. We expect our nonlinear results will translate to other more complicated flow configurations, such as flow due to a steadily moving closed body such as a submarine.

Effects of mould filling on evolution of the solid-liquid interface during solidification

This work presents a numerical analysis of simultaneous mould filling and phase change for solidification in a two-dimensional rectangular cavity. The role of residual flow strength and temperature gradients within the solidifying domain, caused by the filling process, on the evolution of solidification interface are investigated. An implicit volume of fluid (VOF)-based algorithm has been employed for simulating the free surface flows during the filling process, while the model for solidification is based on a fixed-grid enthalpy-based control volume approach. Solidification modeling is coupled with VOF through User Defined Functions developed in the commercial computational fluid dynamics (CFD) code FLUENT 6.3.26. Comparison between results of the conventional analysis without filling effect and those of the present analysis shows that the residual flow resulting from the filling process significantly influences the progress of the solidification interface. A parametric study is also performed with variables such as cooling rate, filling velocity and filling configuration, in order to investigate the coupled effects of the buoyancy-driven flow and the residual flow on the solidification behavior.

Volume conservation issues in incompressible smoothed particle hydrodynamics

A divergence-free velocity field is usually sought in numerical simulations of incompressible fluids. We show that the particle methods that compute a divergence-free velocity field to achieve incompressibility suffer from a volume conservation issue when a finite time-step position update scheme is used. Further, we propose a deformation gradient based approach to arrive at a velocity field that reduces the volume conservation issues in free surface flows and maintains density uniformity in internal flows while retaining the simplicity of first order time updates. (C) 2015 Elsevier Inc. All rights reserved.

An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Turbulent flow structure in experimental laboratory wind-generated gravity waves

This paper is the third part of a report on systematic measurements and analyses of wind-generated water waves in a laboratory environment. The results of the measurements of the turbulent flow on the water side are presented here, the details of which include the turbulence structure, the correlation functions, and the length and velocity scales. It shows that the mean turbulent velocity profiles are logarithmic, and the flows are hydraulically rough. The friction velocity in the water boundary layer is an order of magnitude smaller than that in the wind boundary layer. The level of turbulence is enhanced immediately beneath the water surface due to micro-breaking, which reflects that the Reynolds shear stress is of the order u *w 2. The vertical velocities of the turbulence are related to the relevant velocity scale at the still-water level. The autocorrelation function in the vertical direction shows features of typical anisotropic turbulence comprising a large range of wavelengths. The ratio between the microscale and macroscale can be expressed as λ/Λ=a Re Λ n, with the exponent n slightly different from -1/2, which is the value when turbulence production and dissipation are in balance. On the basis of the wavelength and turbulent velocity, the free-surface flows in the present experiments fall into the wavy free-surface flow regime. The integral turbulent scale on the water side alone underestimates the degree of disturbance at the free surface. © 2012 Elsevier B.V.

Linear stability of miscible two-fluid flow down an incline

We show that miscible two-layer free-surface flows of varying viscosity down an inclined substrate are different in their stability characteristics from both immiscible two-layer flows, and flows with viscosity gradients spanning the entire flow. New instability modes arise when the critical layer of the viscosity transport equation overlaps the viscosity gradient. A lubricating configuration with a less viscous wall layer is identified to be the most stabilizing at moderate miscibility (moderate Peclet numbers). This also is in contrast with the immiscible case, where the lubrication configuration is always destabilizing. The co-existence that we find under certain circumstances, of several growing overlap modes, the usual surface mode, and a Tollmien-Schlichting mode, presents interesting new possibilities for nonlinear breakdown. © 2013 AIP Publishing LLC.

A SINGULAR EXACT SOLUTION FOR NONLINEAR SHALLOW-WATER ON A ROTATING EARTH

In this paper, we present an exact solution for nonlinear shallow water on a rotating planet. It is a kind of solitary waves with always negative wave height and a celerity smaller than linear shallow water propagation speed square-root gh. In fact, it propagates with a speed equal to (1 + a/h) square-root gh(1 + a/h) where a is the negative wave height. The lowest point of the water surface is a singular point where the first order derivative has a discontinuity of the first kind. The horizontal scale of the wave has actually no connection with the water depth.

The hydraulics and resulting bed scour within the vicinity of Submerged Single Span Arch Bridges

As a consequence of climate change there is now a more frequent occurrence of extreme rainfall events where, with higher rates of urbanisation, the built environment has become increasingly affected by flooding.. This is of particular importance in relation to the stability of bridge structures that span rivers and canals etc. In November 2009, the UK and Ireland were subjected to extraordinarily severe weather conditions for several days. The rainfall was logged as the highest level of rainfall ever recorded within the UK, and as a direct consequence, unprecedented flooding occurred in Cumbria. This flooding led to the collapse of three road bridges which were generally 19th century masonry arch bridges, with relatively shallow foundations. In the UK, knowledge of the combined effect of bridge scouring and inundation has been not been particularly widely studied. Research carried out by Hamill et al [1] considered the hydraulic analysis of single arch bridges under flood conditions, but no consideration was given towards the likely damage to these structures due to scouring. Prior to this, Bierry and Delleur [2] produced a classic paper in predicting the discharge downstream of an inundated arch, focussing on predicting afflux as opposed to bridge scour. Further work on backwater effects was carried out by Martin-Vide & Prio [3] in semi-circular arch bridges. Both pressurized and free-surface flows at the bridge were investigated. Flows on a mobile bed in clear-water conditions were compared to those with a rigid bed, but no predictive equation for scour under pressurised conditions was considered. This paper will present initial findings from an experimental investigation into the effects of surcharged flow and subsequent scour within the vicinity of single span arch bridges. Velocities profiles will be shown within the vicinity of the arch, in addition to the depth of clear water scour, for a series of flows and model spans. The data will be presented, where results will be correlated to the most recent predictive equations that are proposed.

Stability of Short-span Bridges Subject to Overtopping during Floods: Case Studies from UK and the USA

As a consequence of increased levels of flooding, largely attributable to urbanization of watersheds (and perhaps climate change, more frequent extreme rainfall events are occurring and threatening existing critical infrastructure. Many of which are short-span bridges over relatively small waterways (e.g., small rivers, streams and canals). Whilst these short-span bridges were designed, often many years ago, to pass relatively minor the then standard return-period floods, in recenttimes the failure incidence of such short-span bridges has been noticeably increasing. This is suggestive of insufficient hydraulic capacity or alternative failure mechanism not envisaged at the time of design e.g. foundation scour or undermining. This paper presen ts, and draws lessons, from bridge failures in Ireland and the USA. For example, in November 2009, the UK and Ireland were subjected to extraordinarily severe weather conditions for several days. The resulting flooding led to the collapse of three UK bridges that were generally 19th century masonry arch bridges, withrelatively shallow foundations. Parallel failure events have been observed in the USA. To date, knowledge of the combined effect of waterway erosion, bridge submergence, and geotechnical collapse has not been adequately studied. Recent research carried out considered the hydraulic analysis of short span bridges under flood conditions, but no consideration was given towards the likely damage to these structures due to erosive coupling of hydraulic and geotechnical factors. Some work has been done to predict the discharge downstream of an inundated arch, focusing onpredicting afflux, as opposed to bridge scour, under both pressurized and free-surface flows, but no ! predictive equation for scour under pressurized conditions was ever considered. The case studies this paper presents will be augmented by the initial findings from the laboratory experiments investigating the effects of surcharged flow and subsequent scour within the vicinity of single span arch bridges. Velocities profiles will be shown within the vicinity of the arch, in addition to the depth of consequent scour, for a series of flows and model spans. The data will be presented and correlated to the most recent predictive equations for submerged contraction and abutment scour. The accuracy of these equations is examined, and the findings used as a basis for developing further studies in relation to short span bridges.

A combined BIE-FE method for the Stokes equations

A numerical algorithm for the biharmonic equation in domains with piecewise smooth boundaries is presented. It is intended for problems describing the Stokes flow in the situations where one has corners or cusps formed by parts of the domain boundary and, due to the nature of the boundary conditions on these parts of the boundary, these regions have a global effect on the shape of the whole domain and hence have to be resolved with sufficient accuracy. The algorithm combines the boundary integral equation method for the main part of the flow domain and the finite-element method which is used to resolve the corner/cusp regions. Two parts of the solution are matched along a numerical ‘internal interface’ or, as a variant, two interfaces, and they are determined simultaneously by inverting a combined matrix in the course of iterations. The algorithm is illustrated by considering the flow configuration of ‘curtain coating’, a flow where a sheet of liquid impinges onto a moving solid substrate, which is particularly sensitive to what happens in the corner region formed, physically, by the free surface and the solid boundary. The ‘moving contact line problem’ is addressed in the framework of an earlier developed interface formation model which treats the dynamic contact angle as part of the solution, as opposed to it being a prescribed function of the contact line speed, as in the so-called ‘slip models’. Keywords: Dynamic contact angle; finite elements; free surface flows; hybrid numerical technique; Stokes equations.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.