scholarly journals Numerical Study of Flow and Particle Deposition in Wall-Flow Filters with Intact or Damaged Exit

Fluids ◽  
2019 ◽  
Vol 4 (4) ◽  
pp. 201 ◽  
Author(s):  
Chris D. Dritselis ◽  
Fotini Tzorbatzoglou ◽  
Marios Mastrokalos ◽  
Onoufrios Haralampous

We examine the time-dependent three-dimensional gas-particle flow in an intact wall-flow filter consisting of channels alternatively plugged at each end and a partially damaged filter in which the rear plugs are removed. Our focus is placed on highlighting the differences in the flow pattern and the deposition process between the two geometries. The Navier–Stokes equations are solved for the fluid flow coupled with a Brinkman/Forchheimmer model in order to simulate the flow in the porous walls and plugs. Discrete particle simulation is utilized to determine the nanoparticle trajectories. Using this scheme, we are able to characterize the main features of the flow fields developing in the intact and damaged filters with respect to the Reynolds number and identify those affecting the transport and deposition of particles that have three representative response times. We present fluid velocity iso-contours, which describe the flow regimes inside the channels, as well as in regions upstream and downstream of them. We provide evidence of local recirculating bubbles at the entrance of the channels and after their exit, whereas back-flow occurs in front of the rear plugs of the intact channels. We show that the flow leaves the channels as strong jets that may break up for certain flow parameters, leading to turbulence with features that depend on the presence of the rear plugs. The removal of the rear plugs affects the flow distribution, which, in turn influences the flow rates along the channels and through the walls. We describe the particle trajectories and the topology of deposited particles and show that particles follow closely the streamlines, which may cross the surface of permeable walls for both flow configurations. The distribution of deposited particles resembles the spatial variation of the through-wall flow rate, exhibiting two peak values at both ends of the intact filter channel, and one local maximum near the entrance of the damaged filter channel that is diminished at the exit. We also investigate in detail the particle deposition on the frontal face and indicate that particle accumulation at the edges of the entrance is favored for particles with low response times in flows with high fluid mass rates for both intact and damaged filters. Finally, we examine the filtration efficiency for the defective channels without rear plugs and show that fewer particles are captured as the Reynolds number is increased. A smaller reduction of the filtration efficiency is also predicted with increasing particle size.

Author(s):  
Azita Soleymani ◽  
Eveliina Takasuo ◽  
Piroz Zamankhan ◽  
William Polashenski

Results are presented from a numerical study examining the flow of a viscous, incompressible fluid through random packing of nonoverlapping spheres at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), spanning a wide range of flow conditions for porous media. By using a laminar model including inertial terms and assuming rough walls, numerical solutions of the Navier-Stokes equations in three-dimensional porous packed beds resulted in dimensionless pressure drops in excellent agreement with those reported in a previous study (Fand et al., 1987). This observation suggests that no transition to turbulence could occur in the range of Reynolds number studied. For flows in the Forchheimer regime, numerical results are presented of the lateral dispersivity of solute continuously injected into a three-dimensional bounded granular bed at moderate Peclet numbers. Lateral fluid dispersion coefficients are calculated by comparing the concentration profiles obtained from numerical and analytical methods. Comparing the present numerical results with data available in the literature, no evidence has been found to support the speculations by others for a transition from laminar to turbulent regimes in porous media at a critical Reynolds number.


2019 ◽  
Vol 29 (10) ◽  
pp. 3908-3937 ◽  
Author(s):  
Younes Menni ◽  
Ahmed Azzi ◽  
Ali J. Chamkha ◽  
Souad Harmand

Purpose The purpose of this paper is to carry out a numerical study on the dynamic and thermal behavior of a fluid with a constant property and flowing turbulently through a two-dimensional horizontal rectangular channel. The upper surface was put in a constant temperature condition, while the lower one was thermally insulated. Two transverse, solid-type obstacles, having different shapes, i.e. flat rectangular and V-shaped, were inserted into the channel and fixed to the top and bottom walls of the channel, in a periodically staggered manner to force vortices to improve the mixing, and consequently the heat transfer. The flat rectangular obstacle was put in the first position and was placed on the hot top wall of the channel. However, the second V-shaped obstacle was placed on the insulated bottom wall, at an attack angle of 45°; its position was varied to find the optimum configuration for optimal heat transfer. Design/methodology/approach The fluid is considered Newtonian, incompressible with constant properties. The Reynolds averaged Navier–Stokes equations, along with the standard k-epsilon turbulence model and the energy equation, are used to control the channel flow model. The finite volume method is used to integrate all the equations in two-dimensions; the commercial CFD software FLUENT along with the SIMPLE-algorithm is used for pressure-velocity coupling. Various values of the Reynolds number and obstacle spacing were selected to perform the numerical runs, using air as the working medium. Findings The channel containing the flat fin and the 45° V-shaped baffle with a large Reynolds number gave higher heat transfer and friction loss than the one with a smaller Reynolds number. Also, short separation distances between obstacles provided higher values of the ratios Nu/Nu0 and f/f0 and a larger thermal enhancement factor (TEF) than do larger distances. Originality/value This is an original work, as it uses a novel method for the improvement of heat transfer in completely new flow geometry.


