Search results

The purpose of this paper is to show how simulation of the flow of particulates and fluids using discrete element modelling (DEM) and smoothed particle dynamics (SPH) particle methods, offer opportunities for better understanding the dynamics of flow processes.

DEM and SPH methods are demonstrated in a broad range of computationally‐demanding applications including comminution, biomedical, geophysical extreme flow events (risk/disaster modelling), eating of food by humans and elite water‐based sports.

DEM is ideally suited to predicting industrial and geophysical applications where collisions between particles are the dominant physics. SPH is highly suited to multi‐physics fluid flow applications in industrial, biophysical and geophysical applications. The advantages and disadvantages of these particle methods are discussed.

Research results are limited by the numerical resolution that can currently be afforded.

The paper demonstrates the use of particle‐based computational methods in a series of high value applications. Enterprises that share interests in these applications will benefit in their product and service development by adopting these methods.

The ability to model disasters provides governments and companies with the opportunity and obligation to use these to render knowable disasters which were previously considered unknowable. The ability to predict the breakdown of food during eating opens up opportunities for the design of superior performing foods with lower salt, sugar and fat that can directly contribute to improved health outcomes and can influence government food regulatory policy.

The paper extends the scale and range of modelling of particle methods for demanding leading‐edge problems, of practical interest in engineering and applied sciences.

Presents three non‐isothermal, time dependent, three dimensional examples having cylindrical geometries to show the significant effort of numerical precision and dissipation on rotating flow predictions. The examples are relevant to turbomachinery design and geophysical studies. Discusses the relationship between numerical precision, numerical dissipation and co‐ordinate system angular velocity. Compares predictions made in stationary and rotating co‐ordinate systems, using contour plots of dimensionless stream function and temperature. Shows that wrong, axisymmetric solutions are predicted if the co‐ordinate system is not selected to minimize relative tangential velocities/Peclet numbers, thereby increasing numerical precision and reducing dissipation.

Particle scale simulation of industrial particle flows using discrete element method (DEM) offers the opportunity for better understanding the flow dynamics leading to improvements in equipment design and operation that can potentially lead to large increases in equipment and process efficiency, throughput and/or product quality. Industrial applications can be characterized as large, involving complex particulate behaviour in typically complex geometries. In this paper, with a series of examples, we will explore the breadth of large scale modelling of industrial processes that is currently possible. Few of these applications will be examined in more detail to show how insights into the fundamentals of these processes can be gained through DEM modelling. Some examples of our collaborative validation efforts will also be described.

The purpose of this paper is to show how particle scale simulation of industrial particle flows using DEM (discrete element method) offers the opportunity for better understanding of the flow dynamics leading to improvements in equipment design and operation.

The paper explores the breadth of industrial applications that are now possible with a series of case studies.

The paper finds that the inclusion of cohesion, coupling to other physics such fluids, and its use in bubbly and reacting flows are becoming increasingly viable. Challenges remain in developing models that balance the depth of the physics with the computational expense that is affordable and in the development of measurement and characterization processes to provide this expanding array of input data required. Steadily increasing computer power has seen model sizes grow from thousands of particles to many millions over the last decade, which steadily increases the range of applications that can be modelled and the complexity of the physics that can be well represented.

The paper shows how better understanding of the flow dynamics leading to improvements in equipment design and operation can potentially lead to large increases in equipment and process efficiency, throughput and/or product quality. Industrial applications can be characterised as large, involving complex particulate behaviour in typically complex geometries. The critical importance of particle shape on the behaviour of granular systems is demonstrated. Shape needs to be adequately represented in order to obtain quantitative predictive accuracy for these systems.

The purpose of this paper is to examine the effects of Darcian drag force and radiation-conduction on unsteady two-dimensional magnetohydrodynamic flow of viscous, electrically conducting and Newtonian fluid over a vertical plate adjacent to a Darcian regime in presence of thermal radiation and transversal magnetic field. A well-tested, numerically stable Crank-Nicolson finite-difference procedure is employed for the conservation equations. Excellent agreement is obtained for numerical solutions with previously published work.

In this investigation, an efficient, accurate, extensively validated and unconditionally stable finite-difference scheme based on the Crank-Nicolson model is developed to solve the governing coupled, non-linear partial differential equations. The accuracy and effectiveness of the method are demonstrated.

Different numerical results are obtained and presented graphically to explain the effect of various physical parameters on the velocity and temperature profiles, local, as well as average, skin friction and Nusselt number. It is found that, with a rise in Darcian drag force, flow velocity and temperature are reduced, but increased for all times. Both average and local skin frictions are reduced considerably with an increase in Darcian drag force, but reversed behavior is observed for the local Nusselt number. Increasing the thermal radiation effects accelerated the flow velocity as well as the fluid temperature and wall local skin friction in a saturated porous medium, but effectively reduced the local Nusselt number and average Nusselt number at the wall. Comparison with previously published works in the limits shows excellent agreement.

The analysis is valid for unsteady, two-dimensional laminar flow of an optically thick no-gray gas, electrically conducting, and Newtonian fluid past an isothermal vertical surface adjacent to the Darcian regime with variable surface temperature. An extension to three-dimensional flow case is left for future work.

Practical interest of such study includes applications in electromagnetic lubrication, boundary cooling, bio-physical systems and in many branches of engineering and science. It is well known that the effect of thermal radiation is important in space technology and high temperature processes. Thermal radiation also plays an important role in controlling heat transfer process in polymer processing industry.

