COMPUTATION OF DILUTE PARTICULATE LAMINAR FLOW OVER A BACKWARD-FACING STEP

Particle-laden flows are calculated for a classical laminar backward- facing step problem. The particle tracks are calculated using a recently developed exponential Lagrangian tracking scheme. The behaviour of the particle-laden flow is considered for various inlet for Reynolds number, Stokes numbers and void fractions. Doping the flow with low-Stokes-number particles has the effect of increasing the inlet inertia of the flow and this increases the strength of the recirculation behind the step. High-Stokes- number particles are dominated by gravitational effects which affect the flow accordingly. Differences between the single-phase flow and the particle-laden flows are therefore dependent on the Stokes number and increase linearly with void fraction.     

A NUMERICAL STUDY OF THE FLOW OVER ELLIPSOIDAL OBJECTS INSIDE A CYLINDRICAL TUBE

A numerical scheme is developed to obtain the flow field around one, two and five ellipsoidal objects inside a cylindrical tube. The scheme uses the Galerkin finite element technique and the primitive variable(u–v–p) formulation. The two-dimensional incompressible Navier–Stokes equations are solved numerically by using the direct mixed interpolation method. A Picard iteration scheme is used for the solution of the resulting system of non-linear algebraic equations. The computer code is verified by checking with known analytical solutions for the flow past a sphere. Results for the shear stress distributions along the ellipsoids, forces and drag coefficients are obtained for different geometric ratios and Reynolds numbers. Some of the intermediate computational results on the velocity fields developed are also reported.     

HIGH-ORDER-ACCURATE DISCRETIZATION STENCIL FOR AN ELLIPTIC EQUATION

The coefficients for a nine-point high-order-accurate discretization scheme for an elliptic equation ∇²u− γ²u=r₀ (∇² is the two-dimensional Laplacian operator) are derived. Examples with Dirichlet and Neumann boundary condtions are considered. In order to demonstrate the high-order accuracy of the method, numerical results are compared with exact solutions.     

METHOD–OF–LINES SOLUTION OF TIME–DEPENDENT TWO–DIMENSIONAL NAVIER–STOKES EQUATIONS

A novel approach to the development of a code for the solution of the time-dependent two-dimensional Navier–Stokes equations is described. The code involves coupling between the method of lines (MOL) for the solution of partial differential equations and a parabolic algorithm which removes the necessity of iterative solution on pressure and solution of a Poisson-type equation for the pressure. The code is applied to a test problem involving the solution of transient laminar flow in a short pipe for an incompressible Newtonian fluid. Comparisons show that the MOL solutions are in good agreement with the previously reported values. The proposed method described in this paper demonstrates the ease with which the Navier–Stokes equations can be solved in an accurate manner using sophisticated numerical algorithms for the solution of ordinary differential equations (ODEs).     

AN UNFACTORED IMPLICIT MOVING MESH METHOD FOR THE TWO-DIMENSIONAL UNSTEADY N–S EQUATIONS

An unfactored implicit time-marching method for the solution of the unsteady two-dimensional Reynolds-averaged thin layer Navier–Stokes equations is presented. The linear system arising from each implicit step is solved by the conjugate gradient squared (CGS) method with preconditioning based on an ADI factorization. The time-marching procedure has been used with a fast transfinite interpolation method to regenerate the mesh at each time step in response to the motion of the aerofoil. The main test cases examined are from the AGARD aeroelastic configurations and involve aerofoils oscillating rigidly in pitch. These test cases have been used to investigate the effect of various parameters, such as CGS tolerance and laminar/turbulent transition location, on the accuracy and efficiency of the method. Comparisons with available experimental data have been made for these cases. In order to illustrate the application of the mesh generator and flow solver to more general flows where the aerofoil deforms, results for an NACA 0012 aerofoil with an oscillating trailing edge flap are also shown.     

超声波燃气表矩形流道的雷诺修正系数仿真研究

在超声波流量计测量技术中, 雷诺修正系数相关的研究对于提高计量精度有重要作用. 为研究矩形流道的雷诺修正系数与雷诺数的关系, 对矩形流道在常温常压流量较小情况下进行仿真, 结果发现: 矩形流道层流状态下的雷诺修正系数与雷诺数呈线性相关. 保持压强、体积流量不变, 在不同温度下进行仿真及拟合, 结果表明: 在不同温度下雷诺修正系数与雷诺数的线性关系依然满足. 在上述实验基础上, 对矩形流道湍流状态下的雷诺修正系数与雷诺数关系进行研究, 通过改变温度、压强和体积流量进行仿真及拟合发现, 矩形流道湍流状态下雷诺修正系数与雷诺数呈非线性相关.     