Author(s):  
Arezou Jafari ◽  
S. Mohammad Mousavi

Numerical study of flow through random packing of non-overlapping spheres in a cylindrical geometry is investigated. Dimensionless pressure drop has been studied for a fluid through the porous media at moderate Reynolds numbers (based on pore permeability and interstitial fluid velocity), and numerical solution of Navier-Stokes equations in three dimensional porous packed bed illustrated in excellent agreement with those reported by Macdonald [1979] in the range of Reynolds number studied. The results compare to the previous work (Soleymani et al., 2002) show more accurate conclusion because the problem of channeling in a duct geometry. By injection of solute into the system, the dispersivity over a wide range of flow rate has been investigated. It is shown that the lateral fluid dispersion coefficients can be calculated by comparing the concentration profiles of solute obtained by numerical simulations and those derived analytically by solving the macroscopic dispersion equation for the present geometry.


2009 ◽  
Vol 640 ◽  
pp. 235-264 ◽  
Author(s):  
SHAOPING QUAN ◽  
DAVID P. SCHMIDT ◽  
JINSONG HUA ◽  
JING LOU

The relaxation and breakup of an elongated droplet in a viscous and initially quiescent fluid is studied by solving the full Navier–Stokes equations using a three-dimensional finite volume method coupled with a moving mesh interface tracking (MMIT) scheme to locate the interface. The two fluids are assumed incompressible and immiscible. The interface is represented as a surface triangle mesh with zero thickness that moves with the fluid. Therefore, the jump and continuity conditions across the interface are implemented directly, without any smoothing of the fluid properties. Mesh adaptations on a tetrahedral mesh are employed to permit large deformation and to capture the changing curvature. Mesh separation is implemented to allow pinch-off. The detailed investigations of the relaxation and breakup process are presented in a more general flow regime compared to the previous works by Stone & Leal (J. Fluid Mech., vol. 198, 1989, p. 399) and Tong & Wang (Phys. Fluids, vol. 19, 2007, 092101), including the flow field of the both phases. The simulation results reveal that the vortex rings due to the interface motion and the conservation of mass play an important role in the relaxation and pinch-off process. The vortex rings are created and collapsed during the process. The effects of viscosity ratio, density ratio and length ratio on the relaxation and breakup are studied. The simulations indicate that the fluid velocity field and the neck shape are distinctly different for viscosity ratios larger and smaller than O(1), and thus a different end-pinching mechanism is observed for each regime. The length ratio also significantly affects the relaxation process and the velocity distributions, but not the neck shape. The influence of the density ratio on the relaxation and breakup process is minimal. However, the droplet evolution is retarded due to the large density of the suspending flow. The formation of a satellite droplet is observed, and the volume of the satellite droplet depends strongly on the length ratio and the viscosity ratio.


1988 ◽  
Vol 110 (4) ◽  
pp. 446-452 ◽  
Author(s):  
C. L. Rumsey

The application of an upwind-biased implicit approximate factorization Navier-Stokes algorithm to the unsteady impulsive start-up flow over a circular cylinder at Reynolds number 1200 is described. The complete form of the compressible Navier-Stokes equations is used, and the algorithm is second-order accurate in both space and time. The drag on the cylinder is computed for early times in the start-up flow. The value of the local maximum drag as well as the time at which it occurs are predicted and compared to another computational result and experiment. The development with time of the shape and size of the separated vortical flow region is computed, as well as the time-variation of several boundary layer parameters and profile shapes. Computations, in general, show excellent agreement with experiment, although the present method predicts a more rapid onset of reversed flow on the cylinder than evidenced in experiment. The effect of grid density on the development of the unsteady flow is also shown.


1996 ◽  
Vol 306 ◽  
pp. 293-323 ◽  
Author(s):  
Vadim Borue ◽  
Steven A. Orszag

High-resolution numerical simulations (with up to 2563 modes) are performed for three-dimensional flow driven by the large-scale constant force fy = F cos(x) in a periodic box of size L = 2π (Kolmogorov flow). High Reynolds number is attained by solving the Navier-Stokes equations with hyperviscosity (-1)h+1Δh (h = 8). It is shown that the mean velocity profile of Kolmogorov flow is nearly independent of Reynolds number and has the ‘laminar’ form vy = V cos(x) with a nearly constant eddy viscosity. Nevertheless, the flow is highly turbulent and intermittent even at large scales. The turbulent intensities, energy dissipation rate and various terms in the energy balance equation have the simple coordinate dependence a + b cos(2x) (with a, b constants). This makes Kolmogorov flow a good model to explore the applicability of turbulence transport approximations in open time-dependent flows. It turns out that the standard expression for effective (eddy) viscosity used in K-[Escr ] transport models overpredicts the effective viscosity in regions of high shear rate and should be modified to account for the non-equilibrium character of the flow. Also at large scales the flow is anisotropic but for large Reynolds number the flow is isotropic at small scales. The important problem of local isotropy is systematically studied by measuring longitudinal and transverse components of the energy spectra and crosscorrelation spectra of velocities and velocity-pressure-gradient spectra. Cross-spectra which should vanish in the case of isotropic turbulence decay only algebraically but somewhat faster than corresponding isotropic correlations. It is verified that the pressure plays a crucial role in making the flow locally isotropic. It is demonstrated that anisotropic large-scale flow may be considered locally isotropic at scales which are approximately ten times smaller than the scale of the flow.


