Mesh adaptation for simulation of unsteady flow with moving immersed boundaries

In this work, an approach for performing mesh adaptation in the numerical simulation of two‐dimensional unsteady flow with moving immersed boundaries is presented. In each adaptation period, the mesh is refined in the regions where the solution evolves or the moving bodies pass and is unrefined in the regions where the phenomena or the bodies deviate. The flow field and the fluid–solid interface are recomputed on the adapted mesh. The adaptation indicator is defined according to the magnitude of the vorticity in the flow field. There is no lag between the adapted mesh and the computed solution, and the adaptation frequency can be controlled to reduce the errors due to the solution transferring between the old mesh and the new one. The preservation of conservation property is mandatory in long‐time scale simulations, so a P1‐conservative interpolation is used in the solution transferring. A nonboundary‐conforming method is employed to solve the flow equations. Therefore, the moving‐boundary flows can be simulated on a fixed mesh, and there is no need to update the mesh at each time step to follow the motion or the deformation of the solid boundary. To validate the present mesh adaptation method, we have simulated several unsteady flows over a circular cylinder stationary or with forced oscillation, a single self‐propelled swimming fish, and two fish swimming in the same or different directions. Copyright © 2012 John Wiley & Sons, Ltd.     

An approach of dynamic mesh adaptation for simulating 3‐dimensional unsteady moving‐immersed‐boundary flows

In this paper, we present an approach of dynamic mesh adaptation for simulating complex 3‐dimensional incompressible moving‐boundary flows by immersed boundary methods. Tetrahedral meshes are adapted by a hierarchical refining/coarsening algorithm. Regular refinement is accomplished by dividing 1 tetrahedron into 8 subcells, and irregular refinement is only for eliminating the hanging points. Merging the 8 subcells obtained by regular refinement, the mesh is coarsened. With hierarchical refining/coarsening, mesh adaptivity can be achieved by adjusting the mesh only 1 time for each adaptation period. The level difference between 2 neighboring cells never exceeds 1, and the geometrical quality of mesh does not degrade as the level of adaptive mesh increases. A predictor‐corrector scheme is introduced to eliminate the phase lag between adapted mesh and unsteady solution. The error caused by each solution transferring from the old mesh to the new adapted one is small because most of the nodes on the 2 meshes are coincident. An immersed boundary method named local domain‐free discretization is employed to solve the flow equations. Several numerical experiments have been conducted for 3‐dimensional incompressible moving‐boundary flows. By using the present approach, the number of mesh nodes is reduced greatly while the accuracy of solution can be preserved.     

Parallel‐in‐time simulation of the unsteady Navier–Stokes equations for incompressible flow

In this paper, we present an application of a parallel‐in‐time algorithm for the solution of the unsteady Navier–Stokes model equations that are of parabolic–elliptic type. This method is based on the alternated use of a coarse global sequential solver and a fine local parallel one. A standard finite volume/finite differences first‐order approach is used for discretization of the unsteady two‐dimensional Navier–Stokes equations. The Taylor vortex decay problem and the confined flow around a square cylinder were selected as unsteady flow examples to illustrate and analyse the properties of the parallel‐in‐time method through numerical experiments. The influence of several parameters on the computing time required to perform a parallel‐in‐time calculation on a PC cluster was verified. Among them we have analysed the influence of the number of processors, the number of iterations for convergence, the resolution of the spatial domain and the influence of the time‐step sizes ratio between the coarse and fine grids. Significant computer time saving was achieved when compared with the single processor computing time, particularly when the spatial dimension of the problem is low and the temporal scale is large. Copyright © 2004 John Wiley & Sons, Ltd.     

基于动态混合网格的不可压非定常流计算方法

鱼类、昆虫等运动速度较低,对它们的数值模拟需要解决不可压问题.虚拟压缩方法通过在连续性方程中加入压强对虚拟时间的偏导数,从而把压力场和速度场耦合起来,解决了不可压缩流的计算问题.基于动态混合网格技术,利用双时间步方法耦合虚拟压缩方法来解决非定常不可压缩流的计算问题.为了加快每一虚拟时间步内的收敛速度,子迭代采用了高效的块LU-SGS方法,并且耦合了基于混合网格的多重网格方法.利用该方法数值模拟了不同雷诺数下的静止圆柱、振荡圆柱的绕流,得到了与实验和他人计算一致的结果.     

A local grid refinement method for three‐dimensional turbulent recirculating flows

