An improved near‐wall modeling for large‐eddy simulation using immersed boundary methods

An improved near‐wall modeling for large‐eddy simulation using the immersed boundary method is proposed. It is shown in this study that the existing near‐wall modeling for the immersed boundary (IB) methods that imposes the velocity boundary condition at the IB node is not sufficient to enforce a correct wall shear stress at the IB node. A new method that imposes a shear stress condition through the modification of the subgrid scale‐eddy viscosity at the IB node is proposed. In this method, the subgrid eddy viscosity at the IB node is modified such that the viscous flux at the face adjacent to the IB node correctly approximates the total shear stress. The method is applied to simulate the fully developed turbulent flows in a plane channel and a circular pipe. It is demonstrated that the new method improves the prediction of the mean velocity and turbulence stresses in comparison with the existing wall modeling based solely on the velocity boundary condition. Copyright © 2015 John Wiley & Sons, Ltd.     

适用于浸入边界法的大涡模拟紊流壁面模型

基于近壁定常剪切应力假设,提出了一种新的适用于浸入边界法的大涡模拟紊流壁面模型。通过引入壁面滑移速度,修正了线性速度剖面计算得到的壁面剪切应力,使之满足Werner-Wengle模型。将其应用于平板紊流和高Re数圆管紊流的数值模拟,对比采用和不采用壁面模型的结果得知,采用此模型的速度剖面与实验值吻合良好,验证了此模型的有效性。研究了不同欧拉/拉格朗日网格相对位置对结果的影响,证明了此模型具有较好的鲁棒性,以及可根据局部流动状态和网格精度自动开闭的特点。     

A Lagrangian method for the treatment of freshly cleared cells in immersed boundary techniques

The present study describes a novel approach treating freshly cleared cells when simulating incompressible flows using the immersed boundary method. The method calculates material derivative in Lagrangian coordinate system to approximate the solution at the freshly cleared cells. In this study, we use a simple test case to compare our proposed method with the commonly used one-dimensional interpolation method in the Eulerian coordinate system.     

An immersed boundary finite difference method for LES of flow around bluff shapes

A three‐dimensional numerical model using large eddy simulation (LES) technique and incorporating the immersed boundary (IMB) concept has been developed to compute flow around bluff shapes. A fractional step finite differences method with rectilinear non‐uniform collocated grid is employed to solve the governing equations. Bluff shapes are treated in the IMB method by introducing artificial force terms into the momentum equations. Second‐order accurate interpolation schemes for all sorts of grid points adjacent to the immersed boundary have been developed to determine the velocities and pressure at these points. To enforce continuity, the methods of imposition of pressure boundary condition and addition of mass source/sink terms are tested. It has been found that imposing suitable pressure boundary condition (zero normal gradient) can effectively reproduce the correct pressure distribution and enforce mass conservation around a bluff shape. The present model has been verified and applied to simulate flow around bluff shapes: (1) a square cylinder and (2) the Tsing Ma suspension bridge deck section model. Complex flow phenomena such as flow separation and vortex shedding are reproduced and the drag coefficient, lift coefficient, and pressure coefficient are calculated and analyzed. Good agreement between the numerical results and the experimental data are obtained. The model is proven to be an efficient tool for flow simulation around bluff bodies in time varying flows. Copyright © 2004 John Wiley & Sons, Ltd.     

A ghost-cell immersed boundary method for large eddy simulation of flows in complex geometries

An efficient ghost-cell immersed boundary (IB) method is proposed for large eddy simulations of three-dimensional incompressible flow in complex geometries. In the framework of finite volume method, the Navier–Stokes equations are integrated using an explicit time advancement scheme on a collocated mesh. Since the IB method is known to generate an unphysical velocity field inside the IB that violates the mass conservation of the cells near the IB, a new IB treatment is devised to eliminate the unphysical velocity generated near the IB and to improve the pressure distribution on the body surface. To validate the proposed method, both laminar and turbulent flow cases are presented. In particular, large eddy simulations were performed to simulate the turbulent flows over a circular cylinder and a sphere at subcritical Reynolds numbers. The computed results show good agreements with the published numerical and experimental data.     

An immersed boundary/wall modeling method for RANS simulation of compressible turbulent flows

