Efficient Large-eddy Simulations of Low Mach Number Flows Using Preconditioning and Multigrid

In the present paper, an implicit time accurate approach combined with multigrid, preconditioning and residual smoothing is used for the large-eddy simulation (LES) of low Mach number flow. In general, due to the restriction imposed on the time step by the physics of the flow, the advantage of an implicit method over an explicit one for LES is not obvious. It is shown that for the test cases considered in this paper, the present approach allows an efficiency gain of a factor 4–7 compared to the use of a purely explicit approach. The efficiency varies according to the test case, grid clustering, physical time step and requested residual drop. Numerical difficulties are catalogued and mitigatory procedures are introduced. Several problems with available experimental and DNS data are employed to verify the efficiency of the method.     

A unified approach to geometrically nonlinear analysis of tapered bonded joints and doublers

A unified approach for approximating the adhesive stresses in a bond line of a tapered bonded joint or doubler is delineated within the framework of a geometrically nonlinear analysis. The approach follows the Goland–Reissner solution method for a single-lap joint and involves a two-step analysis procedure. The approach also allows for the analysis of a tapered bonded joint and doubler with non-identical adherends. In the first step of the procedure, the two adherends are assumed to be rigidly bonded, and the nonlinear moment distribution along the joint is determined. Since the bending moment solution in this step is simple, it will be derived in closed-form using elementary functions. In the second step analysis, only the overlapped area of the joint is considered with the nonlinear bending moments obtained from the first step at the end of the overlap prescribed as one of its boundary conditions. This latter problem is then solved by using the multi-segment method of integration [Kalnins, A., 1964. Analysis of shell of revolutions subjected to symmetrical and non-symmetrical loads. Journal of Applied Mechanics 31, 1355–1365]. In contrast to the original Goland–Reissner solution method [Goland, M., Reissner, E., 1944. The stresses in cemented joints. Journal of Applied Mechanics 11, A17–A27], the second step analysis can be conducted within both geometrically linear theory and an approximate geometrically nonlinear theory.     

Use of two‐step split‐operator approach for 2D shallow water flow computation

The objective of this paper is to present a methodology of using a two‐step split‐operator approach for solving the shallow water flow equations in terms of an orthogonal curvilinear co‐ordinate system. This approach is in fact one kind of the so‐called fractional step method that has been popularly used for computations of dynamic flow. By following that the momentum equations are decomposed into two portions, the computation procedure involves two steps. The first step (dispersion step) is to compute the provisional velocity in the momentum equation without the pressure gradient. The second step (propagation step) is to correct the provisional velocity by considering a divergence‐free velocity field, including the effect of the pressure gradient. This newly proposed method, other than the conventional split‐operator methods, such as the projection method, considers the effects of pressure gradient and bed friction in the second step. The advantage of this treatment is that it increases flexibility, efficiency and applicability of numerical simulation for various hydraulic problems. Four cases, including back‐water flow, reverse flow, circular basin flow and unsteady flow, have been demonstrated to show the accuracy and practical application of the method. Copyright © 1999 John Wiley & Sons, Ltd.     

轴对称冲击有限元一致质量矩阵迭代解

给出高速冲击动力有限元一致质量矩阵解的迭代过程，即把集总质量矩阵解作为初值进行有限次迭代得到满足工程需要的一致质量矩阵近似解．实际算例说明，一致质量近似解较集总质量解改善了对应力波传播过程的分析，而且在高速冲击计算中能给出与实验接近的计算结果．     

An alternative approach to integrating plasticity relations

A new plasticity integration algorithm is proposed based upon observations from the closed form integration of a generalized quadratic yield function over a single time step. The key to the approach is specification of the normal to the plastic flow potential as a function of the current state and strain increment. This uniquely defines the direction of the stress tensor for a convex, non-faceted flow potential. The stress magnitude and plastic strain increment are computed to satisfy the yield function. A non-quadratic, isotropic, associative flow model is coded to demonstrate accuracy and time step convergence following a step change in loading path. The model is used in additional simulations of strain localization in an expanding ring and a perforated plate.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

A three dimensional implicit immersed boundary method with application

Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.     

An Eulerian approach for simulating frictional heating in disc-pad systems

