基于离散伴随方法的透平叶栅气动优化设计

本文研究并给出了基于离散伴随理论和自动微分技术构建离散伴随系统的方法、伴随系统的求解策略以及基于离散伴随方法的透平叶栅气动优化设计流程,建立了相应的优化设计系统。利用该优化系统在无黏环境下,以叶栅通道进出口的熵增率为目标函数、以叶栅通道内的质量流量为约束,对某二维跨音速透平叶栅进行了气动优化设计。与优化前相比,优化后透平叶栅进出口熵增率减少8.82%,质量流量变化幅度小于0.003%。优化结果表明,本文提出的优化系统能够有效改善透平叶栅的气动优化性能,验证了本文提出的基于离散伴随方法的透平叶栅气动优化设计方法的正确性与有效性。     

Adjoint wall functions: A new concept for use in aerodynamic shape optimization

The continuous adjoint method for the computation of sensitivity derivatives in aerodynamic optimization problems of steady incompressible flows, modeled through the k–ε turbulence model with wall functions, is presented. The proposed formulation leads to the adjoint equations along with their boundary conditions by introducing the adjoint to the friction velocity. Based on the latter, an adjoint law of the wall that bridges the gap between the solid wall and the first grid node off the wall is proposed and used during the solution of the system of adjoint (to both the mean flow and turbulence) equations. Any high Reynolds turbulence model, other than the k–ε one used in this paper, could also profit from the proposed adjoint wall function technique. In the examined duct flow problems, where the total pressure loss due to viscous effects is used as objective function, emphasis is laid on the accuracy of the computed sensitivity derivatives, rather than the optimization itself. The latter might rely on any descent method, once the objective function gradient has accurately been computed.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

Aerodynamic design optimization by using a continuous adjoint method

This paper presents the fundamentals of a continuous adjoint method and the applications of this method to the aerodynamic design optimization of both external and internal flows.General formulation of the continuous adjoint equations and the corresponding boundary conditions are derived.With the adjoint method,the complete gradient information needed in the design optimization can be obtained by solving the governing flow equations and the corresponding adjoint equations only once for each cost function,regardless of the number of design parameters.An inverse design of airfoil is firstly performed to study the accuracy of the adjoint gradient and the effectiveness of the adjoint method as an inverse design method.Then the method is used to perform a series of single and multiple point design optimization problems involving the drag reduction of airfoil,wing,and wing-body configuration,and the aerodynamic performance improvement of turbine and compressor blade rows.The results demonstrate that the continuous adjoint method can efficiently and significantly improve the aerodynamic performance of the design in a shape optimization problem.     

Lattice Boltzmann simulations on GPUs with ESPResSo

For the dynamics of macromolecules in solution, hydrodynamic interactions mediated by the solvent molecules often play an important role, although one is not interested in the dynamics of the solvent itself. In computer simulations one can therefore save a large amount of computer time by replacing the solvent with a lattice fluid. The macromolecules are propagated by Molecular Dynamics (MD), while the fluid is governed by the fluctuating Lattice-Boltzmann (LB) equation. We present a fluctuating LB implementation for a single graphics card (GPU) coupled to a MD simulation running on conventional processors (CPUs). Particular emphasis lies on the optimization of the combined code. In our implementation, the LB update is performed in parallel with the force calculation on the CPU, which often completely hides the additional computational cost of the LB. Compared to our parallel LB implementation on a conventional quad-core CPU, the GPU LB is 50 times faster, and we show that a whole commodity cluster with Infiniband interconnnect cannot outperform a single GPU in strong scaling. The presented code is part of the open source simulation package ESPResSo ().     

A Simplified Linearized Lattice Boltzmann Method for Acoustic Propagation Simulation

