No Arabic abstract
Adjoint method is widely used in aerodynamic design because only once solution of flow field is required for adjoint method to obtain the gradients of all design variables. However, the calculation cost of adjoint vector is approximately equal to that of flow computation. In order to accelerate the solution of adjoint vector and improve the adjoint-based optimization efficiency, machine learning for adjoint vector modeling is presented. Deep neural network (DNN) is employed to construct the mapping between the adjoint vector and the local flow variables. DNN can efficiently predict adjoint vector and its generalization is examined by a transonic drag reduction about NACA0012 airfoil. The results indicate that with negligible calculation cost of the adjoint vector, the proposed DNN-based adjoint method can achieve the same optimization results as the traditional adjoint method.
This article presents an original methodology for the prediction of steady turbulent aerodynamic fields. Due to the important computational cost of high-fidelity aerodynamic simulations, a surrogate model is employed to cope with the significant variations of several inflow conditions. Specifically, the Local Decomposition Method presented in this paper has been derived to capture nonlinear behaviors resulting from the presence of continuous and discontinuous signals. A combination of unsupervised and supervised learning algorithms is coupled with a physical criterion. It decomposes automatically the input parameter space, from a limited number of high-fidelity simulations, into subspaces. These latter correspond to different flow regimes. A measure of entropy identifies the subspace with the expected strongest non-linear behavior allowing to perform an active resampling on this low-dimensional structure. Local reduced-order models are built on each subspace using Proper Orthogonal Decomposition coupled with a multivariate interpolation tool. The methodology is assessed on the turbulent two-dimensional flow around the RAE2822 transonic airfoil. It exhibits a significant improvement in term of prediction accuracy for the Local Decomposition Method compared with the classical method of surrogate modeling for cases with different flow regimes.
We extend the impulse theory for unsteady aerodynamics, from its classic global form to finite-domain formulation then to minimum-domain form, and from incompressible to compressible flows. For incompressible flow, the minimum-domain impulse theory raises the finding of Li and Lu (J. Fluid Mech., 712: 598-613, 2012) to a theorem: The entire force with discrete wake is completely determined by only the time rate of impulse of those vortical structures still connecting to the body, along with the Lamb-vector integral thereof that captures the contribution of all the rest disconnected vortical structures. For compressible flow, we find that the global form in terms of the curl of momentum, obtained by Huang (Unsteady Vortical Aerodynamics. Shanghai Jiaotong Univ. Press, 1994), can be generalized to having arbitrary finite domain, but the formula is cumbersome and in general the curl of momentum no longer has discrete structure and hence no minimum-domain theory exists. Nevertheless, as the measure of transverse process only, the unsteady field of vorticity may still have discrete wake. This leads to a minimum-domain compressible vorticity-moment theory in terms of density-weighted vorticity (but it is beyond the classic concept of impulse). These new findings and applications have been confirmed by our numerical experiments. The results not only open an avenue to combine the theory with computation-experiment in wide applications, but also reveals a physical truth that it is no longer necessary to account for all wake vortical structures in computing the force and moment.
Here, through a systematic methodology and the use of high performance computing, we calculate the optimum shape for a wave energy converter under the action of incident waves of (i) monochromatic unidirectional, (ii) monochromatic directional, (iii) polychromatic unidirectional and (iv) polychromatic directional (with both directional symmetry and asymmetry). As a benchmark for our study, without loss of generality, we consider a submerged planar pressure differential wave energy converter, and use Genetic Algorithm to search through a wide range of shapes. A new parametric description of absorber shape based on Fourier decomposition of geometrical shapes is introduced, and for each shape hydrodynamic coefficients are calculated, optimum power take-o? parameters are obtained, and overall efficiency is determined. We show that an optimum geometry of the absorber plate can absorb a significantly higher energy (in some cases a few times higher) when compared to a circular shape of the same area. Specifically, for a unidirectional incident wave, the optimum shape, as expected, is obtained to be the most elongated shape. For directional incident waves, a butterfly-shape is the optimum geometry whose details depend on not only the amplitude and direction of incident wave components, but also the relative phases of those components. For the latter effect, we find an optimally averaged profile through a statistical analysis. Keywords: Wave energy conversion, Shape optimization
A sensitive porosity adjoint method (SPAM) for optimizing the topology of fluid machines has been proposed. A sensitivity function with respect to the porosity has been developed. In the first step of the optimization process, porous media are introduced into the flow regime according to the sensitivity function. Then the optimized porous media are transformed to solid walls. The turbulent flow in porous media is accounted for by a modified eddy-viscosity based turbulence model. Its influence on the adjoint equations is nevertheless neglected, which refers to the so called frozen turbulence assumption. A test case of application in terms of the turbulent rough wall channel flow shows that a considerable reduction of the objective function can be obtained by this method. The transformation from porous media to solid walls may have important effect on the optimization results.
Adjoint-based sensitivity analysis methods are powerful tools for engineers who use flow simulations for design. However, the conventional adjoint method breaks down for scale-resolving simulations like large-eddy simulation (LES) or direct numerical simulation (DNS), which exhibit the chaotic dynamics inherent in turbulent flows. Sensitivity analysis based on least-squares shadowing (LSS) avoids the issues encountered by conventional methods, but has a high computational cost. The following report outlines a new, more computationally efficient formulation of LSS, non-intrusive LSS, and estimates its cost for several canonical flows using Lyapunov analysis.