Author(s):  
S. Bhattacharyya ◽  
D. K. Maiti

Numerical study on the wake behind a square cylinder placed parallel to a wall has been made. Flow has been investigated in the laminar Reynolds number (based on the cylinder length) range. We have studied the flow field for different values of the non-dimensional gap length between cylinder and the wall. The case when the cylinder is placed on the wall has also been considered. The governing unsteady Navier-Stokes equations are discretised through the finite volume method on staggered grid system. A SIMPLER type of algorithm has been used to compute the discretised equations iteratively. Vortex shedding has been found to be influenced by the wall. Vortex shedding suppression occurs beyond a critical value of the gap length. Due to the shear, the drag experienced by the cylinder is found to increase with the reduction of gap length. The flow is found to be steady when the cylinder is placed on the wall at a range of Reynolds number.


Author(s):  
Eiyad Abu-Nada ◽  
AbdelSalam Al-Sarkhi

The current research investigates the effect of Reynolds number and Knudsen number on the coefficient of skin friction and reattachment length for a micro-scale fluid flow over a step mounted on a lower wall of a micro-channel. Five Reynolds numbers are studied Re = 1, 10, 25, 50, and 75 and the Knudsen number is varied from 1×10−3 to 0.1. Finite difference method with non-uniform grid is used to solve the incompressible Navier Stokes equations accompanied with velocity slip boundary condition. As Knudsen number (Kn) decreases the magnitude of modified local shear stress (1/2 Cf Re), on the upper wall of the channel, increases. In the circulation zone behind the step and for the case of high Reynolds number (Re = 50 and Re = 75) the modified local shear stress increases as Knudsen number increases. Results show that the modified total skin friction (1/2 CD Re) decreases as the Knudsen number increases. The modified total skin friction drops significantly with Knudsen number for Kn >= 1 × 10−2. However, (1/2 CD Re) is relatively independent of Knudsen number for Kn < 1 × 10−2. Finally, for 1 × 10−2 < Kn < 0.1, as Knudsen number increase the reattachment length increases.


1998 ◽  
Vol 120 (2) ◽  
pp. 345-353 ◽  
Author(s):  
Chao-Tsung Hsiao ◽  
Laura L. Pauley

The flow over a finite-span hydrofoil creating a tip vortex was numerically studied by computing the full Navier-Stokes equations. A good agreement in pressure distribution and oil flow pattern was achieved between the numerical solution and available experimental data. The steady-state roll-up process of the tip vortex was described in detail from the numerical results. The effect of the angle of attack, the Reynolds number, and the hydrofoil planform on the tip vortex was investigated. The axial and tangential velocities within the tip-vortex core in the near-field wake region were greatly influenced by the angle of attack. A jet-like profile in the axial velocity was found within the tip-vortex core at high angle of attack, while a wake-like profile in the axial velocity was found at low angle of attack. Increasing the Reynolds number was found to increase the maximum axial velocity, but only had a slight impact on the tangential velocity. Finally, a swept hydrofoil planform was found to attenuate the strength of the tip vortex due to the low-momentum boundary layer traveling into the tip vortex on the suction side.


2019 ◽  
Vol 29 (5) ◽  
pp. 1629-1646 ◽  
Author(s):  
Arya Ghiasi ◽  
Seyed Esmaeil Razavi ◽  
Abel Rouboa ◽  
Omid Mahian

Purpose This study aims to investigate the effect of the simultaneous usage of active and passive methods (which in this case are rotational oscillation and attached splitter plate, respectively) on the flow and temperature fields to find an optimum situation which this combination results in heat transfer increment and drag reduction. Design/methodology/approach The method of the solution was based on finite volume discretization of Navier–Stokes equations. A dynamic grid is coupled with the solver by the arbitrary Lagrangian–Eulerian (ALE) formulation for modeling cylinder oscillation. Parametric studies were performed by altering oscillation frequency, splitter plate length and Reynolds number. Findings Oscillation in different frequencies was found to be complicated. Higher frequencies provide more heat transfer, but in the lock-on region, they bring remarkable increment to the drag coefficient. It was observed that simultaneous usage of oscillation and splitter plate may have both positive and negative effects on drag reduction and heat transfer increment. Finally F = 2 and L = 0.5 were chosen as an optimum combination. Originality/value In this study, the laminar incompressible flow and heat transfer from a confined rotationally oscillating circular cylinder with an attached splitter plate are investigated. Parametric studies are performed by changing oscillation frequency, splitter plate length and Reynolds number.


Sign in / Sign up

Export Citation Format

Share Document