In this paper, an immersed boundary (IB) method is developed to simulate compressible turbulent flows governed by the Reynolds‐averaged Navier‐Stokes equations. The flow variables at the IB nodes (interior nodes in the immediate vicinity of the solid wall) are evaluated via linear interpolation in the normal direction to close the discrete form of the governing equations. An adaptive wall function and a 2‐layer wall model are introduced to reduce the near‐wall mesh density required by the high resolution of the turbulent boundary layers. The wall shear stress modified by the wall modeling technique and the no‐penetration condition are enforced to evaluate the velocity at an IB node. The pressure and temperature at an IB node are obtained via the local simplified momentum equation and the Crocco‐Busemann relation, respectively. The SST k ? ω and S‐A turbulence models are adopted in the framework of the present IB approach. For the Shear‐Stress Transport (SST) k ? ω model, analytical solutions in near‐wall region are utilized to enforce the boundary conditions of the turbulence equations and evaluate the turbulence variables at an IB node. For the S‐A model, the turbulence variable at an IB node is calculated by using the near‐wall profile of the eddy viscosity. In order to validate the present IB approach, numerical experiments for compressible turbulent flows over stationary and moving bodies have been performed. The predictions show good agreements with the referenced experimental data and numerical results.     

基于反馈力浸入边界法模拟复杂动边界流动

浸入边界法是模拟流固耦合的重要数值方法之一。本文采用反馈力浸入边界方法,对旋转圆柱和水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟。其中,固体边界采用一系列离散的点近似代替,流体为不可压缩牛顿流体,使用笛卡尔自适应加密网格,利用有限差分法进行求解。固体对流场的作用通过构造适宜的反馈力函数实现。本文首先通过旋转圆柱绕流的计算结果同实验结果进行对比,吻合较好,验证了该计算方法的可靠性。然后针对水电站水力过渡过程中水轮机活动导叶旋转摆动绕流后的动边界流场进行数值模拟,得到导叶动态绕流后的流场分布特性和涡结构的演化特性。     

Moving immersed boundary method

A novel implicit immersed boundary method of high accuracy and efficiency is presented for the simulation of incompressible viscous flow over complex stationary or moving solid boundaries. A boundary force is often introduced in many immersed boundary methods to mimic the presence of solid boundary, such that the overall simulation can be performed on a simple Cartesian grid. The current method inherits this idea and considers the boundary force as a Lagrange multiplier to enforce the no‐slip constraint at the solid boundary, instead of applying constitutional relations for rigid bodies. Hence excessive constraint on the time step is circumvented, and the time step only depends on the discretization of fluid Navier‐Stokes equations, like the CFL condition in present work. To determine the boundary force, an additional moving force equation is derived. The dimension of this derived system is proportional to the number of Lagrangian points describing the solid boundaries, which makes the method very suitable for moving boundary problems since the time for matrix update and system solving is not significant. The force coefficient matrix is made symmetric and positive definite so that the conjugate gradient method can solve the system quickly. The proposed immersed boundary method is incorporated into the fluid solver with a second‐order accurate projection method as a plug‐in. The overall scheme is handled under an efficient fractional step framework, namely, prediction, forcing, and projection. Various simulations are performed to validate current method, and the results compare well with previous experimental and numerical studies.     

Large eddy simulation in a turbulent channel flow using exit boundary conditions

The influence of the exit boundary conditions on the vanishing first derivative of the velocity components and constant pressure on the large eddy simulation of the fully developed turbulent channel flow has been investigated for equidistant and stretched grids at the channel exit. Results show that the chosen exit boundary conditions introduce some small disturbances that are mostly damped by the grid stretching. The difference of rms values between the fully developed turbulent channel flow with periodicity conditions and the fully developed channel flow using inlet and the exit boundary conditions is less than 10% for the equidistant grids and less than 5% for the stretched grids. The chosen boundary conditions are of interest because they may be used in complex problems with back flow at the exit. Copyright © 1999 John Wiley & Sons, Ltd.     

An immersed boundary method wall model for high‐Reynolds‐number channel flow over complex topography

