No Arabic abstract
This paper presents flow simulation results of the EUROLIFT DLR-F11 multi-element wing configuration, obtained with a highly scalable finite element solver, PHASTA. This work was accomplished as a part of the 2nd high lift prediction workshop. In-house meshes were constructed with increasing mesh density for analysis. A solution adaptive approach was used as an alternative and its effectiveness was studied by comparing its results with the ones obtained with other meshes. Comparisons between the numerical solution obtained with unsteady RANS turbulence model and available experimental results are provided for verification and discussion. Based on the observations, future direction for adaptive research and simulations with higher fidelity turbulence models is outlined.
Capillary phenomena are involved in many industrial processes, especially those dealing with composite manufacturing. However, their modelling is still challenging. Therefore, a finite element setting is proposed to better investigate this complex issue. The variational formulation of a liquid-air Stokes system is established, while the solid substrate is described through boundary conditions. Expressing the weak form of Laplaces law over liquid-air, liquid-solid and air-solid interfaces, leads to a natural enforcement of the mechanical equilibrium over the wetting line, without imposing explicitly the contact angle itself. The mechanical problem is discretized by using finite elements, linear both in velocity and pressure, stabilized with a variational multiscale method, including the possibility of enrichment of the pressure space. The moving interface is captured by a Level-Set methodology, combined with a mesh adaptation technique with respect to both pressure and level-set fields. Our methodology can simulate capillary-driven flows in 2D and 3D with the desired precision: droplet spreading, droplet coalescence, capillary rise. In each case, the equilibrium state expected in terms of velocity, pressure and contact angle is reached.
The aerodynamic performance of the high-lift configuration greatly influences the safety and economy of commercial aircraft. Accurately predicting the aerodynamic performance of the high-lift configuration, especially the stall behavior, is important for aircraft design. However, the complex flow phenomena of high-lift configurations pose substantial difficulties to current turbulence models. In this paper, a three-equation k-(v^2)-{omega} turbulence model for the Reynolds-averaged Navier-Stokes equations is used to compute the stall behavior of high-lift configurations. A separated shear layer fixed function is implemented in the turbulence model to better capture the nonequilibrium characteristics of turbulence. Different high-lift configurations, including the two-dimensional multielement NLR7301 and Omar airfoils and a complex full-configuration model (JAXA Standard Model), are numerically tested. The results indicate that the effect of the nonequilibrium characteristics of turbulence is significant in the free shear layer, which is key to accurately predicting the stall behavior of high-lift devices. The modified SPF k-(v^2 )-{omega} model is more accurate in predicting stall behavior than the Spalart-Allmaras, shear stress transport, and original k-(v^2)-{omega} models for the full high-lift configuration. The relative errors in the predicted maximum lift coefficients are within 3% of the experimental data.
A new approach to turbulence simulation, based on a combination of large-eddy simulation (LES) for the whole flow and an array of non-space-filling quasi-direct numerical simulations (QDNS), which sample the response of near-wall turbulence to large-scale forcing, is proposed and evaluated. The technique overcomes some of the cost limitations of turbulence simulation, since the main flow is treated with a coarse-grid LES, with the equivalent of wall functions supplied by the near-wall sampled QDNS. Two cases are tested, at friction Reynolds number Re$_tau$=4200 and 20,000. The total grid node count for the first case is less than half a million and less than two million for the second case, with the calculations only requiring a desktop computer. A good agreement with published DNS is found at Re$_tau$=4200, both in terms of the mean velocity profile and the streamwise velocity fluctuation statistics, which correctly show a substantial increase in near-wall turbulence levels due to a modulation of near-wall streaks by large-scale structures. The trend continues at Re$_tau$=20,000, in agreement with experiment, which represents one of the major achievements of the new approach. A number of detailed aspects of the model, including numerical resolution, LES-QDNS coupling strategy and sub-grid model are explored. A low level of grid sensitivity is demonstrated for both the QDNS and LES aspects. Since the method does not assume a law of the wall, it can in principle be applied to flows that are out of equilibrium.
In a recent paper, Liu, Zhu and Wu (2015, {it J. Fluid Mech.} {bf 784}: 304) present a force theory for a body in a two-dimensional, viscous, compressible and steady flow. In this companion paper we do the same for three-dimensional flow. Using the fundamental solution of the linearized Navier-Stokes equations, we improve the force formula for incompressible flow originally derived by Goldstein in 1931 and summarized by Milne-Thomson in 1968, both being far from complete, to its perfect final form, which is further proved to be universally true from subsonic to supersonic flows. We call this result the textit{unified force theorem}, which states that the forces are always determined by the vector circulation $pGamma_phi$ of longitudinal velocity and the scalar inflow $Q_psi$ of transverse velocity. Since this theorem is not directly observable either experimentally or computationally, a testable version is also derived, which, however, holds only in the linear far field. We name this version the textit{testable unified force formula}. After that, a general principle to increase the lift-drag ratio is proposed.
In the theory of the Navier-Stokes equations, the viscous fluid in incompressible flow is modelled as a homogeneous and dense assemblage of constituent fluid particles with viscous stress proportional to rate of strain. The crucial concept of fluid flow is the velocity of the particle that is accelerated by the pressure and viscous interaction around it. In this paper, by virtue of the alternative constituent micro-finite element, we introduce a set of new intrinsic quantities, called the vortex fields, to characterise the relative orientation between elements and the feature of micro-eddies in the element, while the description of viscous interaction in fluid returns to the initial intuition that the interlayer friction is proportional to the slip strength. Such a framework enables us to reconstruct the dynamics theory of viscous fluid, in which the flowing fluid can be modelled as a finite covering of elements and consequently indicated by a space-time differential manifold that admits complex topological evolution.