Thermal stresses as a result from frictional heating must be considered when designing disc brakes, clutches or other rotating machine components with sliding contact conditions. The rotational symmetry of the disc in these kind of applications makes it possible to model these systems using an Eulerian approach instead of a Lagrangian framework. In this paper such an approach is developed and implemented. The disc is formulated in an Eulerian frame where the convective terms are defined by the angular velocity. By utilizing the Eulerian framework, a node-to-node formulation of the contact interface is obtained, producing most accurate frictional heat power solutions. The energy balance of the interface is postulated by introducing an interfacial temperature. Both frictional power and contact conductances are included in this energy balance. The contact problem is solved by a non-smooth Newton method. By adopting the augmented Lagrangian approach, this is done by rewriting Signorini’s contact conditions to an equivalent semi-smooth equation. The heat transfer in the disc is discretized by a Petrov–Galerkin approach, i.e. the numerical difficulties due to the non-symmetric convective matrix appearing in a pure Galerkin discretization is treated by following the streamline-upwind approach. In such manner a stabilization is obtained by adding artificial conduction along the streamlines. For each time step the thermo-elastic contact problem is first solved for the temperature field from the previous time step. Then, the heat transfer problem is solved for the corresponding frictional power. In such manner a temperature history is obtained sequentially via the trapezoidal rule. In particular the parameter is set such that both the Crank–Nicolson and the Galerkin methods are utilized. The method seems very promising. This is demonstrated by solving a two-dimensional benchmark as well as a real disc brake system in three dimensions.     

A diffuse-interface immersed boundary method for simulation of compressible viscous flows with stationary and moving boundaries

In this paper, a diffuse-interface immersed boundary method (IBM) is proposed for simulation of compressible viscous flows with stationary and moving boundaries. In the method, the solution of flow field and the implementation of boundary conditions are decoupled into two steps by applying the fractional step technique, ie, the predictor step and the corrector step. Firstly, in the predictor step, the intermediate flow field is resolved by a recently developed gas kinetic flux solver (GKFS) without consideration of the solid boundary. The GKFS is a finite volume approach that solves the Navier-Stokes equations for the flow variables at cell centers. In GKFS, the inviscid and viscous fluxes are evaluated as a single entity by reconstructing the local solution of continuous Boltzmann equation. Secondly, in the corrector step, the intermediate flow field is corrected by the present diffuse-interface IBM. During this process, the velocity field is firstly corrected by the implicit boundary condition–enforced IBM so that the no-slip boundary condition can be accurately satisfied. After that, the density correction is made by an iterative approach with the help of the continuity equation. Finally, the correction of the temperature field is made in the same way as that of the velocity field. Good agreements between the present simulations and the reference data in literature demonstrate the reliability of the proposed method.     

非结构重叠网格方法研究及应用

航空、航天和兵器技术等领域的研究中存在大量包含运动边界的流场。非结构重叠网格方法是一种高效的处理动边界问题的新方法。围绕相对运动的每个物体单独生成非结构网格,在网格重叠区域通过搜索和插值完成网格系之间的信息传递,提出了动态八叉树搜索算法,发展了绝对坐标系和相对坐标系相结合的流场求解方式,采用二阶精度Van Leer/Hanel格式和四阶Runge-Kutta法分别进行空间和时间离散,形成了一种新的非结构重叠网格算法。对三维Riemann问题的求解结果与精确解能很好吻合,证明了本文的重叠网格算法具有较好的时空离散精度和插值精度。对7.62mm步枪射击过程进行了数值模拟,描述了弹丸离开膛口后膛口流场的发展过程,与实验结果体现的发展过程较为吻合,验证了本文提出的非结构网格算法体系具有较好的计算性能,是研究含动边界复杂流场的一种有效手段。     

Sequential linearization method for viscous/elastic heterogeneous materials

The paper addresses the problem of suitable approximation of the interaction between phases in heterogeneous materials that exhibit both viscous and elastic properties. A novel approach is proposed in which linearized subproblems for an inhomogeneity–matrix system with viscous or elastic interaction rules are solved sequentially within one incremental step. It is demonstrated that in the case of a self-consistent averaging scheme, an additional accommodation subproblem, besides purely viscous and elastic subproblems, is to be solved in order to estimate the material response satisfactorily. By examples of an isotropic two-phase material it is shown that the proposed approach provides acceptable predictions in comparison with the existing models.     

Implicit time accurate simulation of unsteady flow

Implicit time integration was studied in the context of unsteady shock‐boundary layer interaction flow. With an explicit second‐order Runge–Kutta scheme, a reference solution to compare with the implicit second‐order Crank–Nicolson scheme was determined. The time step in the explicit scheme is restricted by both temporal accuracy as well as stability requirements, whereas in the A‐stable implicit scheme, the time step has to obey temporal resolution requirements and numerical convergence conditions. The non‐linear discrete equations for each time step are solved iteratively by adding a pseudo‐time derivative. The quasi‐Newton approach is adopted and the linear systems that arise are approximately solved with a symmetric block Gauss–Seidel solver. As a guiding principle for properly setting numerical time integration parameters that yield an efficient time accurate capturing of the solution, the global error caused by the temporal integration is compared with the error resulting from the spatial discretization. Focus is on the sensitivity of properties of the solution in relation to the time step. Numerical simulations show that the time step needed for acceptable accuracy can be considerably larger than the explicit stability time step; typical ratios range from 20 to 80. At large time steps, convergence problems that are closely related to a highly complex structure of the basins of attraction of the iterative method may occur. Copyright © 2001 John Wiley & Sons, Ltd.     

