Benham G P, Hewitt I J, Please C P, Bird P A D
1Mathematical Institute, University of Oxford, Andrew Wiles Building, Radcliffe Observatory Quarter, Woodstock Road, Oxford, OX2 6GG UK.
VerdErg Renewable Energy Limited, 6 Old London Rd, Kingston upon Thames, KT2 6QF UK.
J Eng Math. 2018;113(1):65-92. doi: 10.1007/s10665-018-9974-6. Epub 2018 Nov 13.
A simplified model is used to identify the diffuser shape that maximises pressure recovery for several classes of non-uniform inflow. We find that optimal diffuser shapes strike a balance between not widening too soon, as this accentuates the non-uniform flow, and not staying narrow for too long, which is detrimental for wall drag. Three classes of non-uniform inflow are considered, with the axial velocity varying across the width of the diffuser entrance. The first case has inner and outer streams of different speeds, with a velocity jump between them that evolves into a shear layer downstream. The second case is a limiting case when these streams are of similar speed. The third case is a pure shear profile with linear velocity variation between the centre and outer edge of the diffuser. We describe the evolution of the time-averaged flow profile using a reduced mathematical model that has been previously tested against experiments and computational fluid dynamics models. The model consists of integrated mass and momentum equations, where wall drag is treated with a friction factor parameterisation. The governing equations of this model form the dynamics of an optimal control problem where the control is the diffuser channel shape. A numerical optimisation approach is used to solve the optimal control problem and Pontryagin's maximum principle is used to find analytical solutions in the second and third cases. We show that some of the optimal diffuser shapes can be well approximated by piecewise linear sections. This suggests a low-dimensional parameterisation of the shapes, providing a structure in which more detailed and computationally expensive turbulence models can be used to find optimal shapes for more realistic flow behaviour.
使用一个简化模型来确定几种非均匀入流情况下能使压力恢复最大化的扩压器形状。我们发现,最优扩压器形状要在不过早变宽(因为这会加剧流动不均匀性)和不过长时间保持狭窄(这对壁面阻力不利)之间取得平衡。考虑了三类非均匀入流情况,轴向速度在扩压器入口宽度上变化。第一种情况是内流和外流速度不同,它们之间存在速度跃变,下游会演变成一个剪切层。第二种情况是这些流速度相近时的极限情况。第三种情况是在扩压器中心和外边缘之间具有线性速度变化的纯剪切剖面。我们使用一个简化数学模型来描述时间平均流剖面的演变,该模型先前已通过实验和计算流体动力学模型进行了验证。该模型由积分质量和动量方程组成,其中壁面阻力采用摩擦系数参数化处理。此模型的控制方程构成了一个最优控制问题的动力学,其中控制量是扩压器通道形状。采用数值优化方法来求解最优控制问题,并利用庞特里亚金极大值原理在第二种和第三种情况下找到解析解。我们表明,一些最优扩压器形状可以用分段线性段很好地近似。这表明形状可以进行低维参数化,提供一种结构,在其中可以使用更详细且计算成本更高的湍流模型来找到更符合实际流动行为的最优形状。