Time-averaged flow characteristics and mixing properties of double-concentric jets

We examined the flow behaviors and mixing characteristics of double-concentric jets using laser-assisted smoke flow visualization method to analyze typical flow patterns and binary boundary detection technique to investigate jet spread width. Time-averaged velocity vectors, streamline patterns, velocity distributions, turbulence properties, and vorticity contours were analyzed using Particle Image Velocimetry (PIV). Topological flow patterns were analyzed to interpret the vortical flow structures. Mixing properties were investigated using a tracer-gas concentration detection method. Four characteristic modes were observed: annular flow dominated mode, transition mode, central jet dominated mode-low shear, and central jet dominated mode-high shear. The jets’ mixing properties were enhanced by two major phenomena: the merging of annular flow and central jet at the centerline and the large turbulence fluctuations produced in the flow field. The merging of the jets induced stagnation points on the central axis in the annular flow dominated mode, which caused reverse flow on the central axis and drastic turbulence fluctuations of the near field region. When the central jet penetrated the recirculation region in the other three modes, the stagnation points on the central axis and the reverse flow vanished. Therefore, the mixing behaviors were prominently enhanced in the annular flow dominated mode.     

Downstream evolution of junction flow three-component velocity fluctuations through the lens of the diagnostic plot

Experimental measurements of the streamwise and transverse velocity components have been acquired in three spanwise/wall-normal planes in the wakes of both a streamlined ‘wing’ and a bluff ‘wing’ junction. The ‘extended’ diagnostic plot introduced by Alfredsson et al. (2011) (see figure 3 therein) is used as a benchmark to locally evaluate the departure of turbulent wing-body junction flow wakes from ‘equilibrium’ boundary layers. Both obstacles produce a secondary flow of Prandtl’s first kind, which disrupts the equilibrium implied by the universality of the extended diagnostic plot. The plane wake of the obstacle itself (away from the junction) also disrupts this equilibrium. It is found that with downstream development the boundary layer eventually recovers to the base zero-pressure-gradient ‘equilibrium’, and that this recovery process emanates from the near-wall region. The transverse velocity components are also examined in “extended diagnostic” form, revealing that the wall-normal fluctuations recover to the zero-pressure-gradient case near the wall more rapidly than the wall-parallel components.     

On the inflation and deflation dynamics of liquid-filled,hyperelastic balloons

We derive a reduced-order model describing the inflation and deflation dynamics of a liquid-filled hyperelastic balloon, focusing on inviscid laminar flow and the extensional motion of the balloon. We initially study the flow and pressure fields for dictated motion of the solid, which throughout deflation are obtained by solving the potential problem. However, during inflation, flow separation creates a jet within the balloon, requiring a different approach. The analyses of both flow regimes lead to a simple piecewise model, describing the fluidic pressure during inflation and deflation, which is verified by finite element computations. We then use a variational approach to derive the equation describing the interaction between the extensional mode of the balloon and the entrapped fluid, yielding a nonlinear hybrid oscillator equation. Analytical and graphical investigations of the suggested model are presented, shedding light on its static and dynamic behaviour under different operating conditions. Our simplified model and its underlying assumptions are verified utilizing a fully coupled finite element scheme, showing excellent agreement.     

面向仿真数据的电网运行方式可视分析

对真实电网提出故障发生的假设并通过计算检验,观察故障对真实电网的影响。电网仿真能有效发现电网运行中的问题,被广泛应用于电网规划与电网运行方式调整。然而,在实际仿真计算过程中,往往需要依据经验调整电网配置和计算参数,以得到稳定、收敛的电网运行方式。大规模的仿真数据、调整的未知性加之分析的片面性会影响传统人工调整的准确性,为此提出了基于可视化的仿真计算调整方法。分析人员通过与可视分析系统的交互,可以快速认识电网的整体状态,发现问题元件,从而实现高效的仿真调整与计算。最后通过案例分析验证了本可视分析系统的有效性。