A method to compute ship exhaust plumes with waves and wind

Thermal and concentration transport models are implemented in CFDShip‐Iowa version 4.5, a semi‐coupled solver for air/water free surface flow (Int. J. Numer. Meth. Fluids 2008; 58 (6):591–624), to investigate the exhaust plume around ship superstructures. An incompressible, variable density approximation is implemented where the density can change in all governing equations due to temperature variations only. The thermal and concentration models are tested for the cases of steady and unsteady flow with thermal and solution transport in a 2D square cavity, and for a 3D thermal plume in an open environment, showing good agreement between computational results and experimental data. To test the method in an extreme motions condition, the exhaust plume of the ONR Tumblehome model DTMB 5613 is studied, showing complicated vortical structures in air including a pair of counter‐rotating vortices downstream of the stack for cross‐flow, and bended bird‐plume shape in the symmetry plane and varying arc‐shape in axial sections both for temperature and NO_x concentration fields. Effects of smoke exhaust speed and wind speed on the temperature and concentration distributions are studied. Finally, a smoke downwash computation is performed for a ship free to move in 6 degrees of freedom in a sea state 8 condition. Copyright © 2010 John Wiley & Sons, Ltd.     

Multiphase modelling of pouring glass over the spout lip of an industrial float in the flat glass forming process

This paper describes the numerical modelling of the three‐phase flow that is established when pouring molten glass onto a tin bath (float) under a non‐reactive atmosphere. Discharged from the furnace, the glass is spouted over a lip onto the bath, where it floats and spreads out in the form of a long ribbon, leading to the flat glass forming process. Numerically, the physical modelling must consider a multiphase problem with three immiscible phases: glass, tin and air. The simulation of this complex free‐surface flow involves the solution of the Navier–Stokes set of equations for all the phases simultaneously, using a volume‐of‐fluid formulation that introduces a marker function convected by the flow to identify each phase. The evolution of the interphases is tracked over time with the implementation of a continuous surface force algorithm. A general purpose, well‐tested commercial code, FLUENT, is employed for the computations. Firstly, a two‐dimensional model considering the symmetry plane of the float is developed to fit accurate physical and numerical parameters. The high complexity of the interphases in addition to great differences between the physical properties of the phases has required extensive tests to ensure the consistency and accuracy of the solutions. Afterwards, a complete three‐dimensional model is built to simulate the pouring process in a real geometry and predict the behaviour of the industrial facility when the operating conditions are modified. Typical flow phenomena inside the tin bath, such as tin currents, wet back flow or even the equilibrium thickness, are also obtained as a promising result of the numerical modelling. Copyright © 2008 John Wiley & Sons, Ltd.     

A lattice Boltzmann model for diffusion of binary gas mixtures that includes diffusion slip

This paper describes the development of a lattice Boltzmann (LB) model for a binary gas mixture, and applications to channel flow driven by a density gradient with diffusion slip occurring at the wall. LB methods for single component gases typically use a non‐physical equation of state in which the relationship between pressure and density varies according to the scaling used. This is fundamentally unsuitable for extension to multi‐component systems containing gases of differing molecular masses. Substantial variations in the species densities and pressures may exist even at low Mach numbers; hence, the usual linearized equation of state for small fluctuations is unsuitable. Also, existing methods for implementing boundary conditions do not extend easily to novel boundary conditions, such as diffusion slip. The new model developed for multi‐component gases avoids the pitfalls of some other LB models. A single computational grid is shared by all the species, and the diffusivity is independent of the viscosity. The Navier–Stokes equation for the mixture and the Stefan–Maxwell diffusion equation are both recovered by the model. Diffusion slip, the non‐zero velocity of a gas mixture at a wall parallel to a concentration gradient, is successfully modelled and validated against a simple one‐dimensional model for channel flow. To increase the accuracy of the scheme, a second‐order numerical implementation is needed. This may be achieved using a variable transformation method that does not increase the computational time. Simulations were carried out on hydrogen and water diffusion through a narrow channel for varying total pressure and concentration gradients. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of mesh deformation on the quality of large eddy simulations