A local grid refinement method is presented and applied to a three‐dimensional turbulent recirculating flow. It is based on the staggered grid arrangement. The computational domain is covered by block‐structured subgrids of different refinement levels. The exchange of information between the subgrids is fully conservative and all grids are treated implicitly. This allows for a simultaneous solution of one variable in all grids. All variables are stored in one‐dimensional arrays. The solver selected for the solution of the discretised finite difference equations is the preconditioned bi‐conjugate gradient (Bi‐CG) method. For the case examined (turbulent flow around a surface‐mounted cube), it was found that the latter method converges faster than the line solver. The locally refined mesh improved the accuracy of the pressure distribution on cube faces compared with a coarse mesh and yielded the same results as a fine single mesh, with a 62% gain in computer time. Copyright © 1999 John Wiley & Sons, Ltd.     

A stable and accurate projection method on a locally refined staggered mesh

In this paper, a robust projection method on a locally refined mesh is proposed for two‐ and three‐dimensional viscous incompressible flows. The proposed method is robust not only when the interface between two meshes is located in a smooth flow region but also when the interface is located in a flow region with large gradients and/or strong unsteadiness. In numerical simulations, a locally refined mesh saves many grid points in regions of relatively small gradients compared with a uniform mesh. For efficiency and ease of implementation, we consider a two‐level blocked structure, for which both of the coarse and fine meshes are uniform Cartesian ones individually. Unfortunately, the introduction of the two‐level blocked mesh results in an important but difficult issue: coupling of the coarse and fine meshes. In this paper, by properly addressing the issue of the coupling, we propose a stable and accurate projection method on a locally refined staggered mesh for both two‐ and three‐dimensional viscous incompressible flows. The proposed projection method is based on two principles: the linear interpolation technique and the consistent discretization of both sides of the pressure Poisson equation. The proposed algorithm is straightforward owing to the linear interpolation technique, is stable and accurate, is easy to extend from two‐ to three‐dimensional flows, and is valid even when flows with large gradients cross the interface between the two meshes. The resulting pressure Poisson equation is non‐symmetric on a locally refined mesh. The numerical results for a series of exact solutions for 2D and 3D viscous incompressible flows verify the stability and accuracy of the proposed projection method. The method is also applied to some challenging problems, including turbulent flows around particles, flows induced by impulsively started/stopped particles, and flows induced by particles near solid walls, to test the stability and accuracy. Copyright © 2010 John Wiley & Sons, Ltd.     

Unstructured pressure‐correction solver based on a consistent discretization of the Poisson equation

A finite volume incompressible flow solver is presented for three‐dimensional unsteady flows based on an unstructured tetrahedral mesh, with collocation of the flow variables at the cell vertices. The solver is based on the pressure‐correction method, with an explicit prediction step of the momentum equations followed by a Poisson equation for the correction step to enforce continuity. A consistent discretization of the Poisson equation was found to be essential in obtaining a solution. The correction step was solved with the biconjugate gradient stabilized (Bi‐CGSTAB) algorithm coupled with incomplete lower–upper (ILU) preconditioning. Artificial dissipation is used to prevent the formation of instabilities. Flow solutions are presented for a stalling airfoil, vortex shedding past a bridge deck and flow in model alveoli. Copyright © 2000 John Wiley & Sons, Ltd.     

Residual‐based turbulence models for moving boundary flows: hierarchical application of variational multiscale method and three‐level scale separation

This paper presents residual‐based turbulence models for problems with moving boundaries and interfaces. The method is developed via a hierarchical application of variational multiscale ideas and the models are cast in an arbitrary Lagrangian–Eulerian (ALE) frame to accommodate the deformation of domain boundaries. An overlapping additive decomposition of velocity and pressure fields into coarse and fine scale components leads to coarse and fine scale mixed‐field problems. The problem governing fine scales is subjected to a further decomposition of the fine scale velocity into overlapping components termed as fine scales level I and level II. In turn, in the bottom‐up integration of scales, the model for level II fine scales serves to stabilize the problem governing level I fine scales, and model for level I fields yields the turbulence models. From the computational perspective, the coarse scales are represented in terms of the standard Lagrange shape functions, whereas level I and level II scales are represented via quadratic and fourth order polynomial bubbles, respectively. Because of the bubble functions approach employed in the consistently derived fine scale models, the resulting method is free of any embedded or tunable parameters. The proposed turbulence models share a common feature with the LES models in that the largest scales in the flow are numerically resolved, whereas the subgrid scales are modeled. The method is applied to flow around a plunging airfoil at Re = 40,000, and results are compared with experimental and numerical data published in the literature. Also presented are the results for the plunging airfoil at Re = 60,000 to show the robustness and range of applicability of the method. Copyright © 2013 John Wiley & Sons, Ltd.     

A hybrid building‐block and gridless method for compressible flows

