Revisiting the carrier-phase source-term formulation in mixed Lagrangian-Eulerian particle models

The formulation of the carrier-phase momentum and enthalpy source terms in mixed Lagrangian-Eulerian models of particle-laden flows is frequently reported inaccurately. Under certain circumstances, this can lead to erroneous implementations, which violate physical laws. A particle- rather than carrier-based approach is suggested for a consistent treatment of these terms.

Viscoelastic flow through a planar contraction using a semi-Lagrangian finite volume method

A semi-Lagrangian finite volume scheme for solving viscoelastic flow problems is presented. A staggered grid arrangement is used in which the dependent variables are located at different mesh points in the computational domain. The convection terms in the momentum and constitutive equations are treated using a semi-Lagrangian approach in which particles on a regular grid are traced backwards over a single time-step. The method is applied to the 4 : 1 planar contraction problem for an Oldroyd B fluid for both creeping and inertial flow conditions. The development of vortex behaviour with increasing values of We is analyzed.

Conservative semi-Lagrangian finite volume schemes

Semi-Lagrangian finite volume schemes for the numerical approximation of linear advection equations are presented. These schemes are constructed so that the conservation properties are preserved by the numerical approximation. This is achieved using an interpolation procedure based on area-weighting. Numerical results are presented illustrating some of the features of these schemes.

A semi-Lagrangian finite volume method for Newtonian contraction flows

A new finite volume method for solving the incompressible Navier--Stokes equations is presented. The main features of this method are the location of the velocity components and pressure on different staggered grids and a semi-Lagrangian method for the treatment of convection. An interpolation procedure based on area-weighting is used for the convection part of the computation. The method is applied to flow through a constricted channel, and results are obtained for Reynolds numbers, based on half the flow rate, up to 1000. The behavior of the vortex in the salient corner is investigated qualitatively and quantitatively, and excellent agreement is found with the numerical results of Dennis and Smith [Proc. Roy. Soc. London A, 372 (1980), pp. 393-414] and the asymptotic theory of Smith [J. Fluid Mech., 90 (1979), pp. 725-754].

A mixed Eulerian-Lagrangian approach for metal extrusion

In this paper a mixed Eulerian-Lagrangian approach for the modelling metal extrusion processes is presented. The approach involves the solution of non-Newtonian fluid flow equations in an Eulerian context, using a free-surface algorithm to track the behaviour of the workpiece and its extrusion. The solid mechanics equations associated with the tools are solved in Lagangrian context. Thermal interactions between the workpiece are modelled and a fluid-structure interaction technique is employed to model the effect of the fluid traction load imposed by the workpiece on the tools. Two extrusion test cases are investigated and the results obtained show the potential of the model with regard to representing the physics of the process and the simulation time.

An augmented Lagrangian contact algorithm employing a vertex-based finite volume method

Abstract not available

A semi-Lagrangian finite volume method for solving viscoelastic flow problems

Recently, there has been considerable interest in solving viscoelastic problems in 3D particularly with the improvement in modern computing power. In many applications the emphasis has been on economical algorithms which can cope with the extra complexity that the third dimension brings. Storage and computer time are of the essence. The advantage of the finite volume formulation is that a large amount of memory space is not required. Iterative methods rather than direct methods can be used to solve the resulting linear systems efficiently.

On the modelling of the thermal interactions between a spray curtain and an impinging cold gas cloud

A mixed Lagrangian-Eulerian model of a Water Curtain barrier is presented. The heat, mass and momentum processes are modelled in a Lagrangian framework for the dispersed phase and in an Eulerian framework for the carrier phase. The derivation of the coupling source terms is illustrated with reference to a given carrier phase cell. The turbulent character of the flow is treated with a single equation model, modified to directly account for the influence of the particles on the flow. The model is implemented in the form of a 2 D incompressible Navier Stokes solver, coupled to an adaptive Rung Kutta method for the Lagrangian sub-system. Simulations of a free standing full cone water spray show satisfactory agreement with experiment. Predictions of a Water Curtain barrier impacted by a cold gas cloud point to markedly different flow fields for the upward and downward configurations, which could influence the effectiveness of chemical absorption in the liquid phase.

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.

Multiphysics modeling and its application to the casting process

A modeling strategy is presented to solve the governing equations of fluid flow, temperature (with solidification), and stress in an integrated manner. These equations are discretized using finite volume methods on unstructured grids, which provide the capability to represent complex domains. Both the cell-centered and vertex-based forms of the finite volume discretization procedure are explained, and the overall integrated solution procedure using these techniques with suitable solvers is detailed. Two industrial processes, based on the casting of metals, are used to demonstrate the capabilities of the resultant modeling framework. This manufacturing process requires a high degree of coupling between the governing physical equations to accurately predict potential defects. Comparisons between model predictions and experimental observations are given.

A computational model for defect prediction in shape castings based on the interaction of free surface flow, heat transfer, and solidification phenomena

High-integrity castings require sophisticated design and manufacturing procedures to ensure they are essentially macrodefect free. Unfortunately, an important class of such defects—macroporosity, misruns, and pipe shrinkage—are all functions of the interactions of free surface flow, heat transfer, and solidication in complex geometries. Because these defects arise as an interaction of the preceding continuum phenomena, genuinely predictive models of these defects must represent these interactions explicitly. This work describes an attempt to model the formation of macrodefects explicitly as a function of the interacting continuum phenomena in arbitrarily complex three-dimensional geometries. The computational approach exploits a compatible set of finite volume procedures extended to unstructured meshes. The implementation of the model is described together with its testing and a measure of validation. The model demonstrates the potential to predict reliably shrinkage macroporosity, misruns, and pipe shrinkage directly as a result of interactions among free-surface fluid flow, heat transfer, and solidification.

Model sensitivity and uncertainty analysis using roadside air quality measurements

Most of the air quality modelling work has been so far oriented towards deterministic simulations of ambient pollutant concentrations. This traditional approach, which is based on the use of one selected model and one data set of discrete input values, does not reflect the uncertainties due to errors in model formulation and input data. Given the complexities of urban environments and the inherent limitations of mathematical modelling, it is unlikely that a single model based on routinely available meteorological and emission data will give satisfactory short-term predictions. In this study, different methods involving the use of more than one dispersion model, in association with different emission simulation methodologies and meteorological data sets, were explored for predicting best CO and benzene estimates, and related confidence bounds. The different approaches were tested using experimental data obtained during intensive monitoring campaigns in busy street canyons in Paris, France. Three relative simple dispersion models (STREET, OSPM and AEOLIUS) that are likely to be used for regulatory purposes were selected for this application. A sensitivity analysis was conducted in order to identify internal model parameters that might significantly affect results. Finally, a probabilistic methodology for assessing urban air quality was proposed.

The application of free-surface techniques to metal forming processes

In this paper the use of free-surface techniques, within the framework of a finite volume methodology, are investigated for the simulation of metal forming processes. In such processes, for example extrusion and forging, a workpiece is subjected to large scale deformation to create the product's shape. The use of Eulerian free-surface techniques to predict this final shape offers the advantage, over the traditionally used Lagrangian finite element method, of not requiring remmeshing. Two free-surface techniques to predict this final shape offers the advantage, over the traditionally used Lagrangian finite element method, of not requiring remesingh. Two free-surface techniques are compared by modelling a typical example of this type of process - non-Newtonian extrusion of an aluminium workpiece through a conical die.