The influence of mesh motion on the quality of large eddy simulation (LES) was studied in the present article. A three‐dimensional, turbulent pipe flow (Re_τ=360) was considered as a test case. Simulations with both stretching and static meshes were carried out in order to understand how mesh motion affects the turbulence statistics. The spatial filtering of static and moving mesh direct numerical simulation (DNS) data showed how an ideal LES would perform, while the comparison of DNS cases with static and moving meshes revealed that no significant numerical errors arise from the mesh motion when the simulation is fully resolved. The comparison of the filtered fields of the DNS with a moving mesh with the corresponding LES fields revealed different responses to mesh motion from different numerical approaches. A straightforward test was applied in order to verify that the moving mesh works consistently in LES: when the mesh is stretched in the streamwise direction, the moving mesh results should be in between the two extremal resolutions between which the mesh is stretched. Numerical investigations using four different LES approaches were carried out. In addition to the Smagorinsky model, three implicit LES approaches were used: linear interpolation (non‐dissipative), the Gamma limiter (dissipative), and the scale‐selective discretisation (slightly dissipative). The results indicate that while the Smagorinsky and the scale‐selective discretisation approaches produce results consistent with the resolution of the non‐static mesh, the implicit LES with linear interpolation or the Gamma scheme do not. Copyright © 2016 John Wiley & Sons, Ltd.     

Modeling Transverse Dispersion and Variable Density Flow in Porous Media

A two-dimensional numerical model is used to study the nonlinear behavior of density gradients on transverse dispersion. Numerical simulations are conducted using d ³ f, a computer code for simulation of density-dependent flow in porous media. Considering a density-stratified horizontal flow in a heterogeneous porous media, a series of simulations is carried out to examine the effect of the density gradient on macro-scale transverse dispersivity. Changing salt concentration significantly affects fluid properties. This physical behavior of the fluid involves a non-linearity in modeling the interaction between salt and fresh water. It is concluded that the large-scale transport properties for high density flow deviate significantly from the tracer case due to the spatial variation of permeability, described by statistical parameters, at the local-scale. Indeed, the presence of vertical flow velocities induced by permeability variations is responsible for the reduction of the mixing zone width in the steady state in the case of a high density gradient. Uncertainties in the model simulations are studied in terms of discretization errors, boundary conditions, and convergence of ensemble averaging. With respect to the results, the gravity number appears to be the controlling parameter for dispersive flux. In addition, the applicability and limitations of the nonlinear model of Hassanizadeh (1990) and Hassanizadeh and Leijnse (1995) (Adv Water Resour 18(4):203–215, 1995) in heterogeneous porous media are investigated. We found that the main cause of the nonlinear behavior of dispersion, which is the interaction between density contrast and vertical velocity, needs to be explicitly accounted for in a macro-scale model.     

Global flow instability in a lid‐driven cavity

The stability of flow in a lid‐driven cavity is investigated using an accurate numerical technique based on a hybrid scheme with spectral collocation and high‐order finite differences. A global stability analysis is carried out and critical parameters are identified for various aspect ratios. It is found that while there is reasonable agreement with the literature for the critical parameters leading to loss of stability for the square cavity, there are significant discrepancies for cavities of aspect ratios 1.5 and 2. Simulations of the linearized unsteady equations confirm the results from the global stability analysis for aspect ratios A = 1, 1.5 and A = 2. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical calculations of erosion in an abrupt pipe contraction of different contraction ratios

Erosion predictions in a pipe with abrupt contraction of different contraction ratios for the special case of two‐phase (liquid and solid) turbulent flow with low particle concentration are presented. A mathematical model based on the time‐averaged governing equations of 2‐D axi‐symmetric turbulent flow is used for the calculations of the fluid velocity field (continuous phase). The particle‐tracking model of the solid particles is based on the solution of the governing equation of each particle motion taking into consideration the effect of particle rebound behaviour. Models of erosion were used to predict the erosion rate in mg/g. The effect of Reynolds number and flow direction with respect to the gravity was investigated for three contraction geometries considering water flow in a carbon steel pipe. The results show that the influence of the contraction ratio on local erosion is very significant. However, this influence becomes insignificant when the average erosion rates over the sudden contraction area are considered. The results also indicate the significant influence of inlet velocity variations. The influence of buoyancy is significant for the cases of low velocity of the continuous flow. A threshold velocity below which erosion may be neglected was indicated. Copyright © 2004 John Wiley & Sons, Ltd.     

