Development of projection and artificial compressibility methodologies using the approximate factorization technique

Predictions for two-dimensional, steady, incompressible flows under both laminar and turbulent conditions are presented. The standard k-? turbulence model is used for the turbulent flows. The computational method is based on the approximate factorization technique. The coupled approach is used to link the equations of motion and the turbulence model equations. Mass conservation is enforced by either the pseudocompressibility method or the pressure correction method. Comparison of the two methods shows a superiority of the pressure correction method. Second- and fourth-order artifical dissipation terms are used in order to achieve good convergence and to handle the turbulence model equations efficiently. Several internal and external test cases are investigated, including attached and separated flows.     

Numerical solutions of the shallow water equations with discontinuous bed topography

A simple scheme is developed for treatment of vertical bed topography in shallow water flows. The effect of the vertical step on flows is modelled with the shallow water equations including local energy loss terms. The bed elevation is denoted with z_b^‐ for the left and z_b⁺ for the right values at each grid point, hence exactly representing a discontinuity in the bed topography. The surface gradient method (SGM) is generalized to reconstruct water depths at cell interfaces involving a vertical step so that the fluxes at the cell interfaces can accurately be calculated with a Riemann solver. The scheme is verified by predicting a surge crossing a step, a tidal flow over a step and dam‐break flows on wet/dry beds. The results have shown good agreements compared with analytical solutions and available experimental data. The scheme is efficient, robust, and may be used for practical flow calculations. Copyright © 2002 John Wiley & Sons, Ltd.     

DSMC／EPSM混合算法研究

将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC/EPSM混合算法。运用混合算法模拟了马赫数等于5时超音速竖板绕流及马赫数等于4时超音速平板绕流，并将结果与DSMC算法的结果进行比较，证明了DSMC／EPSM混合算法的有效性，同时将EPSM算法与DSMC算法的效率进行了比较。     

Optimization of unsteady incompressible Navier–Stokes flows using variational level set method

This paper presents the optimization of unsteady Navier–Stokes flows using the variational level set method. The solid–liquid interface is expressed by the level set function implicitly, and the fluid velocity is constrained to be zero in the solid domain. An optimization problem, which is constrained by the Navier–Stokes equations and a fluid volume constraint, is analyzed by the Lagrangian multiplier based adjoint approach. The corresponding continuous adjoint equations and the shape sensitivity are derived. The level set function is evolved by solving the Hamilton–Jacobian equation with the upwind finite difference method. The optimization method can be used to design channels for flows with or without body forces. The numerical examples demonstrate the feasibility and robustness of this optimization method for unsteady Navier–Stokes flows.Copyright © 2012 John Wiley & Sons, Ltd.     

Prediction of turbulent wall shear flows directly from wall

The fully elliptic Reynolds-averaged Navier–Stokes equations have been used together with Lam and Bremhorst's low-Reynolds-number model, Chen and Patel's two-layer model and a two-point wall function method incorporated into the standard k-? model to predict channel flows and a backward-facig step flow. These flows enable the evaluation of the performance of different near-wall treatments in flows involving streamwise and normal pressure gradients, flows with separation and flows with non-equilibrium turbulence characteristics. Direct numerical simulation (DNS) of a channel flow with Re =3200 further provides the detailed budgets of each modelling term of the k and ?-transport equations. Comparison of model results with DNS data to evaluate the performance of each modelling term is also made in the present study. It is concluded that the low-Reynolds-number model has wider applicability and performs better than the two-layer model and wall function approaches. Comparison with DNS data further shows that large discrepancies exist between the DNS budgets and the modelled production and destruction terms of the ? equation. However, for simple channel flow the discrepancies are similar in magnitude but opposite in sign, so they are cancelled by each other. This may explain why, even when employing such an inaccurately modelled ?-equation, one can still predict satisfactorily some simple turbulent flows.     

NUMERICAL STUDY OF PERIODICAL FLOWS OF PIEZOELECTRIC VALVELESS MICROPUMP FOR BIOCHIPS

