No‐slip consistent immersed boundary particle tracking to simulate impaction filtration in porous media

In this paper, we present a new method for simulating the motion of a disperse particle phase in a carrier gas through porous media. We assume a sufficiently dilute particle‐laden flow and compute, independently of the disperse phase, the steady laminar fluid velocity using the immersed boundary method. Given the velocity of the carrier gas, the equations of motion for the particles experiencing the Stokes drag force are solved to determine their trajectories. The ‘no‐slip consistent’ particle tracking algorithm avoids possible numerical filtration of very small particles due to the nonzero velocity field at the solid–fluid interface introduced by the immersed boundary method. This physically consistent tracking allows a reliable estimation of the filtration efficiency of porous filters due to inertial impaction. We illustrate and test our new approach for model porous media consisting of a structured array of aligned rectangular fibers, arranged in line and staggered. In the staggered geometry, the effect of the residual velocity at the solid–fluid interface is significant for particles with low inertia. Without adopting the developed no‐slip consistent numerical method, an artificial numerical filtration is observed, which becomes dominant for small enough particles. For both the in line and the staggered geometries, the filtration rate depends quite strongly and non monotonically on the particle inertia. This is expressed most clearly in the staggered arrangement in which a very strong increase in the filtration efficiency is observed at a well‐defined critical droplet size, corresponding to a qualitative change in the dominant particle paths in the porous medium. Copyright © 2013 John Wiley & Sons, Ltd.     

Computing the flow around a moving profile by numerically integrating the Navier–Stokes equations

A vorticity velocity formulation is proposed for the solution of the equations for viscous flow around a moving profile. A non-inertial reference frame is used and the velocities are computed from a Poincaré integral formula. The studies are directed towards the need to understand helicopter blade aerodynamics. Worked examples are given which validate the method and programme for laminar flows, at least for low Reynolds numbers. © 1998 John Wiley & Sons, Ltd.     

A multi‐energy‐level lattice Boltzmann model for the compressible Navier–Stokes equations

In this paper, we propose a new lattice Boltzmann model for the compressible Navier–Stokes equations. The new model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. As the 25‐bit model, we obtained the equilibrium distribution functions and the compressible Navier–Stokes equations with the second accuracy of the truncation errors. The numerical examples show that the model can be used to simulate the shock waves, contact discontinuities and supersonic flows around circular cylinder. The numerical results are compared with those obtained by traditional method. Copyright © 2007 John Wiley & Sons, Ltd.     

The Lagrangian approach of advective term treatment and its application to the solution of Navier—Stokes equations

A one-dimensional transport test applied to some conventional advective Eulerian schemes shows that linear stability analyses do not guarantee the actual performances of these schemes. When adopting the Lagrangian approach, the main problem raised in the numerical treatment of advective terms is a problem of interpolation or restitution of the transported function shape from discrete data. Several interpolation methods are tested. Some of them give excellent results and these methods are then extended to multi-dimensional cases. The Lagrangian formulation of the advection term permits an easy solution to the Navier-Stokes equations in primitive variables V, p, by a finite difference scheme, explicit in advection and implicit in diffusion. As an illustration steady state laminar flow behind a sudden enlargement is analysed using an upwind differencing scheme and a Lagrangian scheme. The importance of the choice of the advective scheme in computer programs for industrial application is clearly apparent in this example.     

Variable density flow in porous media. A study by means of pore level numerical simulations

This paper is devoted to the modelling of isothermal low Reynolds and Mach numbers transient compressible flow through porous media. Traditionally, this type of flow at the macroscopic level is described by the classical Darcy's law combined with a mass balance that includes the transient term. This model is called the ‘classic model’. The aim of this paper is to explore the validity of this classic model. To this end, the flow of an ideal gas is considered within two‐dimensional model porous media. The flow is due to the imposed pressure variations at the outlet of the fluid domain. At the microscopic level, the flow is computed by solving the full compressible Navier–Stokes equations in two dimensions. Special attention is given to the outlet boundary conditions. The analysis is based on the comparison between the macroscopic data, obtained on the one hand by spatially averaging the microscopic results, and on the other hand by solving the problem directly at the macroscopic level. Situations for which a good agreement is found between the two series of data and situations for which discrepancies are observed are exhibited. These various behaviours are discussed in terms of the various time scales controlling the flow and are explained by analysing the flow structure at pore level. Copyright © 1999 John Wiley & Sons, Ltd.     