Heat and mass transfer analogy for condensation of humid air in a vertical channel

This study examines energy transport associated with liquid film condensation in natural convection flows driven by differences in density due to temperature and concentration gradients. The condensation problem is based on the thin-film assumptions. The most common compositional gradient, which is encountered in humid air at ambient temperature is considered. A steady laminar Boussinesq flow of an ideal gas–vapor mixture is studied for the case of a vertical parallel plate channel. New correlations for the latent and sensible Nusselt numbers are established, and the heat and mass transfer analogy between the sensible Nusselt number and Sherwood number is demonstrated. Received on 15 November 1999     

A comparative study of high‐order variable‐property segregated algorithms for unsteady low Mach number flows

In this paper, five different algorithms are presented for the simulation of low Mach flows with large temperature variations, based on second‐order central‐difference or fourth‐order compact spatial discretization and a pressure projection‐type method. A semi‐implicit three‐step Runge–Kutta/Crank–Nicolson or second‐order iterative scheme is used for time integration. The different algorithms solve the coupled set of governing scalar equations in a decoupled segregate manner. In the first algorithm, a temperature equation is solved and density is calculated from the equation of state, while the second algorithm advances the density using the differential form of the equation of state. The third algorithm solves the continuity equation and the fourth algorithm solves both the continuity and enthalpy equation in conservative form. An iterative decoupled algorithm is also proposed, which allows the computation of the fully coupled solution. All five algorithms solve the momentum equation in conservative form and use a constant‐ or variable‐coefficient Poisson equation for the pressure. The efficiency of the fourth‐order compact scheme and the performances of the decoupling algorithms are demonstrated in three flow problems with large temperature variations: non‐Boussinesq natural convection, channel flow instability, flame–vortex interaction. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical investigation of fully developed laminar flow in irregular annuli: Square‐circular ducts

Numerical study is carried out on fully developed laminar axial flow of incompressible Newtonian fluid in irregular annuli of different cross‐sectional geometry. Accurate prediction of the friction factor of this flow was obtained using high‐order finite element method. Investigations are performed for the annular flow between rectangular and circular ducts under all possible arrangements. The governing equations are solved by an efficient finite element technique. It was found that a higher annular area‐ratio will lead to a monotonic increase in (f Re) value in the case of regular annuli, and will lead to an increase followed by a decrease in (f Re) value in the case of irregular annuli. Furthermore, this study illustrates that irregular annuli have lower (f Re) value than regular annuli, and that the square‐in‐circle case (SC‐case) has the lowest (f Re) value whereas the circle‐in‐circle case (CC‐case) has the highest (f Re) value. The finite element technique developed in this study was validated against the already reported numerical results in the literature, and the results presented in this work show excellent agreement. Copyright © 2009 John Wiley & Sons, Ltd.     

Investigation of the specific mass flow rate distribution in pipes supplied with a pulsating flow

A pulsating flow is typical of inlet and exhaust pipes of internal combustion engines and piston compressors. Unsteady flow phenomena are especially important in the case of turbocharged engines, because dynamic effects occurring in the exhaust pipe can affect turbine operation conditions and performance.One of the basic parameters describing the unsteady flow is a transient mass flow rate related to the instantaneous flow velocity, which is usually measured by means of hot-wire anemometers. For the flowing gas, it is more appropriate to analyze the specific mass flow rate φ_m = ρv, which takes into account also variations in the gas density. In order to minimize the volume occupied by measuring devices in the control section, special double-wire sensors for the specific mass flow rate (CTA) and temperature (CCT) measurement were applied. The article describes procedures of their calibration and measurement. Different forms of calibration curves are analyzed as well in order to match the approximation function to calibration points. Special attention is paid to dynamic phenomena related to the resonance occurring in a pipe for characteristic frequencies depending on the pipe length. One of these phenomena is a reverse flow, which makes it difficult to interpret properly the recorded CTA signal. Procedures of signal correction are described in detail. To verify the measurements, a flow field investigation was carried out by displacing probes radially and determining the profiles of the specific mass flow rate under the conditions of a steady and pulsating flow. The presence and general features of a reverse flow, which was identified experimentally, were confirmed by 1-D unsteady flow calculations.     