High‐Reynolds‐number channel flows regularly encounter topographies composed of multiple length scales and that protrude into the boundary layer. Physically, the presence of immersed obstacles leads to increased velocity gradients, turbulence production, and manifestation of wakes. Considerable challenges are associated with numerically describing the presence of obstacles in channel flows. Common approaches include generation of a computational mesh that is uniquely designed for the flow and obstacle, the immersed boundary method, and terrain‐following coordinates. There are challenges and limitations associated with each of these techniques. Specification of boundary conditions representing the perimeter of solid obstacles is a primary challenge of the immersed boundary method. In this document, a simplistic canopy stress‐like wall model is used to impose boundary conditions. The model isolates aerodynamically relevant local frontal areas through evaluation of the gradient of the topographic height field. The gradient of the height field describes both the surface‐normal direction and the frontal area, making it ideal for detecting areas on which the flow impinges. The model is tested in numerical simulations of turbulent half‐channel flow over topographies with different obstacles affixed–right prisms, rectangular prisms, ellipsoidal mounds, and sinusoids. In all cases, the performance is strong relative to datasets presented in the literature. Results are finally presented for numerical simulation of flow over complex synthetic fractal‐like topography and a synthetic city. These results show interesting trends in how the turbulent multiscale flow field responds to multiscale topography. Copyright © 2012 John Wiley & Sons, Ltd.     

A DNS study of jet control with microjets using an immersed boundary method

In this work, a microjet arrangement to control a turbulent jet is studied by means of direct numerical simulation. A customised numerical strategy was developed to investigate the interactions between the microjets and the turbulent jet. This approach is based on an improved immersed boundary method in order to reproduce realistically the control device while being compatible with the accuracy and the parallel strategy of the in-house code Incompact3d. The 16 converging microjets, so-called fluidevrons, lead to an increase of the turbulent kinetic energy in the near-nozzle region through an excitation at small scale caused by the interaction between the fluidevrons and the main jet. As a consequence, very intense unstable ejections are produced from the centre of the jet toward its surrounding. Further downstream, the turbulent kinetic energy levels are lower with a lengthening of the potential core compared to a natural jet, in agreement with experimental results.     

An immersed boundary method for simulation of inviscid compressible flows

In this paper, an immersed boundary method for simulating inviscid compressible flows governed by Euler equations is presented. All the mesh points are classified as interior computed points, immersed boundary points (interior points closest to the solid boundary), and exterior points that are blanked out of computation. The flow variables at an immersed boundary point are determined via the approximate form of solution in the direction normal to the wall boundary. The normal velocity is evaluated by applying the no‐penetration boundary condition, and therefore, the influence of solid wall in the inviscid flow is taken into account. The pressure is computed with the local simplified momentum equation, and the density and the tangential velocity are evaluated by using the constant‐entropy relation and the constant‐total‐enthalpy relation, respectively. With a local coordinate system, the present method has been extended easily to the three‐dimensional case. The present work is the first endeavor to extend the idea of hybrid Cartesian/immersed boundary approach to compressible inviscid flows. The tedious task of handling multi‐valued points can be eliminated, and the overshoot resulting from the extrapolation for the evaluation of flow variables at exterior points can also be avoided. In order to validate the present method, inviscid compressible flows over fixed and moving bodies have been simulated. All the obtained numerical results show good agreement with available data in the literature. Copyright © 2014 John Wiley & Sons, Ltd.     

防风网透流风空气动力学特性大涡数值模拟研究

基于有限体积法建立不可压缩粘性流体运动的大涡模拟模型,采用Smagorinsky-Lilly亚格子模型,并引入浸入边界法(IBM)实现无滑移固壁边界条件,对雷诺数30~30000之间防风网透流风进行模拟研究。基于模拟结果,提出蝶型防风网透流风存在4个典型分区结构,流场中存在由蝶型形态引起的大尺度分层剪切流动,加强流体动能耗散。透流风在雷诺数300时发生层流至湍流的转捩,而在雷诺数增长至3000以上时,湍流充分发展,纵向流速脉动强度可达70%。防风网整体空气阻力远大于单个孔口射流阻力的线性叠加,射流间的相互作用以及大尺度的分层剪切结构大大增加流体阻力损失,这为通过优化孔口布置和网板形态来节省材料提供了科学依据。     

Numerical modeling of flow past a volumeless and thin rigid body using direct forcing immersed boundary method