An iterative method for the computation of the response of linearised Navier–Stokes equations to harmonic forcing and application to forced cylinder wakes

The central aim of this paper is the development and application of an efficient, iterative methodology for the computation of the perturbation fields induced by harmonic forcing of the linearised Navier–Stokes equations. The problem is formulated directly in the frequency domain, and the resulting system of equations is solved iteratively until convergence. The method is easily implemented to any implicit code that can solve iteratively the steady‐state Navier–Stokes equations. In this paper, it is applied to investigate the flow around a static cylinder with pulsating approaching flow and a cylinder undergoing forced stream‐wise oscillations. All terms of the perturbation kinetic energy equation are computed, and it is shown that perturbations grow by extracting energy from two sources: the underlying base flow field and the externally provided energy that maintains the imposed oscillation. The periodic drag force acting on the cylinder is also computed, and it is demonstrated that Morrison's equation is a simple model that can estimate with good accuracy the amplitude and phase of this force with respect to the approaching flow. The perturbation fields induced by periodic inlet flow (static cylinder) and forced stream‐wise cylinder oscillation are closely related: the velocity fields are identical in the appropriate reference frames, and a simple expression is derived, which links the pressures in the two flow cases. Copyright © 2014 John Wiley & Sons, Ltd.     

Comment on ‘generalized potential flow theory and direct calculation of velocities applied to the numerical solution of the Navier–Stokes and the Boussinesq equations’

In a recent paper a generalized potential flow theory and its application to the solution of the Navier–Stokes equation are developed.¹ The purpose of this comment is to show that the analysis presented in that paper is in general not correct. We note that the theoretical development of Reference 1 is in fact an extension—although not cited—of some work first done by Hawthorne for steady inviscid flow.² Hawthorne's solution is correct, and his analysis, which we briefly describe, provides a useful introduction to this note.     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

Lattice Boltzmann modeling of pore‐scale fluid flow through idealized porous media

In order to capture the hydro‐mechanical impacts on the solid skeleton imposed by the fluid flowing through porous media at the pore‐scale, the flow in the pore space has to be modeled at a resolution finer than the pores, and the no‐slip condition needs to be enforced at the grain–fluid interface. In this paper, the lattice Boltzmann method (LBM), a mesoscopic Navier–Stokes solver, is shown to be an appropriate pore‐scale fluid flow model. The accuracy and lattice sensitivity of LBM as a fluid dynamics solver is demonstrated in the Poiseuille channel flow problem (2‐D) and duct flow problem (3‐D). Well‐studied problems of fluid creeping through idealized 2‐D and 3‐D porous media (J. Fluid Mech. 1959; 5 (2):317–328, J. Fluid Mech. 1982; 115 :13–26, Int. J. Multiphase Flow 1982; 8 (4):343–360, Phys. Fluids A 1989; 1 (1):38–46, Int. J. Numer. Anal. Meth. Geomech. 1999; 23 :881–904, Int. J. Numer. Anal. Meth. Geomech. 2010; DOI: 10.1002/nag.898, Int. J. Multiphase Flow 1982; 8 (3):193–206) are then simulated using LBM to measure the friction coefficient for various pore throats. The simulation results agree well with the data reported in the literature. The lattice sensitivity of the frictional coefficient is also investigated. Copyright © 2010 John Wiley & Sons, Ltd.     

On the immersed boundary‐lattice Boltzmann simulations of incompressible flows with freely moving objects

For simulating freely moving problems, conventional immersed boundary‐lattice Boltzmann methods encounter two major difficulties of an extremely large flow domain and the incompressible limit. To remove these two difficulties, this work proposes an immersed boundary‐lattice Boltzmann flux solver (IB‐LBFS) in the arbitrary Lagragian–Eulerian (ALE) coordinates and establishes a dynamic similarity theory. In the ALE‐based IB‐LBFS, the flow filed is obtained by using the LBFS on a moving Cartesian mesh, and the no‐slip boundary condition is implemented by using the boundary condition‐enforced immersed boundary method. The velocity of the Cartesian mesh is set the same as the translational velocity of the freely moving object so that there is no relative motion between the plate center and the mesh. This enables the ALE‐based IB‐LBFS to study flows with a freely moving object in a large open flow domain. By normalizing the governing equations for the flow domain and the motion of rigid body, six non‐dimensional parameters are derived and maintained to be the same in both physical systems and the lattice Boltzmann framework. This similarity algorithm enables the lattice Boltzmann equation‐based solver to study a general freely moving problem within the incompressible limit. The proposed solver and dynamic similarity theory have been successfully validated by simulating the flow around an in‐line oscillating cylinder, single particle sedimentation, and flows with a freely falling plate. The obtained results agree well with both numerical and experimental data. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of multifluid flows with large time steps using the particle finite element method