A simplified linearized lattice Boltzmann method (SLLBM) suitable for the simulation of acoustic waves propagation in fluids was proposed herein. Through Chapman–Enskog expansion analysis, the linearized lattice Boltzmann equation (LLBE) was first recovered to linearized macroscopic equations. Then, using the fractional-step calculation technique, the solution of these linearized equations was divided into two steps: a predictor step and corrector step. Next, the evolution of the perturbation distribution function was transformed into the evolution of the perturbation equilibrium distribution function using second-order interpolation approximation of the latter at other positions and times to represent the nonequilibrium part of the former; additionally, the calculation formulas of SLLBM were deduced. SLLBM inherits the advantages of the linearized lattice Boltzmann method (LLBM), calculating acoustic disturbance and the mean flow separately so that macroscopic variables of the mean flow do not affect the calculation of acoustic disturbance. At the same time, it has other advantages: the calculation process is simpler, and the cost of computing memory is reduced. In addition, to simulate the acoustic scattering problem caused by the acoustic waves encountering objects, the immersed boundary method (IBM) and SLLBM were further combined so that the method can simulate the influence of complex geometries. Several cases were used to validate the feasibility of SLLBM for simulation of acoustic wave propagation under the mean flow.     

基于离散伴随方法的透平叶栅反设计

研究构建了基于离散伴随方法的叶轮机械叶栅气动反设计系统,将离散伴随系统从无黏环境扩展到了黏性环境,编程实现了黏性离散伴随求解器;改善了叶栅参数化方式并重新编程实现了叶栅参数化程序,解决了叶栅参数化过程中叶栅尾缘附近区域型线波动的问题。利用该系统对某二维跨声速透平叶栅在给定叶型壁面目标压力分布的情况下,通过构造目标函数将叶栅反设计问题转化为气动优化设计问题,成功进行了气动压力反设计。结果证明本文建立的叶栅反设计系统能够有效进行压力反设计,验证了本文建立的基于离散伴随方法叶轮机械叶栅气动反设计方法的正确性与有效性。     

耦合不可压流场输运方程的格子Boltzmann方法研究

基于格子Boltzmann方法,提出了求解耦合不可压缩流场输运方程的一种改进数值方法. 该方法使用格子Boltzmann方法求解流场方程,并根据流场格子模型的密度分布函数构建了输运方程的二阶离散格式. 通过二维平板通道流场输运系统验证了该方法的有效性.数值结果表明,该方法可以有效地减少计算过程中出现的非物理耗散, 并克服了传统模型所需巨大存储量的缺点.     

An adjoint method for shape optimization in unsteady viscous flows

A new method for shape optimization for unsteady viscous flows is presented. It is based on the continuous adjoint approach using a time accurate method and is capable of handling both inverse and direct objective functions. The objective function is minimized or maximized subject to the satisfaction of flow equations. The shape of the body is parametrized via a Non-Uniform Rational B-Splines (NURBS) curve and is updated by using the gradients obtained from solving the flow and adjoint equations. A finite element method based on streamline-upwind Petrov/Galerkin (SUPG) and pressure stabilized Petrov/Galerkin (PSPG) stabilization techniques is used to solve both the flow and adjoint equations. The method has been implemented and tested for the design of airfoils, based on enhancing its time-averaged aerodynamic coefficients. Interesting shapes are obtained, especially when the objective is to produce high performance airfoils. The effect of the extent of the window of time integration of flow and adjoint equations on the design process is studied. It is found that when the window of time integration is insufficient, the gradients are most likely to be erroneous.     

Lattice Boltzmann Model for the Incompressible Navier–Stokes Equation

In this paper a lattice Boltzmann (LB) model to simulate incompressible flow is developed. The main idea is to explicitly eliminate the terms of o(M ²), where M is the Mach number, due to the density fluctuation in the existing LB models. In the proposed incompressible LB model, the pressure p instead of the mass density ρ is the independent dynamic variable. The incompressible Navier–Stokes equations are derived from the incompressible LB model via Chapman–Enskog procedure. Numerical results of simulations of the plane Poiseuille flow driven either by pressure gradient or a fixed velocity profile at entrance as well as of the 2D Womersley flow are presented. The numerical results are found to be in excellent agreement with theory.     

A study on the inclusion of body forces in the lattice Boltzmann BGK equation to recover steady-state hydrodynamics