IntroductionMicro-pump is one of the most important fluid-driving elements in biochips. Workingprinciple of micropumps is quite different from traditional macro-scale pumps. Currently,valveless piezoelectric micropump is one of the most popular micropumps…     

Concurrent use of finite element and finite volume methods for shallow water flows in locally 1‐D channel networks

A numerical model is developed for shallow water equation in locally 1‐D channel networks. The model concurrently uses the standard Galerkin finite element method for the continuity equation and the finite volume method with an upwind scheme for the momentum equation. The surface gradient method is consistently employed. A minimum treatment is given for channel junctions so that application to multiply connected channels do not require any special consideration The model is capable of computing different types of transcritical flows, wet and dry flows, and flows with complex source terms. Standardized test problems and laboratory experimental data are used for verifying the model. Applicability of the models is validated in a multiply connected channel network draining hydromorphic farmlands located in a West African savanna, and Manning's roughness coefficient is identified, so that the steady solution is consistent with field observations. Unsteady simulation demonstrates that the model is capable of stably reproducing shifts of hydraulic jumps in flows of sub‐millimeter water depths. Copyright © 2011 John Wiley & Sons, Ltd.     

A lattice Boltzmann model for simulation of compressible flows

This paper presents a new model of lattice Boltzmann method for full compressible flows. On the basis of multi‐speed model, an extra potential energy distribution function is introduced to recover the full compressible Navier–Stokes equations with a flexible specific‐heat ratio and Prandtl number. The Chapman–Enskog expansion of the kinetic equations is performed, and the two‐dimension‐seventeen‐velocity density equilibrium distribution functions are obtained. The governing equations are discretized using the third order monotone upwind scheme for scalar conservation laws finite volume scheme. The van Albada limiter is used to avoid spurious oscillations. In order to verify the accuracy of this double‐distribution‐function model, the Riemann problems, Couette flows, and flows around a NACA0012 airfoil are simulated. It is found that the proposed lattice Boltzmann model is suitable for compressible flows, even for strong shock wave problem, which has an extremely large pressure ratio, 100,000. Copyright © 2014 John Wiley & Sons, Ltd.     

Simulation of 2D external viscous flows by means of a domain decomposition method using an influence matrix technique

Two-dimensional external viscous flows are numerically approximated by means of a domain decomposition method which combines a vortex method and a finite difference method. The vortex method is used in the flow region which is dominated by convective effects, whereas the finite difference method is used in the flow region where viscous diffusion effects are dominant. An influence matrix technique combined with the uniformity condition of the pressure is used to enforce the tangential velocity boundary condition. Comparisons between numerical and experimental data show that the method is well adapted for simulating two-dimensional flows.     

Simulation of incompressible viscous flows around moving objects by a variant of immersed boundary‐lattice Boltzmann method

A variant of immersed boundary‐lattice Boltzmann method (IB‐LBM) is presented in this paper to simulate incompressible viscous flows around moving objects. As compared with the conventional IB‐LBM where the force density is computed explicitly by Hook's law or the direct forcing method and the non‐slip condition is only approximately satisfied, in the present work, the force density term is considered as the velocity correction which is determined by enforcing the non‐slip condition at the boundary. The lift and drag forces on the moving object can be easily calculated via the velocity correction on the boundary points. The capability of the present method for moving objects is well demonstrated through its application to simulate flows around a moving circular cylinder, a rotationally oscillating cylinder, and an elliptic flapping wing. Furthermore, the simulation of flows around a flapping flexible airfoil is carried out to exhibit the ability of the present method for implementing the elastic boundary condition. It was found that under certain conditions, the flapping flexible airfoil can generate larger propulsive force than the flapping rigid airfoil. Copyright © 2009 John Wiley & Sons, Ltd.     

Coupling finite difference method with finite particle method for modeling viscous incompressible flows