Multifluids are those fluids in which their physical properties (viscosity or density) vary internally and abruptly forming internal interfaces that introduce a large nonlinearity in the Navier–Stokes equations. For this reason, standard numerical methods require very small time steps in order to solve accurately the internal interface position. In a previous paper, the authors developed a particle‐based method (named particle finite element method (PFEM)) based on a Lagrangian formulation and FEM for solving the fluid mechanics equations for multifluids. PFEM was capable of achieving accurate results, but the limitation of small time steps was still present. In this work, a new strategy concerning the time integration for the analysis of multifluids is developed allowing time steps one order of magnitude larger than the previous method. The advantage of using a Lagrangian solution with PFEM is shown in several examples. All kind of heterogeneous fluids (with different densities or viscosities), multiphase flows with internal interfaces, breaking waves, and fluid separation may be easily solved with this methodology without the need of small time steps. Copyright © 2014 John Wiley & Sons, Ltd.     

Linear and non‐linear simulations of feedback control in plane Poiseuille flow

This paper examines the performance of optimal linear quadratic state and output feedback controllers in stabilizing two‐dimensional perturbations in a plane Poiseuille flow. The synthesis of the controllers is based on a linearized model of the flow using a new set of interpolating polynomials in the wall‐normal direction, which automatically satisfy the homogeneous Dirichlet and Neumann boundary conditions at the walls and eliminate spurious eigenvalues. The controllers are implemented into a non‐linear Navier–Stokes solver, which is modified to compute the evolution of the flow perturbations. Two cases are examined, one with small initial disturbances that do not violate the linearity assumptions and the other with much larger disturbances that trigger the non‐linear convection terms. For the smallest disturbances, the solver accurately reproduced the results of the linear simulations of open‐ and closed‐loop systems. The simulations for the larger disturbances without control showed a rapid initial growth but the flow soon reached a saturated state in agreement with previous findings in the literature. The large initial growth is a consequence of the non‐normal nature of the system dynamics. The state feedback and output feedback controllers were able to reduce significantly the perturbation energy. For the larger disturbances, the energy calculated from the state variables is well below the energy evaluated by direct integration of the velocity field. This is probably due to the non‐linear terms transferring energy to harmonics of the considered wavenumber, which are not sensed by the linear controller. Copyright © 2008 John Wiley & Sons, Ltd.     

Numerical simulations of laminar flow over a 3D backward-facing step

A numerical investigation of laminar flow over a three-dimensional backward-facing step is presented with comparisons with detailed experimental data, available in the literature, serving to validate the numerical results. The continuity constraint method, implemented via a finite element weak statement, was employed to solve the unsteady three-dimensional Navier–Stokes equations for incompressible laminar isothermal flow. Two-dimensional numerical simulations of this step geometry underestimate the experimentally determined extent of the primary separation region for Reynolds numbers Re greater than 400. It has been postulated that this disagreement between physical and computational experiments is due to the onset of three-dimensional flow near Re ≈ 400. This paper presents a full three-dimensional simulation of the step geometry for 100⩽ Re⩽ 800 and correctly predicts the primary reattachment lengths, thus confirming the influence of three-dimensionality. Previous numerical studies have discussed possible instability modes which could induce a sudden onset of three-dimensional flow at certain critical Reynolds numbers. The current study explores the influence of the sidewall on the development of three-dimensional flow for Re greater than 400. Of particular interest is the characterization of three-dimensional vortices in the primary separation region immediately downstream of the step. The complex interaction of a wall jet, located at the step plane near the sidewall, with the mainstream flow reveals a mechanism for the increasing penetration (with increasing Reynolds number) of three-dimensional flow structures into a region of essentially two-dimensional flow near the midplane of the channel. The character and extent of the sidewall-induced flow are investigated for 100⩽Re⩽ 800. © 1997 John Wiley & Sons, Ltd.     

Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.     

Multigrid solution of the Navier–Stokes equations on highly stretched grids