Direct numerical simulation of flow over a forward-facing step – Flow structure and aeroacoustic source regions

A direct numerical simulation of the flow over a forward-facing step at a Reynolds number of 8000 based on the step height is presented. Calculations were performed using second-order finite volume discretisation in space on co-located meshes. A hybrid calculation approach based on Lighthill’s acoustic analogy is explained. Results of the simulation are intended to be used as a database for the validation of different discretisation schemes for the flow computation and simulation approaches for the calculation of sound radiation using a hybrid approach. Turbulent statistics are presented along with aeroacoustic source regions. Strong and weak forms of the aeroacoustic source term are presented and compared. For visualization purposes, the strong form is more suitable, whereas for the calculation of sound radiation both forms can be used. From the visualization of the aeroacoustic sources, it can be seen that they mainly concentrate on the region of the leading edge of the step and the shear layer close to the step.     

Detailed parametric design methodology for hydrodynamics of liquid–solid circulating fluidized bed using design of experiments

A design-of-experiments methodology is used to develop a statistical model for the prediction of the hydrodynamics of a liquid–solid circulating fluidized bed. To illustrate the multilevel factorial design approach, a step by step methodology is taken to study the effects of the interactions among the independent factors considered on the performance variables. A multilevel full factorial design with three levels of the two factors and five levels of the third factor has been studied. Various statistical models such as the linear, two-factor interaction, quadratic, and cubic models are tested. The model has been developed to predict responses, viz., average solids holdup and solids circulation rate. The validity of the developed regression model is verified using the analysis of variance. Furthermore, the model developed was compared with an experimental dataset to assess its adequacy and reliability. This detailed statistical design methodology for non-linear systems considered here provides a very important tool for design and optimization in a cost-effective approach.     

概率功能度量求解的共轭梯度步长调节法

功能度量法(PMA)由于其稳定高效的特点,适用于概率结构优化设计中概率约束的评定。PMA中改进均值法常用于求解概率功能度量,针对其求解高度非线性功能函数时出现周期振荡和混沌等不收敛现象,提出了一种新的共轭梯度步长调节法(CGS)。该方法基于RMIL共轭搜索方向和自适应步长调节策略提出,新的共轭搜索方向在保证收敛性的前提下加速了迭代进程,而自适应步长调节策略无需了解功能函数凹凸性及非线性程度等先验信息,无需确定步长的合适取值。通过限定步长准则自动选取初始步长,并随迭代过程不断调节,直至最终收敛。多个算例表明,与其他求解方法相比,本文的共轭梯度步长调节法更加高效且稳健。     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

分层流体中内孤立波在台阶上的反射和透射

基于匹配渐近展开和格林函数的方法，研究了两层流体系统中内孤立波在台阶地形上透射、 反射及其分裂的演化特征. 通过保角变换和求解奇异Fredholm积分方程，获得了反映地形 效应对Boussinesq方程影响的约化边界条件，藉此建立了KdV演化方程的``初值'问题， 根据散射反演理论获得了反射波和透射波的解析表达式. 分析结果表明：上下流体层的厚度 比、密度比以及台阶高度对于反射和透射波振幅及其分裂具有显著的影响. 尤其当上层流体 厚度小于下层厚度时，由于存在临界点，在其附近反射波的幅值随台阶高度的演化由单调增 变为单调减，透射波的幅值由单调减变为单调增；上台阶的反射波与入射波反相，其最大幅 值可达到入射波的数倍；此外，下台阶反射波也可发展为单支孤立波，它区别于单层流体中 反射波仅为衰减的振荡波列.     

3-D numerical simulation of fracture processes in heterogeneous brittle materials

By using the lattice model combined with finite element methods and statistical techniques, a numerical approach is developed to establish mechanical models of three-dimensional heterogeneous brittle materials. A special numerical code is introduced, in which a lattice model and statistical approaches are used to simulate the initial heterogeneity of material properties. The size of displacement-load step is adaptively determined so that only few elements would fail in each load step. When the tensile principal strain in an element exceeds the ultimate strain of this element, the element is considered broken and its Young's modulus is set to be very low. Some important behaviors of heterogeneous brittle materials are indicated using this code. Load-displacement curves and figures of three-dimensional fracture patterns are also numerically obtained, which are similar to those observed in laboratory tests.