A hybrid approach to couple finite difference method (FDM) with finite particle method (FPM) (ie, FDM-FPM) is developed to simulate viscous incompressible flows. FDM is a grid-based method that is convenient for implementing multiple or adaptive resolutions and is computationally efficient. FPM is an improved smoothed particle hydrodynamics (SPH), which is widely used in modeling fluid flows with free surfaces and complex boundaries. The proposed FDM-FPM leverages their advantages and is appealing in modeling viscous incompressible flows to balance accuracy and efficiency. In order to exchange the interface information between FDM and FPM for achieving consistency, stability, and convergence, a transition region is created in the particle region to maintain the stability of the interface between two methods. The mass flux algorithm is defined to control the particle creation and deletion. The mass is updated by N-S equations instead of the interpolation. In order to allow information exchange, an overlapping zone is defined near the interface. The information of overlapping zone is obtained by an FPM-type interpolation. Taylor-Green vortices and lid-driven shear cavity flows are simulated to test the accuracy and the conservation of the FDM-FPM hybrid approach. The standing waves and flows around NACA airfoils are further simulated to test the ability to deal with free surfaces and complex boundaries. The results show that FDM-FPM retains not only the high efficiency of FDM with multiple resolutions but also the ability of FPM in modeling free surfaces and complex boundaries.     

A FINITE VOLUME SOLVER FOR THE SHALLOW WATER EQUATIONS IN VARYING BED ENVIRONMENTS

Abstract

A numerical scheme for solving the shallow-water equations is presented. An analogy is made between flows governed by shallow-water equations and the Euler system of equations used in gas dynamics. An emphasis is placed on the difference presented by the bathymetry in hydraulic systems. The discretization of the governing equations is based on Roe's flux difference-splitting solver, initially developed for solving inviscid compressible flows. The spatial discretization is handled within a finite-volume context by using triangles or quadrilaterals as the basic control-volume cells. This approach enables an easy and flexible treatment of general geometries. A development of the boundary conditions tailored for the current scheme is given. Fundamental validation tests are presented.     

Simulation of supersonic and hypersonic flows

A simple convection algorithm for simulation of time-dependent supersonic and hypersonic flows of a perfect but viscous gas is described. The algorithm is based on conservation and convection of mass, momentum and energy in a grid of rectangular cells. Examples are given for starting flow in a shock tube and oblique shocks generated by a wedge at Mach numbers up to 30·4. Good comparisons are achieved with well-known perfect gas flows.     

Effects of slip velocity and surface acoustic wave on the laminar flow stability

The stability of slip flows when a surface acoustic wave (SAW) propagating along the walls of a microchannel in the laminar flow regime is investigated. The governing equation which was derived by considering the weakly nonlinear coupling between the deformable wall and streaming slip flow is linearized and then the eigenvalue problem is solved by a numerical code together with the associated interface and slip velocity boundary conditions. The value of the critical Reynolds number was found to be near 1,441 for a Knudsen number being 0.001 (associated with a physical parameter K ₀ characterizing the SAW effect) which is much smaller than the static-wall case for conventional pressure-driven flows.     

Simulation of three‐dimensional flows over moving objects by an improved immersed boundary–lattice Boltzmann method

An improved immersed boundary–lattice Boltzmann method (IB–LBM) developed recently [28] was applied in this work to simulate three‐dimensional (3D) flows over moving objects. By enforcing the non‐slip boundary condition, the method could avoid any flow penetration to the wall. In the developed IB–LBM solver, the flow field is obtained on the non‐uniform mesh by the efficient LBM that is based on the second‐order one‐dimensional interpolation. As a consequence, its coefficients could be computed simply. By simulating flows over a stationary sphere and torus [28] accurately and efficiently, the proposed IB–LBM showed its ability to handle 3D flow problems with curved boundaries. In this paper, we further applied this method to simulate 3D flows around moving boundaries. As a first example, the flow over a rotating sphere was simulated. The obtained results agreed very well with the previous data in the literature. Then, simulation of flow over a rotating torus was conducted. The capability of the improved IB–LBM for solving 3D flows over moving objects with complex geometries was demonstrated via the simulations of fish swimming and dragonfly flight. The numerical results displayed quantitative and qualitative agreement with the date in the literature. Copyright © 2011 John Wiley & Sons, Ltd.     

An immersed‐boundary finite‐volume method for direct simulation of flows with suspended paramagnetic particles