A hybrid building‐block Cartesian grid and gridless method is presented to compute unsteady compressible flows for complex geometries. In this method, a Cartesian mesh based on a building‐block grid is used as a baseline mesh to cover the computational domain, while the boundary surfaces are represented using a set of gridless points. This hybrid method combines the efficiency of a Cartesian grid method and the flexibility of a gridless method for the complex geometries. The developed method is used to compute a number of test cases to validate the accuracy and efficiency of the method. The numerical results obtained indicate that the use of this hybrid method leads to a significant improvement in performance over its unstructured grid counterpart for the time‐accurate solution of the compressible Euler equations. An overall speed‐up factor from six to more than one order of magnitude and a saving in storage requirements up to one order of magnitude for all test cases in comparison with the unstructured grid method are demonstrated. Copyright © 2008 John Wiley & Sons, Ltd.     

A posteriori bounds for linear functional outputs of Crouzeix–Raviart finite element discretizations of the incompressible Stokes problem

A finite element technique is presented for the efficient generation of lower and upper bounds to outputs which are linear functionals of the solutions to the incompressible Stokes equations in two space dimensions. The finite element discretization is effected by Crouzeix–Raviart elements, the discontinuous pressure approximation of which is central to this approach. The bounds are based upon the construction of an augmented Lagrangian: the objective is a quadratic ‘energy’ reformulation of the desired output, the constraints are the finite element equilibrium equations (including the incompressibility constraint), and the inter‐sub‐domain continuity conditions on velocity. Appealing to the dual max–min problem for appropriately chosen candidate Lagrange multipliers then yields inexpensive bounds for the output associated with a fine‐mesh discretization. The Lagrange multipliers are generated by exploiting an associated coarse‐mesh approximation. In addition to the requisite coarse‐mesh calculations, the bound technique requires the solution of only local sub‐domain Stokes problems on the fine mesh. The method is illustrated for the Stokes equations, in which the outputs of interest are the flow rate past and the lift force on a body immersed in a channel. Copyright © 2000 John Wiley & Sons, Ltd.     

An efficient semi‐implicit subgrid method for free‐surface flows on hierarchical grids

We present a new modelling strategy for improving the efficiency of computationally intensive flow problems in environmental free‐surface flows. The approach combines a recently developed semi‐implicit subgrid method with a hierarchical grid solution strategy. The method allows the incorporation of high‐resolution data on subgrid scale to obtain a more accurate and efficient hydrodynamic model. The subgrid method improves the efficiency of the hierarchical grid method by providing better solutions on coarse grids. The method is applicable to both steady and unsteady flows, but we particularly focus on river flows with steady boundary conditions. There, the combined hierarchical grid–subgrid method reduces the computational effort to obtain a steady state with factors up to 43. For unsteady models, the method can be used for efficiently generating accurate initial conditions on high‐resolution grids. Additionally, the method provides automatic insight in grid convergence. We demonstrate the efficiency and applicability of the method using a schematic test for the vortex shedding around a circular cylinder and a real‐world river case study. Copyright © 2015 John Wiley & Sons, Ltd.     

On high‐resolution schemes for solving unsteady compressible two‐phase dilute viscous flows

A high‐resolution numerical scheme based on the MUSCL–Hancock approach is developed to solve unsteady compressible two‐phase dilute viscous flow. Numerical considerations for the development of the scheme are provided. Several solvers for the Godunov fluxes are tested and the results lead to the choice of an exact Riemann solver adapted for both gaseous and dispersed phases. The accuracy of the scheme is proven step by step through specific test cases. These simulations are for one‐phase viscous flows over a flat plate in subsonic and supersonic regimes, unsteady flows in a low‐pressure shock tube, two‐phase dilute viscous flows over a flat plate and, finally, two‐phase unsteady viscous flows in a shock tube. The results are compared with well‐established analytical and numerical solutions and very good agreement is achieved. Copyright © 1999 John Wiley & Sons, Ltd.     

Relaxation system based sub-grid scale modelling for large eddy simulation of Burgers' equation

An implicit sub-grid scale model for large eddy simulation is presented by utilising the concept of a relaxation system for one dimensional Burgers' equation in a novel way. The Burgers' equation is solved for three different unsteady flow situations by varying the ratio of relaxation parameter (ε) to time step. The coarse mesh results obtained with a relaxation scheme are compared with the filtered DNS solution of the same problem on a fine mesh using a fourth-order CWENO discretisation in space and third-order TVD Runge-Kutta discretisation in time. The numerical solutions obtained through the relaxation system have the same order of accuracy in space and time and they closely match with the filtered DNS solutions.     

Discontinuity-capturing finite element computation of unsteady flow with adaptive unstructured mesh