Three‐dimensional finite element modelling of coupled free/porous flows: applications to industrial and environmental flows

Conjunctive modelling of free/porous flows provides a powerful and cost‐effective tool for designing industrial filters used in the process industry and also for quantifying surface–subsurface flow interactions, which play a significant role in urban flooding mechanisms resulting from sea‐level rise and climate changes. A number of well‐established schemes are available in the literature for simulation of such regimes; however, three‐dimensional (3D) modelling of such flow systems still presents numerical and practical challenges. This paper presents the development of a fully 3D, transient finite element model for the prediction and quantitative analyses of the hydrodynamic behaviour encountered in industrial filtrations and environmental flows represented by coupled flows. The weak‐variational formulation in this model is based on the use of C⁰ continuous equal‐order Lagrange polynomial functions for velocity and pressure fields represented by 3D hexahedral finite elements. A mixed UVWP finite element scheme based on the standard Galerkin technique satisfying the Ladyzhenskaya–Babuska–Brezzi stability criterion through incorporation of an artificial compressibility term in the continuity equation has been employed for the solution of coupled partial differential equations. We prove that the discretization generates unified stabilization for both the Navier–Stokes and Darcy equations and preserves the geometrical flexibility of the computational grids. A direct node‐linking procedure involving the rearrangement of the global stiffness matrix for the interface elements has been developed by the authors, which is utilized to couple the governing equations in a single model. A variety of numerical tests are conducted, indicating that the model is capable of yielding theoretically expected and accurate results for free, porous and coupled free/porous problems encountered in industrial and environmental engineering problems representing complex filtration (dead‐end and cross‐flow) and interacting surface–subsurface flows. The model is computationally cost‐effective, robust, reliable and easily implementable for practical design of filtration equipments, investigation of land use for water resource availability and assessment of the impacts of climatic variations on environmental catastrophes (i.e. coastal and urban floods). The model developed in this work results from the extension of a multi‐disciplinary project (AEROFIL) primarily sponsored by the European aerospace industries for development of a computer simulation package (Aircraft Cartridge Filter Analysis Modelling Program), which was successfully utilized and deployed for designing hydraulic dead‐end filters used in Airbus A380.Copyright © 2012 John Wiley & Sons, Ltd.     

Modified Lagrangian vortex method with improved boundary conditions for water waves past a thin bottom‐standing barrier

Herein, the modified Lagrangian vortex method (LVM), a hybrid analytical‐numerical algorithm per se, is devised to simulate the process of vortex formation and shedding from the sharp edge of a zero‐thickness vertical plate under linear water‐wave attack. Application of the Helmholtz decomposition facilitates a convenient switch between the inviscid‐ and viscous‐flow models, thereby enabling easy incorporation of vorticity effects into the potential‐flow calculations for the viscous‐dominated region. In evaluating the potential‐flow component, making good use of the quickly convergent technique with singular basis functions, correctly capturing the singular behavior in velocity fields near the tip of the plate, leads to a considerable reduction of computational burdens and to 12‐decimal‐place accuracy. The viscous correction is carried out via the meshless LVM with improved boundary conditions. Comparisons with previously published results show good agreement. Simulations of vortex generation and evolution illuminate the ability of the present method, and provide a supplement to pertinent experimental works. The hybrid scheme proposed herein allows flexibility for the former LVM and convenience in the code development. Such a compromise fits particularly well for the high‐resolution modeling of sharp‐edged vortex shedding without heavy numerical developments. Copyright © 2014 John Wiley & Sons, Ltd.     

Simultaneous effects of induced magnetic field and heat and mass transfer on the peristaltic motion of second‐order fluid in a channel

