scholarly journals Numerical Simulations of the Impact and Spreading of a Particulate Drop on a Solid Substrate

2012 ◽  
Vol 2012 ◽  
pp. 1-10
Author(s):  
Hyun Jun Jeong ◽  
Wook Ryol Hwang ◽  
Chongyoup Kim

We present two-dimensional numerical simulations of the impact and spreading of a droplet containing a number of small particles on a flat solid surface, just after hitting the solid surface, to understand particle effects on spreading dynamics of a particle-laden droplet for the application to the industrial inkjet printing process. The Navier-Stokes equation is solved by a finite-element-based computational scheme that employs the level-set method for the accurate interface description between the drop fluid and air and a fictitious domain method for suspended particles to account for full hydrodynamic interaction. Focusing on the particle effect on droplet spreading and recoil behaviors, we report that suspended particles suppress the droplet oscillation and deformation, by investigating the drop deformations for various Reynolds numbers. This suppressed oscillatory behavior of the particulate droplet has been interpreted with the enhanced energy dissipation due to the presence of particles.

2019 ◽  
Vol 88 (2) ◽  
pp. 20901 ◽  
Author(s):  
Mouloud Driouche ◽  
Tahar Rezoug ◽  
Mohammed El Ganaoui

The substrate melting can significantly improve the properties of plasma spray coatings. Indeed the adhesion of the projected particles to the substrate can be ameliorated by the substrate melting. In this article, a numerical model is developed to study the dynamics of fluids and heat transfer with liquid/solid phase change during impact of a fully melted alumina particle on an aluminum solid substrate, taking into account of the substrate melting. The model is based on solving the Navier-Stokes and energy equations with liquid / solid phase change. These equations are coupled with the fluid of volume method (VOF), to follow the free surface of the particle during its spreading and solidification. The finite volume method is used to discretize the equations in a 2D axisymmetric domain. A comparison with the published experimental results was carried out to validate this numerical model. Simulations were performed for different initial droplet diameters to study its effect on droplet spreading as well as on substrate melting. It has been observed that the substrate melting begins before the droplet spreads completely; the substrate melting reaches its maximum when the droplet is close to its total solidification. Droplet spreading and substrate melting are more important for large sizes droplets.


2007 ◽  
Vol 573 ◽  
pp. 171-190 ◽  
Author(s):  
A. DIPANKAR ◽  
T. K. SENGUPTA ◽  
S. B. TALLA

Vortex shedding behind a cylinder can be controlled by placing another small cylinder behind it, at low Reynolds numbers. This has been demonstrated experimentally by Strykowski & Sreenivasan (J. Fluid Mech. vol. 218, 1990, p. 74). These authors also provided preliminary numerical results, modelling the control cylinder by the innovative application of boundary conditions on some selective nodes. There are no other computational and theoretical studies that have explored the physical mechanism. In the present work, using an over-set grid method, we report and verify numerically the experimental results for flow past a pair of cylinders. Apart from providing an accurate solution of the Navier–Stokes equation, we also employ an energy-based receptivity analysis method to discuss some aspects of the physical mechanism behind vortex shedding and its control. These results are compared with the flow picture developed using a dynamical system approach based on the proper orthogonal decomposition (POD) technique.


1981 ◽  
Vol 103 (4) ◽  
pp. 785-790 ◽  
Author(s):  
J. H. Masliyah ◽  
K. Nandakumar

The Navier-Stokes equation in a rotating frame of reference is solved numerically to obtain the flow field for a steady, fully developed laminar flow of a Newtonian fluid in a twisted tube having a square cross-section. The macroscopic force and energy balance equations and the viscous dissipation term are presented in terms of variables in a rotating reference frame. The computed values of friction factor are presented for dimensionless twist ratios, (i.e., length of tube over a rotation of π radians normalized with respect to half the width of tube) of 20, 10, 5 and 2.5 and for Reynolds numbers up to 2000. The qualitative nature of the axial velocity profile was observed to be unaffected by the swirling motion. The secondary motion was found to be most important near the wall.


Author(s):  
Carl E. Rathmann