The new capability has been added as the numerical method for modeling volumeless and thin rigid bodies to the direct forcing immersed boundary (DFIB) method. The DFIB approach is based on adding a virtual force to the Navier–Stokes equations of incompressible flow to account for the interaction between the fluid and structures. The volume of a solid function (VOS) identifies the stationary or moving solid structures in a given fluid domain. A new VOS-based algorithm was developed to identify thin, rigid structure boundary points in fluid flow and ensure that the fluid cannot cross through the boundary of a thin rigid structure while moving or stationary. The DFIB method was first validated in a three-dimensional (3D) turbulent flow over a circular cylinder. The large-eddy simulation simulated the turbulent flow scales. The proposed algorithm was tested using a 3D turbulent flow past a stationary and rotating Savonius wind turbine that functions as a thin, rigid body. The validation results showed that the selected DFIB approach, combined with the novel algorithm, could simulate a thin, volumeless, rigid structure that is stationary and rotating in incompressible turbulent flows. The current method is also applicable for two-way fluid-structure interaction problems.     

基于浸入边界法的鱼体自主游动的数值模拟

对游动或飞行生物自主运动特性的深入研究,可促进仿生学的进一步发展。本文以"C"型游动鱼作为研究对象,建立了自主游动的柔性鱼模型。此模型较为真实地反映了鱼自主游动时鱼体内力(由鱼体肌肉收缩提供)、鱼体运动和外界流体之间的耦合作用。基于传统的反馈力方法和混合有限元浸入边界法对鱼的自主游动进行了数值模拟。分析了鱼自主游动启动阶段和巡游阶段流场特性及鱼体运动特征。模拟结果表明,受到鱼体自身组织结构和外界流场作用,鱼游动时通过呈"C"型和类"S"型的不断转换,以获取能量,实现鱼体自主游动。     

A simple collision algorithm for arbitrarily shaped objects in particle-resolved flow simulation using an immersed boundary method

In the present study, we proposed a simple collision algorithm, which can be handled arbitrarily shaped objects, for flow solvers using the immersed boundary method (IBM) based on the level set and ghost cell methods. The proposed algorithm can handle the collision of the arbitrarily shaped object with little additional computational costs for the collision calculation because collision detection and calculation are performed using the level set function and image point, which are incorporated into the original IBM solver. The proposed algorithm was implemented on the solid-liquid IBM flow solver and validated by simulations of the flow over an isolated cylinder and sphere. Also, grid and time step size sensitivity on the total energy conservation of objects were investigated in cylinder-cylinder, cylinder-red-blood-cells-shaped (RBC-shaped) objects, sphere-sphere, and sphere-flat plate interaction problems. Through validation, good agreement with previous studies, grid and time step size convergence, and sufficient total energy conservation were confirmed. As a demonstration, the drafting, kissing, and tumbling processes were computed, and it was confirmed that the present result by the proposed method is similar to the previous computations. In addition, particle-laden flow in a channel including obstacles with collision and adhesion phenomena and the interaction of cylinders and wavy-wall were computed. The results of these simulations reveal the capability of solving a flow containing arbitrarily shaped moving objects with collision phenomena by a simple proposed method.     

一种投影浸入边界法的快速直接求解方法

根据投影浸入边界法分步投影求解的特点,同时针对压力泊松方程离散后的大型稀疏线性方程组是非奇异非对称的特点,结合开源函数库UMFPACK,在传递线性方程组的系数矩阵和右端向量时,采用函数库Eigen将系数矩阵的数据结构改写优化,大大降低了存储空间,实现对高维大型稀疏线性方程组的快速求解,同时求解保持良好的稳定性。本文首先利用一具有解析解的数值算例验证了求解泊松方程数值方法的准确性和网格依赖性,进而利用VC++编写投影浸入边界法的数值计算程序,以单圆柱绕流为基准数值算例,通过与其他文献和实验结果的对比,验证了投影浸入边界法数值计算结果的可靠性,并进一步分析了不同雷诺数下圆柱绕流的流场结构特征和尾涡结构的动态演化过程。     

Implicit large eddy simulation of a supersonic flat-plate boundary layer flow by weighted compact nonlinear scheme

The performance of implicit large eddy simulation (ILES) of a supersonic flat-plate turbulent boundary layer flow by weighted compact nonlinear scheme (WCNS) has been investigated. In view of features of WCNS and ILES, it was expected that ILES by WCNS could be an efficient approach to perform LES of supersonic turbulent flows. The flowfield calculated by WCNS was of lower turbulent intensity compared with an explicit LES data obtained by a numerical scheme of the same order of accuracy on a computational grid of similar resolution. It was concluded that the numerical dissipation inherent in WCNS is so large that applying WCNS to ILES of this flowfield is inefficient compared with explicit LES.