In this article, we investigate the influence of heat and mass transfer on the peristaltic flow of magnetohydrodynamic second‐order fluid in a channel when the induced magnetic field effects are present. Problem formulation in a wave frame of reference is presented. The governing nonlinear analysis is carried out under the assumption of small wave number. Explicit expressions of the pressure gradient, the stream function, the magnetic force function, the axial induced magnetic field, the current density distribution, the temperature, and the concentration distribution are derived. The effects of embedded parameters are also examined. Copyright © 2011 John Wiley & Sons, Ltd.     

Thermocapillary migration of liquid droplets in a temperature gradient in a density matched system

An experimental investigation of thermocapillary flow in droplets of a vegetable oil (partially hydrogenated soybean oil) immersed in silicone oil was conducted in a test cell with a heated top wall and a cooled bottom wall. The liquids are nearly immiscible and have equal densities at a temperature below the room temperature, thus providing a simulation of low-gravity conditions by reducing the buoyancy forces. The interfacial tension between the two oils was measured in the temperature range 20–50°C using a capillary tube and d/dT was determined to be negative. Droplets ranging in sizes from 3 mm to 1 cm diameter were injected into the silicone oil. The vertical temperature profile in the bulk liquid (silicone oil) produces temperature variations along the interface which induce variations in the interfacial tension. The flow inside the droplet driven by the resulting interfacial shear stresses was observed using a laser light-sheet flow visualization technique. The flow direction is consistent with the sign of d/dT. The observed maximum surface velocities are compared to the theoretical predictions of Young et al. (1959).For short times after injection, the droplets were driven by this flow towards the hot wall above the matched-density temperature until the droplets reached a point where the forces due to the flow and buoyancy were equal. After longer times, the droplets moved to the cold side due to suspected density changes caused by mass transfer from the droplets to the silicone oil. This was confirmed by tests under isothermal conditions, where it was observed that droplets of all sizes fell to the cold bottom eventually.Thus, even though the thermocapillary flow inside the droplets persisted for long times in spite of the mass transfer, the migration of droplets towards the hot side was eventually affected by uncontrolled buoyancy forces resulting from density changes due to mass transfer. While additional liquids are being tried, it is suggested from the present experience that reduced gravity experiments will probably be necessary to provide unambiguous data for the migration of droplets.     

Similarity solutions for the flow behind an exponential shock in a non-ideal gas

Similarity solutions for the flow of a non-ideal gas behind a strong exponential shock driven out by a piston (cylindrical or spherical) moving with time according to an exponential law are obtained. Similarity solutions exist only when the surrounding medium is of constant density. Solutions are obtained, in both the cases, when the flow between the shock and the piston is isothermal or adiabatic. It is found that the assumption of zero temperature gradient brings a profound change in the density distribution as compare to that of the adiabatic case. Effects of the non-idealness of the gas on the flow-field between the shock and the piston are investigated. The variations of density-ratio across the shock and the location of the piston with the parameter of non-idealness of the gas are also obtained.     

Numerical simulation of heat,mass and momentum transport behaviours in directionally solidifying alloy castings under electromagnetic fields using an extended Direct‐SIMPLE scheme

A Direct‐SIMPLE scheme is further extended for numerical simulations of heat and species mass transfer, and liquid flow behaviours in solidification processes of shaped alloy castings under various electromagnetic (EM) fields, based on a binary continuum model (a solidification transport phenomena/processes (STP)‐based dendrite solidification model). Unlike in a SIMPLE scheme, no iterative computations are needed to achieve the final pressure and velocity corrections in the extended Direct‐SIMPLE procedure, therefore extremely high computational efforts can be avoided. Using three different types of model alloys, pseudo‐binary In718 base‐4.85 wt% Nb, γ(TiAl)‐55 at% Al and Al‐4.5 wt% Cu systems, sample computations for solving strongly coupled solidification transport phenomena in directionally solidifying shaped castings under static and harmonic EM‐fields of different strengths are carried out, to demonstrate the feasibility and efficient calculation performance of the present model and numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

The effects of control‐volume distributed multi‐point flux approximation (CVD‐MPFA) on upscaling—A study using the CVD‐MPFA schemes