When the lattice Boltzmann (LB) method is used to solve hydrodynamic problems containing a body force term varying in space and/or time, its modelling at the mesoscopic scale must be verified in terms of consistency in order to avoid the appearance of non-hydrodynamic error terms at the macroscopic scale. In the present work it is shown that the modelling of spatially varying steady body force terms in the LB equation must be different from the time-dependent case, when a steady-state flow solution is sought. For that, the Chapman-Enskog analysis is used to derive the LB body force model for the LB BGK equations in a steady-state flow problem. The theoretical findings are supported by numerical tests performed on two different 2D steady-state laminar flows driven by spatially varying body forces with known analytical solutions.     

Lattice Boltzmann method on quadtree grids for simulating fluid flow through porous media: A new automatic algorithm

During the past two decades, the lattice Boltzmann (LB) method has been introduced as a class of computational fluid dynamic methods for fluid flow simulations. In this method, instead of solving the Navier Stocks equation, the Boltzmann equation is solved to simulate the flow of a fluid. This method was originally developed based on uniform grids. However, in order to model complex geometries such as porous media, it can be very slow in comparison with other techniques such as finite differences and finite elements. To eliminate this limitation, a number of studies have aimed to formulate the lattice Boltzmann on the unstructured grids. This paper deals with simulating fluid flow through a synthetic porous medium using the LB method and on the quadtree grid structure. To this end, the LB method was used on nonuniform grids coupled with a technique for image reconstruction which resulted in the quadtree grids for simulation of fluid flow through porous media. Accuracy and efficiency of this algorithm is compared against the conventional LB method based on uniform grids. While the decrease in computational time in the proposed LB method on nonuniform grids is found to be significant regarding the size of the initial and reconstructed images, the same level of accuracy is obtained when compared with the conventional LB method on uniform grids.     

Error analysis and correction for Lattice Boltzmann simulated flow conductance in capillaries of different shapes and alignments

We study the relative error in conductance calculations, for simulated flow of a single component single phase fluid through a capillary in three dimensions, by the Lattice Boltzmann (LB) method with bounce-back boundary conditions. The relative error with respect to analytical results for capillary cross-sections of circular, triangular and square shapes are calculated as a function of the cross-section diameter, a, and for different alignment of the cross-section relative to the underlying lattice grid. It is shown, when the shapes are not aligned perfectly to the lattice, that the relative error decreases systematically with the size, a, as ～1/a when a is evaluated by mapping the computed cross-sectional area, in terms of the enclosed number of grid points, to the respective geometrical shapes concerned. For perfectly aligned geometries, viz. the square capillary aligned to the LB lattice grid or rotated with its side along the diagonal of the LB grid, the relative error decreases as ～1/a². A simple method is suggested to locate the boundary wall depending on its orientation relative to the grid, such that the exact conductance of the new shape matches the LB computed conductance.     

膜生物反应器内生化反应过程的格子Boltzmann模拟

采用介观尺度的格子Boltzmann数值算法研究平板式膜生物反应器内生物膜结构对生化反应及流动传输的影响.利用四参数随机生成法重构不同多孔结构的生物膜,计算中耦合多块模型获得局部详细信息且提高计算效率.结果表明：生物膜结构稳定条件下,孔隙率增大有利于主流区与生物膜内的物质传输,可以提高底物降解效率;孔隙率一定条件下,改变生物膜的孔隙结构分布可以加强传质,强化生化反应过程,提高底物降解效率.     

Local-in-time adjoint-based method for design optimization of unsteady flows

We present a new local-in-time discrete adjoint-based methodology for solving design optimization problems arising in unsteady aerodynamic applications. The new methodology circumvents storage requirements associated with the straightforward implementation of a global adjoint-based optimization method that stores the entire flow solution history for all time levels. This storage cost may quickly become prohibitive for large-scale applications. The key idea of the local-in-time method is to divide the entire time interval into several subintervals and to approximate the solution of the unsteady adjoint equations and the sensitivity derivative as a combination of the corresponding local quantities computed on each time subinterval. Since each subinterval contains relatively few time levels, the storage cost of the local-in-time method is much lower than that of the global methods, thus making the time-dependent adjoint optimization feasible for practical applications. Another attractive feature of the new technique is that the converged solution obtained with the local-in-time method is a local extremum of the original optimization problem. The new method carries no computational overhead as compared with the global implementation of adjoint-based methods. The paper presents a detailed comparison of the global- and local-in-time adjoint-based methods for design optimization problems governed by the unsteady compressible 2-D Euler equations.     