Relaxation-based multigrid solvers for the steady incompressible Navier–Stokes equations are examined to determine their computational speed and robustness. Four relaxation methods were used as smoothers in a common tailored multigrid procedure. The resulting solvers were applied to three two-dimensional flow problems, over a range of Reynolds numbers, on both uniform and highly stretched grids. In all cases the L₂ norm of the velocity changes is reduced to 10^?6 in a few 10's of fine-grid sweeps. The results of the study are used to draw conciusions on the strengths and weaknesses of the individual relaxation methods as well as those of the overall multigrid procedure when used as a solver on highly stretched grids.     

Effect of vertical grid variability on a free surface flow model

A previously developed numerical model that solves the incompressible, non‐hydrostatic, Navier–Stokes equations for free surface flow is analysed on a non‐uniform vertical grid. The equations are vertically transformed to the σ‐coordinate system and solved in a fractional step manner in which the pressure is computed implicitly by correcting the hydrostatic flow field to be divergence free. Numerical consistency, accuracy and efficiency are assessed with analytical methods and numerical experiments for a varying vertical grid discretization. Specific discretizations are proposed that attain greater accuracy and minimize computational effort when compared to a uniform vertical discretization. Published in 2007 by John Wiley & Sons, Ltd.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

Determination of permeability in fibrous porous media using the lattice Boltzmann method with application to PEM fuel cells

The lattice Boltzmann method (LBM) is used to simulate the flow through an idealized proton exchange membrane fuel cell (PEMFC) porous transport layer (PTL) geometry generated using a Monte Carlo method. Using the calculated flow field, Darcy's law is applied and the permeability is calculated. This process is applied in both through‐ and in‐plane directions of the paper as both of these permeability values are important in computational fluid dynamics models of PEMFCs. It is shown that the LBM can be used to determine permeability in a random porous media by solving the flow in the microstructure of the material. The permeability in the through‐ and in‐plane directions is shown to be different and the anisotropic nature of the geometry creates anisotropic permeability. It is also found that fiber arrangement plays a large role in the permeability of the PTL. New correlations are presented for in‐ and though‐plane permeabilities of fibrous porous media with (0.6<ε<0.8). Copyright © 2008 John Wiley & Sons, Ltd.     

Fast solvers for models of ICEO microfluidic flows

We demonstrate the performance of a fast computational algorithm for modeling the design of a microfluidic mixing device. The device uses an electrokinetic process, induced charge electroosmosis (J. Fluid Mech. 2004; 509 ), by which a flow through the device is driven by a set of polarizable obstacles in it. Its design is realized by manipulating the shape and orientation of the obstacles in order to maximize the amount of fluid mixing within the device. The computation entails the solution of a constrained optimization problem in which function evaluations require the numerical solution of a set of partial differential equations: a potential equation, the incompressible Navier–Stokes equations, and a mass‐transport equation. The most expensive component of the function evaluation (which must be performed at every step of an iteration for the optimization) is the solution of the Navier–Stokes equations. We show that by using some new robust algorithms for this task (SIAM J. Sci. Comput. 2002; 24 :237–256; J. Comput. Appl. Math. 2001; 128 :261–279), based on certain preconditioners that take advantage of the structure of the linearized problem, this computation can be done efficiently. Using this computational strategy, in conjunction with a derivative‐free pattern search algorithm for the optimization, applied to a finite element discretization of the problem, we are able to determine optimal configurations of microfluidic devices. Copyright © 2009 John Wiley & Sons, Ltd.     

On numerical errors and turbulence modeling in tip vortex flow prediction

The accuracy of tip vortex flow prediction in the near‐field region is investigated numerically by attempting to quantify the shortcomings of the turbulence models and the flow solver. In particular, some turbulence models can produce a ‘numerical diffusion’ that artificially smears the vortex core. Low‐order finite differencing techniques of the convective and pressure terms of the Navier–Stokes equations and inadequate grid density and distribution can also produce the same adverse effect. The flow over a wing and the near‐wake with the wind tunnel walls included was simulated using 2.5 million grid points. Two subset problems, one using a steady, three‐dimensional analytical vortex, and the other, a vortex obtained from experiment and propagated downstream, were also devised in order to make the study of vortex preservation more tractable. The method of artificial compressibility is used to solve the steady, three‐dimensional, incompressible Navier–Stokes equations. Two one‐equation turbulence models (Baldwin–Barth and Spalart–Allmaras turbulence models), have been used with the production term modified to account for the stabilizing effect of the nearly solid body rotation in the vortex core. Finally, a comparison between the computed results and experiment is presented. Published in 1999 by John Wiley & Sons, Ltd.