Modeling of viscous flows in two-dimensional turbomachinery cascade via viscous-inviscid interaction method

[Abstract]: A two-dimensional time-accurate time-marching viscous flow solver employing the viscous-inviscid interaction method suitable for turbomachinery applications is described. The inviscid main flow solver uses the second-order accurate cellvertex finite-volume spatial discretisation and fourth-order accurate Runge-Kutta temporal integration. The viscous effect due to boundary layer development on the blade surfaces and wakes are modelled using an independent one-dimensional boundary layer subroutine capable of modelling laminar, transition and fully turbulent flows. The solver has been applied to subsonic, transonic and supersonic flow in a cascade of nozzle blades. The results are compared with the experimental data and they showed very good agreemen

Rheology of Lamellar Liquid Crystals in Two and Three Dimensions: A Simulation Study

We present large scale computer simulations of the nonlinear bulk rheology of lamellar phases (smectic liquid crystals) at moderate to large values of the shear rate (Peclet numbers 10-100), in both two and three dimensions. In two dimensions we find that modest shear rates align the system and stabilise an almost regular lamellar phase, but high shear rates induce the nucleation and proliferation of defects, which in steady state is balanced by the annihilation of defects of opposite sign. The critical shear rate at onset of this second regime is controlled by thermodynamic and kinetic parameters; we offer a scaling analysis that relates the critical shear rate to a critical "capillary number" involving those variables. Within the defect proliferation regime, the defects may be partially annealed by slowly decreasing the applied shear rate; this causes marked memory effects, and history-dependent rheology. Simulations in three dimensions show instead shear-induced ordering even at the highest shear rates studied here. This suggests that the critical shear rate shifts markedly upward on increasing dimensionality. This may in part reflect the reduced constraints on defect motion, allowing them to find and annihilate each other more easily. Residual edge defects in the 3D aligned state mostly point along the flow velocity, an orientation impossible in two dimensions.Comment: 18 pages, 12 figure

Direct numerical simulation of backward-facing step flow at Ret = 395 and expansion ratio 2

Backward-facing step (BFS) constitutes a canonical configuration to study wallbounded flows subject to massive expansions produced by abrupt changes in geometry. Recirculation flow regions are common in this type of flow, driving the separated flow to its downstream reattachment. Consequently, strong adverse pressure gradients arise through this process, feeding flow instabilities. Therefore, both phenomena are strongly correlated as the recirculation bubble shape defines how the flow is expanded, and how the pressure rises. In an incompressible flow, this shape depends on the Reynolds value and the expansion ratio. The influence of these two variables on the bubble length is widely studied, presenting an asymptotic behaviour when both parameters are beyond a certain threshold. This is the usual operating point of many practical applications, such as in aeronautical and environmental engineering. Several numerical and experimental studies have been carried out regarding this topic. The existing simulations considering cases beyond the above-mentioned threshold have only been achieved through turbulence modelling, whereas direct numerical simulations (DNS) have been performed only at low Reynolds numbers. Hence, despite the great importance of achieving this threshold, there is a lack of reliable numerical data to assess the accuracy of turbulence models. In this context, a DNS of an incompressible flow over a BFS is presented in this paper, considering a friction Reynolds number (Reτ) of 395 at the inflow and an expansion ratio 2. Finally, the elongation of the Kelvin–Helmholtz instabilities along the shear layer is also studied.Postprint (published version

Wind and boundary layers in Rayleigh-Benard convection. Part 2: boundary layer character and scaling

The effect of the wind of Rayleigh-Benard convection on the boundary layers is studied by direct numerical simulation of an L/H=4 aspect-ratio domain with periodic side boundary conditions for Ra={10^5, 10^6, 10^7} and Pr=1. It is shown that the kinetic boundary layers on the top- and bottom plate have some features of both laminar and turbulent boundary layers. A continuous spectrum, as well as significant forcing due to Reynolds stresses indicates undoubtedly a turbulent character, whereas the classical integral boundary layer parameters -- the shape factor and friction factor (the latter is shown to be dominated by the pressure gradient) -- scale with Reynolds number more akin to laminar boundary layers. This apparent dual behavior is caused by the large influence of plumes impinging onto and detaching from the boundary layer. The plume-generated Reynolds stresses have a negligible effect on the friction factor at the Rayleigh numbers we consider, which indicates that they are passive with respect to momentum transfer in the wall-parallel direction. However, the effect of Reynolds stresses cannot be neglected for the thickness of the kinetic boundary layer. Using a conceptual wind model, we find that the friction factor C_f should scale proportional to the thermal boundary layer thickness as C_f ~ lambda_Theta, while the kinetic boundary layer thickness lambda_u scales inversely proportional to the thermal boundary layer thickness and wind Reynolds number lambda_u ~ lambda_Theta^{-1} Re^{-1}. The predicted trends for C_f and \lambda_u are in agreement with DNS results

Wind and boundary layers in Rayleigh-Benard convection. Part 2: boundary layer character and scaling

The effect of the wind of Rayleigh-Benard convection on the boundary layers is studied by direct numerical simulation of an L/H=4 aspect-ratio domain with periodic side boundary conditions for Ra={10^5, 10^6, 10^7} and Pr=1. It is shown that the kinetic boundary layers on the top- and bottom plate have some features of both laminar and turbulent boundary layers. A continuous spectrum, as well as significant forcing due to Reynolds stresses indicates undoubtedly a turbulent character, whereas the classical integral boundary layer parameters -- the shape factor and friction factor (the latter is shown to be dominated by the pressure gradient) -- scale with Reynolds number more akin to laminar boundary layers. This apparent dual behavior is caused by the large influence of plumes impinging onto and detaching from the boundary layer. The plume-generated Reynolds stresses have a negligible effect on the friction factor at the Rayleigh numbers we consider, which indicates that they are passive with respect to momentum transfer in the wall-parallel direction. However, the effect of Reynolds stresses cannot be neglected for the thickness of the kinetic boundary layer. Using a conceptual wind model, we find that the friction factor C_f should scale proportional to the thermal boundary layer thickness as C_f ~ lambda_Theta, while the kinetic boundary layer thickness lambda_u scales inversely proportional to the thermal boundary layer thickness and wind Reynolds number lambda_u ~ lambda_Theta^{-1} Re^{-1}. The predicted trends for C_f and \lambda_u are in agreement with DNS results

A numerical method for computing unsteady 2-D boundary layer flows

A numerical method for computing unsteady two-dimensional boundary layers in incompressible laminar and turbulent flows is described and applied to a single airfoil changing its incidence angle in time. The solution procedure adopts a first order panel method with a simple wake model to solve for the inviscid part of the flow, and an implicit finite difference method for the viscous part of the flow. Both procedures integrate in time in a step-by-step fashion, in the course of which each step involves the solution of the elliptic Laplace equation and the solution of the parabolic boundary layer equations. The Reynolds shear stress term of the boundary layer equations is modeled by an algebraic eddy viscosity closure. The location of transition is predicted by an empirical data correlation originating from Michel. Since transition and turbulence modeling are key factors in the prediction of viscous flows, their accuracy will be of dominant influence to the overall results