A family of flux‐continuous, locally conservative, finite‐volume schemes has been developed for solving the general geometry‐permeability tensor (petroleum reservoir‐simulation) pressure equation on structured and unstructured grids and are control‐volume distributed (textit Comput. Geo. 1998; 2 :259–290; Comput. Geo. 2002; 6 :433–452). The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir‐simulation schemes (two‐point flux approximation) when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization (Int. J. Numer. Meth. Fluids 2006; 51 :1177–1203). Improved convergence (for two‐ and three‐dimensional formulation) using the quadrature parameterization has been observed for the family of flux‐continuous control‐volume distributed multi‐point flux approximation (CVD‐MPFA) schemes (Ph.D. Thesis, University of Wales, Swansea, U.K., 2007). In this paper family of flux‐continuous (CVD‐MPFA) schemes are used as a part of numerical upscaling procedure for upscaling the fine‐scale grid information (permeability) onto a coarse grid scale. A series of data‐sets (SPE, 2001) are tested where the upscaled permeability tensor is computed on a sequence of grid levels using the same fixed range of quadrature points in each case. The refinement studies presented involve:

• (i) Refinement comparison study: In this study, permeability distribution for cells at each grid level is obtained by upscaling directly from the fine‐scale permeability field as in standard simulation practice.
• (ii) Refinement study with renormalized permeability: In this refinement comparison, the local permeability is upscaled to the next grid level hierarchically, so that permeability values are renormalized to each coarser level. Hence, showing only the effect of increased grid resolution on upscaled permeability, compared with that obtained directly from the fine‐scale solution.
• (iii) Refinement study with invariant permeability distribution: In this study, a classical mathematical convergence test is performed. The same coarse‐scale underlying permeability map is preserved on all grid levels including the fine‐scale reference solution.

The study is carried out for the discretization of the scheme in physical space. The benefit of using specific quadrature points is demonstrated for upscaling in this study and superconvergence is observed. Copyright © 2010 John Wiley & Sons, Ltd.     

Predictions of the excess pressure drop of Boger fluids through a 2:1:2 contraction–expansion geometry using non-equilibrium molecular dynamics

Non-equilibrium molecular dynamics are used to generate the flow of polymer solutions, specifically of Boger fluids, through a planar 2:1:2 contraction–expansion geometry. The solvent molecules are represented by Lennard–Jones particles, while linear molecules are described by spring-monomers with a finite extensible non-linear elastic spring potential. The equations for Poiseuille flow are solved using a multiple time-scale algorithm extended to non-equilibrium situations. Simulations are performed at constant temperature using Nose–Hoover dynamics. At simulation conditions, changes in concentration show no significant effect on molecular conformation, velocity profiles, and stress fields, while variations in the Deborah number have a strong influence on fluid response. Increasing the magnitude of the Deborah number (De), larger deformation rates are developed in the flow region. For a Deborah number of one, the non-dimensional pressure drop presents values lower than the correspondent Newtonian case. However, for large Deborah numbers, the pressure drop increases above the Newtonian reference. An effective excess pressure drop above the Newtonian value is predicted for Boger fluids along this geometry.     

Effects of Knudsen number and geometry on gaseous flow and heat transfer in a constricted microchannel

A flow and heat transfer numerical simulation is performed for a 2D laminar incompressible gas flow through a constricted microchannel in the slip regime with constant wall temperature. The effects of rarefaction, creeping flow, first order slip boundary conditions and hydrodynamically/thermally developing flow are assumed. The effects of Knudsen number and geometry on thermal and hydrodynamic characteristics of flow in a constricted microchannel are explored. SIMPLE algorithm in curvilinear coordinate is used to solve the governing equations including continuity, energy and momentum with the temperature jump and velocity slip conditions at the solid walls in discretized form. The resulting velocity and temperature profiles are then utilized to obtain the microchannel C _f Re and Nusselt number as a function of Knudsen number and geometry. The results show that Knudsen number has declining effect on the C _f Re and Nusselt number in the constricted microchannel. In addition, the temperature jump on wall and slip velocity increase with increasing Knudsen number. Moreover, by decreasing the throttle area, the fluid flow characteristics experience more intense variations in the constricted region. To verify the code a comparison is carried out with available results and good agreement is achieved.