Numerical results generated by a highly efficient finite‐difference method (originated by Keller for aerodynamical flows at the California Institute of Technology in 1970), and a robust double shooting Runge‐Kutta‐Merson scheme are presented for the boundary layer equations representing the convection flow of a viscous incompressible fluid past a hot vertical flat plate embedded in a non‐Darcy porous medium. Viscous dissipation due to mechanical work is included in the temperature field equation. The computations for both solution techniques are compared at the leading edge (ξ = 0.0) and found to be in excellent agreement. The effects of the viscous heating parameter (Ec), thermal conductivity ratio (λ) and a Darcy porous parameter (Re/GrDa) on the fluid velocities, temperatures, local shear stress and wall heat transfer rate are discussed with applications to geothermal and industrial flows.

This study aims to deal with the laminar flow owing to rough rotating disk in the existence of vertical magnetic field and partial slip effects. The aim is to resolve heat transfer problem in the existence of non-linear radiative flux and thermal slip effects. The study also analyzes the mass transfer process when the flow field contains chemically reacting species.

Modified von-Kármán transformations are applied to change the conservation equations into similar forms. The transformed equations are treated by a convenient shooting method and by contemporary built in routine bvp4c of MATLAB.

The numerical solutions are used to address the role of main ingredients of the problem, namely, wall roughness, radiation and chemical reaction on the flow fields.

Temperature profiles are considerably affected by a parameter measuring wall to ambient temperature ratio. Furthermore, behavior of concentration field is highly influenced by the reaction rate of the diffusing species.

The concept of non-linear radiation in chemically reactive flow over a rotating disk is just introduced here.

This paper aims to intend for investigating the influence of an inclined magnetic field on the peristaltic flow of a couple stress fluid through an inclined channel with heat and mass transfer.

Long wavelength and low Reynolds number methodology is actualized for simplifying the highly nonlinear equations. Mathematical expressions of axial velocity, pressure gradient and volume flow rate are obtained. Pressure rise, frictional force and pumping phenomenon are portrayed and symbolized graphically. Exact and numerical solutions have been carried out. The computed results are presented graphically for various embedded parameters. Temperature and concentration profile are also scrutinized and sketched.

Results from the current study concluded that the fluid motion can be enhanced by increasing the inclination of both the magnetic field and the channel.

The elemental characteristics of this analysis is a complete interpretation of the influence of couple stress parameter and inclination of magnetic field on the velocity, pressure gradient, pressure rise and frictional forces.

The purpose of this paper is to consider unsteady free convection flow of a dissipative fluid past an exponentially accelerated infinite vertical porous plate in the presence of Newtonian heating and mass diffusion.

The problem is governed by coupled non-linear partial differential equations with appropriate boundary conditions. A Galerkin finite element numerical solution is developed to solve the resulting well-posed two-point boundary value problem. It is a powerful, stable technique which provides excellent convergence and versatility in accommodating coupled systems of ordinary and partial differential equations.

It is found that the skin friction coefficient increases with increases in either of the Eckert number, thermal Grashof number, mass Grashof number or time whereas it decreases with increases in either of the suction parameter, Schmidt number or the acceleration parameter for both air and water. The skin friction coefficient is also found to decrease with increases in the values of the Prandtl number. In addition, it is found that the rate of heat transfer increases with an increase in the suction parameter and decreases with an increase in the Eckert number for both air and water. Lastly, it is found that the rate of heat transfer increases with increasing values of the Prandtl number and decreases with increasing time for all values of the Prandtl number.

The present study has considered only Newtonian fluids. Future studies will address non-Newtonian liquids.

A very useful source of information for researchers on the subject of free convective flow over the surface when the rate of heat transfer from the surface is proportional to the local surface temperature.

This paper is relatively original and illustrates the effects of viscous dissipation on free convective flow past an exponentially accelerated infinite vertical porous plate with Newtonian heating and mass diffusion.

Modelling accurately the transient behaviour of natural convection flow in enclosures been a challenging task because of a variety of numerical errors which have limited achieving the higher order temporal accuracy. A fourth-order accurate finite difference method in both space and time is proposed to overcome these numerical errors and accurately model the transient behaviour of natural convection flow in enclosures using vorticity–streamfunction formulation.

Fourth-order wide stencil formula with appropriate one-sided difference extrapolation technique near the boundary is used for spatial discretisation, and classical fourth-order Runge–Kutta scheme is applied for transient term discretisation. The proposed method is applied on two transient case studies, i.e. convection–diffusion of a Gaussian Pulse and Taylor Vortex flow having analytical solution.

Error magnitude comparison and rate of convergence analysis of the proposed method with these analytical solutions establish fourth-order accuracy and prove the ability of the proposed method to truly capture the transient behaviour of incompressible flow. Also, to test the transient natural convection flow behaviour, the algorithm is tested on differentially heated square cavity at high Rayleigh number in the range of 10³-10⁸, followed by studying the transient periodic behaviour in a differentially heated vertical cavity of aspect ratio 8:1. An excellent comparison is obtained with standard benchmark results.

The developed method is applied on 2D enclosures; however, the present methodology can be extended to 3D enclosures using velocity–vorticity formulations which shall be explored in future.

The proposed methodology to achieve fourth-order accurate transient simulation of natural convection flows is novel, to the best of the authors’ knowledge. Stable fourth-order vorticity boundary conditions are derived for boundary and external boundary regions. The selected case studies for comparison demonstrate not only the fourth-order accuracy but also the considerable reduction in error magnitude by increasing the temporal accuracy. Also, this study provides novel benchmark results at five different locations within the differentially heated vertical cavity of aspect ratio 8:1 for future comparison studies.