A discontinuity-capturing scheme of finite element method (FEM) is proposed. The unstructured-grid technique combined with a new type of adaptive mesh approach is developed for both compressible and incompressible unsteady flows, which exhibits the capability of capturing the shock waves and/or thin shear layers accurately in an unsteady viscous flow at high Reynolds number. In particular, a new testing variable, i.e., the disturbed kinetic energyE, is suggested and used in the adaptive mesh computation, which is universally applicable to the capturing of both shock waves and shear layers in the inviscid flow and viscous flow at high Reynolds number. Based on several calculated examples, this approach has been proved to be effective and efficient for the calculations of compressible and incompressible flows. The project supported by the National Natural Science Foundation of China (10125210), the Hundred-Talent Programme of the Chinese Academy of Sciences and the Innovation Project of the Chinese Academy of Sciences (KJCX-SW-L04, KJCX2-SW-L2)     

Transient behavior of an axisymmetric turbulent jet

The prediction of the response of unsteady flows submitted to external excitation is a real challenge for the optimization of industrial processes. As the jet flow is a very basic turbulent flow related to mixing and entrainment phenomena via turbulent structure dynamic, we investigate the transient behavior of an axisymmetric jet submitted to a large and sudden decrease of its ejection velocity. The non stationary flow evolution is studied experimentally. Measurements along the jet axis based on pure ensemble averaging show clearly the convective motion of the perturbation and the adaptation of the local interaction to the local jet time scale. A transverse investigation in the non stationary region show that the mean flow and its turbulence is deeply affected during the local velocity decrease.     

三维非均匀不稳定渗流方程的二维不均匀粗化算法

在无源汇条件下,根据流过某一个横截面的流体流量等于流过这一横截面内所有精细网格的流体流量之和这一特点提出了粗化网格等效渗透率的计算方法。在粗化区内,利用直接解法求解二维渗流方程,再用这些解合成粗化网格的三维合成解,并由合成解计算粗化网格的等效渗透率。根据精度的要求采用了不均匀网格粗化,在流体流速大的区域采用精细网格。利用所得等效渗透率计算了粗化网格的某三维非均匀不稳定渗流场的压降解,结果表明三维非均匀不稳定渗流方程的二维不均匀粗化解非常逼近采用精细网格的解,但计算的速度比采用精细网格提高了80倍。     

Validation of adaptive unstructured hexahedral mesh computations of flow around a wind turbine airfoil

In this paper we investigate local adaptive refinement of unstructured hexahedral meshes for computations of the flow around the DU91 wind turbine airfoil. This is a 25% thick airfoil, found at the mid‐span section of a wind turbine blade. Wind turbine applications typically involve unsteady flows due to changes in the angle of attack and to unsteady flow separation at high angles of attack. In order to obtain reasonably accurate results for all these conditions one should use a mesh which is refined in many regions, which is not computationally efficient. Our solution is to apply an automated mesh adaptation technique. In this paper we test an adaptive refinement strategy developed for unstructured hexahedral meshes for steady flow conditions. The automated mesh adaptation is based on local flow sensors for pressure, velocity, density or a combination of these flow variables. This way the mesh is refined only in those regions necessary for high accuracy, retaining computational efficiency. A validation study is performed for two cases: attached flow at an angle of 6° and separated flow at 12°. The results obtained using our adaptive mesh strategy are compared with experimental data and with results obtained with an equally sized non‐adapted mesh. From these computations it can be concluded that for a given computing time, adapted meshes result in solutions closer to the experimental data compared to non‐adapted meshes for attached flow. Finally, we show results for unsteady computations. Copyright © 2005 John Wiley & Sons, Ltd.     

Convergence of steady and unsteady formulations for inviscid hovering rotor solutions

An upwind Euler solver is presented, and applied to multibladed lifting hovering rotor flow. These flows can be simulated as a steady case, in a blade‐fixed rotating co‐ordinate system. However, forward flight simulation will always require an unsteady solution. Hence, as a stepping stone in the development of a forward flight simulation tool, both explicit steady and implicit unsteady simulations of the same hovering case are presented. Convergence of the two approaches is examined and compared, in terms of residual history, cost, and solution evolution, as a means of both validating the unsteady formulation and considering implications for forward flight simulation. Consideration of the solution evolution and wake capturing shows that for hovering rotor cases, the unsteady and steady solutions are the same, but the unsteady solution is more expensive in terms of CPU time. It is also shown that for hover, the fewer real time‐steps taken per revolution the more efficient the implicit scheme is. However, this is a characteristic of the case, which results in smooth solution variation between time steps. It is also demonstrated that for rotary flow simulation, the global residual is not a useful quantity to assess convergence. The residual reaches a very low (constant in the implicit case) value while the solution is still evolving. Copyright © 2003 John Wiley & Sons, Ltd.