In the present study, we have proposed an immersed‐boundary finite‐volume method for the direct numerical simulation of flows with inertialess paramagnetic particles suspended in a nonmagnetic fluid under an external magnetic field without the need for any model such as the dipole–dipole interaction. In the proposed method, the magnetic field (or force) is described by the numerical solution of the Maxwell equation without current, where the smoothed representation technique is employed to tackle the discontinuity of magnetic permeability across the particle–fluid interface. The flow field, on the other hand, is described by the solution of the continuity and momentum equations, where the discrete‐forcing‐based immersed‐boundary method is employed to satisfy the no‐slip condition at the interface. To validate the method, we performed numerical simulations on the two‐dimensional motion of two and three paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and then compared the results with the existing finite‐element and semi‐analytical solutions. Comparison shows that the proposed method is robust in the direct simulation of such magnetic particulate flows. This method can be extended to more general flows without difficulty: three‐dimensional particulate flows, flows with a great number of particles, or flows under an arbitrary external magnetic field. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulation of Three-Dimensional Bluff-Body Flows Using the Vortex Particle and Boundary Element Methods

Recent contributions to the 3-D vortex method for bluff-body flows are presented. The numerical method--a vortex method combined with a boundary element method--is briefy reviewed. It is applied to direct numerical simulation (DNS) of the flow past a sphere (Re= 300, 500 and 1000). The on-going work to extend the method towards vortex-based large-eddy simulation (LES) for high Reynolds number flows is also presented. Preliminary results for the flow past a hemisphere are discussed.     

Parallel large eddy simulation of turbulent flow around MIRA model using linear equal‐order finite element method

A parallel large eddy simulation code that adopts domain decomposition method has been developed for large‐scale computation of turbulent flows around an arbitrarily shaped body. For the temporal integration of the unsteady incompressible Navier–Stokes equation, fractional 4‐step splitting algorithm is adopted, and for the modelling of small eddies in turbulent flows, the Smagorinsky model is used. For the parallelization of the code, METIS and Message Passing Interface Libraries are used, respectively, to partition the computational domain and to communicate data between processors. To validate the parallel architecture and to estimate its performance, a three‐dimensional laminar driven cavity flow inside a cubical enclosure has been solved. To validate the turbulence calculation, the turbulent channel flows at Re_τ = 180 and 1050 are simulated and compared with previous results. Then, a backward facing step flow is solved and compared with a DNS result for overall code validation. Finally, the turbulent flow around MIRA model at Re = 2.6 × 10⁶ is simulated by using approximately 6.7 million nodes. Scalability curve obtained from this simulation shows that scalable results are obtained. The calculated drag coefficient agrees better with the experimental result than those previously obtained by using two‐equation turbulence models. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical solutions of incompressible Euler and Navier-Stokes equations by efficient discrete singular convolution method

An efficient discrete singular convolution (DSC) method is introduced to the numerical solutions of incompressible Euler and Navier-Stokes equations with periodic boundary conditions. Two numerical tests of two-dimensional Navier-Stokes equations with periodic boundary conditions and Euler equations for doubly periodic shear layer flows are carried out by using the DSC method for spatial derivatives and fourth-order Runge-Kutta method for time advancement, respectively. The computational results show that the DSC method is efficient and robust for solving the problems of incompressible flows, and has the potential of being extended to numerically solve much broader problems in fluid dynamics. The project supported by the National Natural Science Foundation of China (No.19902010).     

微尺度气体流动

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和Millikan的小球阻力系数的测量,实验结果揭示了稀薄气体效应即尺度效应对气体运动的重要影响.由于流动特征长度很小,微尺度气流经常处于滑流区甚至过渡领域,流动的相似参数为Knudsen数和Mach数.因此可以考虑利用相似准则,通过增大几何尺寸、减小压力的途径,解决微机电系统实验观测遇到的困难.为解决直接模拟Monte Carlo方法分析微机电系统中低速稀薄气流遇到的统计涨落困难,我们提出了信息保存法(IP),该方法能够有效克服统计散布,并已成功用于多种微尺度气流.