For well over 150 years now, theoreticians and practitioners have been developing and teaching students easily visualized models of fluid behavior that distinguish between the laminar and turbulent fluid regimes. Because of an emphasis on applications, perhaps insufficient attention has been paid to actually understanding the mechanisms by which fluids transition between these regimes. Summarized in this paper is the product of four decades of research into the sources of these mechanisms, at least one of which is a direct consequence of the non-linear terms of the Navier-Stokes equation. A scheme utilizing chaotic dynamic effects that become dominant only for sufficiently high Reynolds numbers is explored. This paper is designed to be of interest to faculty in the engineering, chemistry, physics, biology and mathematics disciplines as well as to practitioners in these and related applications.


2011 ◽  
Vol 685 ◽  
pp. 461-494 ◽  
Author(s):  
Alain Merlen ◽  
Christophe Frankiewicz

AbstractThe flow around a cylinder rolling or sliding on a wall was investigated analytically and numerically for small Reynolds numbers, where the flow is known to be two-dimensional and steady. Both prograde and retrograde rotation were analytically solved, in the Stokes regime, giving the values of forces and torque and a complete description of the flow. However, solving Navier–Stokes equation, a rotation of the cylinder near the wall necessarily induces a cavitation bubble in the nip if the fluid is a liquid, or compressible effects, if it is a gas. Therefore, an infinite lift force is generated, disconnecting the cylinder from the wall. The flow inside this interstice was then solved under the lubrication assumptions and fully described for a completely flooded interstice. Numerical results extend the analysis to higher Reynolds number. Finally, the effect of the upstream pressure on the onset of cavitation is studied, giving the initial location of the phenomenon and the relation between the upstream pressure and the flow rate in the interstice. It is shown that the flow in the interstice must become three-dimensional when cavitation takes place.


2014 ◽  
Vol 24 (09) ◽  
pp. 1450110 ◽  
Author(s):  
Rodica Borcia ◽  
Michael Bestehorn

We use a phase field model which couples the generalized Navier–Stokes equation (including the Korteweg stress tensor) with the continuity equation for studying nonlinear pattern formation on the surface of a liquid film under (linear and circular) lateral harmonic vibrations at the solid substrate. First, we prove the thermodynamic consistency of our phase field model. Next, we present computer simulations in three spatial dimensions. We illustrate nonequilibrium patterns at the instability onset, confirming in this way the results recently reported in Phys. Rev. E 88, 023025 (2013). The lateral profiles of the deflected surface are compared with those reported in J. Fluid Mech. 686, 409 (2011) for Faraday instability excited by vertical harmonic vibrations at the bottom plate.


1993 ◽  
Vol 256 ◽  
pp. 615-646 ◽  
Author(s):  
Paolo Orlandi ◽  
Roberto Verzicco

Accurate numerical simulations of vortex rings impinging on flat boundaries revealed the same features observed in experiments. The results for the impact with a free-slip wall compared very well with previous numerical simulations that used spectral methods, and were also in qualitative agreement with experiments. The present simulation is mainly devoted to studying the more realistic case of rings interacting with a no-slip wall, experimentally studied by Walker et al. (1987). All the Reynolds numbers studied showed a very good agreement between experiments and simulations, and, at Rev > 1000 the ejection of a new ring from the wall was seen. Axisymmetric simulations demonstrated that vortex pairing is the physical mechanism producing the ejection of the new ring. Three-dimensional simulations were also performed to investigate the effects of azimuthal instabilities. These simulations have confirmed that high-wavenumber instabilities originate in the compression phase of the secondary ring within the primary one. The large instability of the secondary ring has been explained by analysis of the rate-of-strain tensor and vorticity alignment. The differences between passive scalars and the vorticity field have been also investigated.


2017 ◽  
Vol 62 (1) ◽  
pp. 16 ◽  
Author(s):  
Younes Menni ◽  
Ahmed Azzi

A computational fluid dynamic analysis of thermal and aerodynamic fields for an incompressible steady-state flow of a Newtonian fluid through a two-dimensional horizontal rectangular section channel with upper and lower wall-attached, vertical, staggered, transverse, cascaded rectangular-triangular (CRT), solid-type baffles is carried out in the present paper using the Commercial, Computational Fluid Dynamics, software FLUENT. The flow model is governed by the Reynolds averaged Navier-Stokes (RANS) equations with the SST k-ω turbulence model and the energy equation. The finite volume method (FVM) with the SIMPLE-discretization algorithm is applied for the solution of the problem. The computations are carried out in the turbulent regime for different Reynolds numbers. In this study, thermo-aeraulic fields, dimensionless axial profiles of velocity, skin friction coefficients, local and average heat transfer coefficients, and thermal enhancement factor were investigated, at constant surface temperature condition along the heated upper wall of the channel, for all the geometry under investigation and chosen for various stations. The impact of the cascaded rectangular-triangular geometry of the baffle on the thermal and dynamic behavior of air is shown and this in comparing the data of this obstacle type with those of the simple flat rectangular-shaped baffle. This CFD analysis can be a real application in the field of heat exchangers, solar air collectors, and electronic equipments.