Polar Coordinate Lattice Boltzmann Kinetic Modeling of Detonation Phenomena

A novel polar coordinate lattice Boltzmann kinetic model for detonation phenomena is presented and applied to investigate typical implosion and explosion processes. In this model, the change of discrete distribution function due to local chemical reaction is dynamically coupled into the modified lattice Boltzmann equation which could recover the Navier-Stokes equations, including contribution of chemical reaction, via the Chapman-Enskog expansion. For the numerical investigations, the main focuses are the nonequilibrium behaviors in these processes. The system at the disc center is always in its thermodynamic equilibrium in the highly symmetric case. The internal kinetic energies in different degrees of freedom around the detonation front do not coincide. The dependence of the reaction rate on the pressure, influences of the shock strength and reaction rate on the departure amplitude of the system from its local thermodynamic equilibrium are probed.     

On the Three-Dimensional Central Moment Lattice Boltzmann Method

A three-dimensional (3D) lattice Boltzmann method based on central moments is derived. Two main elements are the local attractors in the collision term and the source terms representing the effect of external and/or self-consistent internal forces. For suitable choices of the orthogonal moment basis for the three-dimensional, twenty seven velocity (D3Q27), and, its subset, fifteen velocity (D3Q15) lattice models, attractors are expressed in terms of factorization of lower order moments as suggested in an earlier work; the corresponding source terms are specified to correctly influence lower order hydrodynamic fields, while avoiding aliasing effects for higher order moments. These are achieved by successively matching the corresponding continuous and discrete central moments at various orders, with the final expressions written in terms of raw moments via a transformation based on the binomial theorem. Furthermore, to alleviate the discrete effects with the source terms, they are treated to be temporally semi-implicit and second-order, with the implicitness subsequently removed by means of a transformation. As a result, the approach is frame-invariant by construction and its emergent dynamics describing fully 3D fluid motion in the presence of force fields is Galilean invariant. Numerical experiments for a set of benchmark problems demonstrate its accuracy.     

Adjoint sensitivity computations for an embedded-boundary Cartesian mesh method

We present a new approach for the computation of shape sensitivities using the discrete adjoint and flow-sensitivity methods on Cartesian meshes with general polyhedral cells (cut-cells) at the wall boundaries. By directly linearizing geometric constructors of the cut-cells, an efficient and robust computation of shape sensitivities is achieved for problems governed by the Euler equations. The accuracy of the linearization is verified by the use of a model problem with an exact solution. Verification studies show that the convergence rate of gradients is second-order for design variables that do not alter the boundary shape, and is reduced to first-order for shape design problems. The approach is applied to several three-dimensional problems, including inverse design and shape optimization of a re-entry capsule in hypersonic flow. The results show that reliable approximations of the gradient are obtained in all cases. The approach is well-suited for geometry control via computer-aided design, and is especially effective for conceptual design studies with complex geometry where fast turn-around time is required.     

Lattice Boltzmann simulation of solid particles suspended in fluid

The lattice Boltzmann method, an alternative approach to solving a fluid flow system, is used to analyze the dynamics of particles suspended in fluid. The interaction rule between the fluid and the suspended particles is developed for real suspensions where the particle boundaries are treated as no-slip impermeable surfaces. This method correctly and accurately determines the dynamics of single particles and multi-particles suspended in the fluid. With this method, computational time scales linearly with the number of suspensions,N, a significant advantage over other computational techniques which solve the continuum mechanics equations, where the computational time scales asN ³. Also, this method solves the full momentum equations, including the inertia terms, and therefore is not limited to low particle Reynolds number.