2012 ◽  
Author(s):  
El–Harbawi M. ◽  
Mustapha S. ◽  
Idris A. ◽  
Jameel A.T. ◽  
T.G. Chuah

Tenaga lebihan bebas Gibb per unit luas (ΔG) bagi suatu saput nipis yang tulen, tanpa cas pada satu sokongan hanya diterbitkan daripada interaksi antara molekul apolar (Lifshitz - van der Waals) dan polar (acid-bes). Interaksi polar dapat dinyatakan secara berasingan sebagai tekanan penghidratan, interaksi hidrofobik dan intekrasi asid-bes. Jumlah tenaga lebihan bebas Gibb per unit luas (ΔG) bagi satu saput nipis pada satu lapisan bergantung kepada ketebalan saput, dan angkali penyebaran apolar dan polar untuk sistem tersebut. Saput dimodelkan sebagai cecair Newtonian 2-dimensi dengan ketumpatan tetap, ρ, dan kelikatan, μ, mengalir di atas satu permukaan mengufuk. Ketebalan purata saput cecair, h0, dianggapkan cukup nipis supaya kesan graviti boleh diabaikan dan melambung atas oleh gas pasif dan menyambung secara lateral hingga infiniti (Model dua-dimensi). Daya badan yang disebut dalam persamaan Navier-Stokes adalah diubahsuaikan oleh kandungan interaksi antara molekul lebihan (daya apolar dan polar) antara saput bendalir dan permukaan pepejal bergantung kepada daya-daya apolar dan polar. Persamaan Navier-Stokes telah diubahsuai dengan keadaan sempadan berkaitan diselesaikan bawah kaedah pendekatan gelombong panjang untuk memperolehi persamaan perkembangan tak-linear bagi saput antara permukaan. Daya apolar dan polar telah didapati memainkan peranan pencirian atas saput nipis dan kesan utama pada sifat tenaga lebihan bebas, kadar penubuhan, kadar penubuhan maksimum, gelombong neutral, nombor gelombang yang dominan, jarak gelombong dominan dan masa memecah. Oleh demikian, teori linear adalah kurang sesuai untuk menyatakan kestabilan pencirian selaput. Kata kunci: Daya apolar, daya polar, kestabilan linear, kadar pertumbuhan, masa memecah The total excess free energy per unit area (ΔG) of a pure, uncharged thin film on a support is solely derived from the apolar (Lifshitz - van der Waals) and polar (acid-base) intermolecular interactions. Polar interactions are variously described as the hydration pressure, hydrophobic interaction and acid base interaction. The total free excess energy (per unit area) of a thin film on a substrate depends on the film thickness, and the apolar and polar spreading coefficients for the system. The film is modelled as a two-dimensional Newtonian liquid of constant density, ρ and viscosity, μ, flowing on a horizontal plane. The liquid film of mean thickness, h0, is assumed to be thin enough to neglect the gravity effect and bounded above by a passive gas and laterally extends to infinity (two-dimensional model). The body force term in the Navier-Stokes equation is modified by the inclusion of excess intermolecular interactions (apolar and polar forces) between fluid film and the solid surface owing to apolar and polar forces. The modified Navier-Stokes equation with associated boundary conditions is solved under long wave approximation method to obtain a nonlinear equation of evolution of the film interface. The apolar and polar forces were found to play the dominant role in characteristic of thin films and the main effect on the behavior of the excess free energy, growth rate, maximum growth rate, neutral wave, dominant wavenumber, dominant wavelength and rupture time. Hence, the linear theory is inadequate to describe the stability characteristics of films. Key words: Apolar force, polar force, linear stability, growth rate, rupture time


2015 ◽  
Vol 99 ◽  
pp. 15-25 ◽  
Author(s):  
Zhe Hu ◽  
Wenyong Tang ◽  
Hongxiang Xue ◽  
Xiaoying Zhang ◽  
Jinting Guo

Sign in / Sign up

Export Citation Format

Share Document