×

Nonintrusive continuum sensitivity analysis for fluid applications. (English) Zbl 1453.76175

Summary: Continuum Sensitivity Analysis (CSA) provides an analytic method of computing derivatives for structures, fluids and fluid-structure-interaction problems with respect to shape or value parameters. CSA does not require mesh sensitivity to calculate local shape sensitivity, thereby contributing to its computational efficiency. CSA involves solving linear sensitivity equations with the corresponding sensitivity boundary conditions. Spatial gradients of flow variables are required to construct the sensitivity boundary conditions, the accuracy of which significantly affects the flow shape derivatives. A method called Spatial Gradient Reconstruction (SGR) is used to accurately compute the flow spatial gradients that are further used to set the sensitivity boundary conditions. Here we present a nonintrusive CSA formulation for calculating the material derivatives of flow variables with respect to shape design parameters. Nonintrusive CSA means CSA without any modification to the black-box program used for the CFD flow analysis. Moreover, nonintrusive CSA is flow solver agnostic, in that, knowledge of the solver’s discretization is not needed. An example of two-dimensional flow over a NACA0012 airfoil demonstrates the application of nonintrusive CSA to steady flow problems involving Euler (compressible inviscid) equations over unstructured grids with finite volume spatial discretization. Accuracy was studied by investigating convergence for a family of six high-quality grids, ranging from four thousand to four million cells. Factors such as the accuracy of the flow spatial gradients, flow sensitivity boundary conditions, weak implementation of the boundary conditions in the finite volume framework, and the use of an approximate flux Jacobian matrix, all of which affect the accuracy of the material derivatives, are discussed. A novelty of this work is that the sensitivity analysis is done nonintrusively using FLUENT and SU2 software to establish the use of black-box codes for obtaining flow derivatives using the CSA approach.

MSC:

76M30 Variational methods applied to problems in fluid mechanics
76M12 Finite volume methods applied to problems in fluid mechanics
76N15 Gas dynamics (general theory)

Software:

FLUENT; TFad; SU2
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] AIAA, 4th aiaa cfd drag prediction workshop. San Antonio, TX, June 2009.
[2] Albring, T.; Sagebaum, M.; Gauger, N. R., A Consistent and Robust Discrete Adjoint Solver for the SU2 Framework - Validation and Application, 77-86 (2016), Springer International Publishing: Springer International Publishing Cham
[3] Anderson, W. K.; Bonhaus, D. L., An implicit upwind algorithm for computing turbulent flows on unstructured grids, Comput. Fluids, 23, 1, 1-21 (1994) · Zbl 0806.76053
[4] Becker, K.; Ashcroft, G., A comparative study of gradient reconstruction methods for unstructured meshes with application to turbomachinery flows, (52nd AIAA Aerospace Sciences Meeting, National Harbor, MD, USA (January 2014))
[5] Borggaard, J.; Burns, J., A Sensitivity Equation Approach to Shape Optimization in Fluid Flows (1994), Langley Research Center, Nasa contractor report 191598
[6] Borggaard, J.; Burns, J., A pde sensitivity equation method for optimal aerodynamic design, J. Comput. Phys., 136, 366-384 (1997) · Zbl 0903.76064
[7] Charlot, L.; Etienne, S.; Pelletier, D., A continuous lagrangian sensitivity equation method for incompressible flow, J. Comput. Phys., 231, 5989-6011 (2012) · Zbl 1277.76022
[8] Chen, P. C.; Zhang, Z., Store ejection loads analysis using an euler solver, (54th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference (April 2013), AIAA: AIAA Boston, MA)
[9] Choi, K. K.; Kim, N.-H., Structural Sensitivity Analysis and Optimization (2005), Springer Science+Business media: Springer Science+Business media New York
[10] Choudhary, A.; Roy, C. J., Verification and Validation for Multiphase Flows, 1-37 (2018), Springer Singapore: Springer Singapore Singapore
[11] Cross, D. M.; Canfield, R. A., Continuum shape sensitivity with spatial gradient reconstruction of nonlinear aeroelastic gust response, (14th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference, No. AIAA 2012-55597, September, 17-19 2012, Indianapolis, Indiana (2012))
[12] Cross, D. M.; Canfield, R. A., Local continuum shape sensitivity with spatial gradient reconstruction, Struct. Multidiscip. Optim., 50, 6, 975-1000 (December 2014)
[13] Cross, D. M.; Canfield, R. A., Convergence study of local continuum sensitivity method using spatial gradient reconstruction, AIAA J. (2016)
[14] Diskin, B.; Thomas, J. L., Notes on accuracy of finite-volume discretization schemes on irregular grids, Appl. Numer. Math., 60, 3, 224-226 (Mar. 2010)
[15] Duvigneau, R.; Pelletier, D., On accurate boundary conditions for a shape sensitivity equation method, Int. J. Numer. Methods Fluids, 50, 2, 147-164 (2006) · Zbl 1162.76360
[16] Gay, D. M., Semiautomatic differentiation for efficient gradient computations, (Bücker, H. M.; Corliss, G.; Hovland, P.; Naumann, U.; Norris, B., Automatic Differentiation: Applications, Theory, and Implementations. Automatic Differentiation: Applications, Theory, and Implementations, Lecture Notes in Computational Science and Engineering. (2005), Springer), 147-158 · Zbl 1270.65012
[17] Gobal, K.; Grandhi, R. V.; Kolonay, R. M., Continuum sensitivity analysis for structural shape design variables using finite-volume method, AIAA J., 53, 2, 347-355 (2015)
[18] Godfrey, A. G.; Cliff, E. M., Sensitivity equations for turbulent flows, (39th AIAA Aerospace Sciences Meeting & Exhibit. No. 2001-1060. Reno, NV (January 8-11 2001))
[19] Hicks, R. M.; Henne, P. A., Wing design by numerical optimization, J. Aircr., 15, 7, 407-412 (1978)
[20] Hirsch, C., Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows (1990), John Wiley and Sons · Zbl 0742.76001
[21] Hogan, R. J., Fast reverse-mode automatic differentiation using expression templates in c++, ACM Trans. Math. Softw., 40, 4, 26:1-26:16 (Jul. 2014)
[22] Jameson, A., Solution of the euler equations for two dimensional transonic flow by a multigrid method, Appl. Math. Comput., 13, 3, 327-355 (1983) · Zbl 0545.76065
[23] Kulkarni, M. D., Continuum sensitivity analysis using boundary velocity formulation for shape derivatives (2016), Ph.D. thesis
[24] Kulkarni, M. D.; Canfield, R. A.; Patil, M. J., Nonintrusive continuum sensitivity analysis for aerodynamic shape optimization, (15th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference. No. AIAA 2014-2043 (June 2014), AIAA: AIAA Atlanta, GA)
[25] Kulkarni, M. D.; Canfield, R. A.; Patil, M. J., Nonintrusive continuum sensitivity analysis for aerodynamic shape optimization, (16th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference, No. AIAA 2014-3237 (June 2015), AIAA: AIAA Dallas, TX)
[26] Kulkarni, M. D.; Canfield, R. A.; Patil, M. J.; Alyahak, E. J., Integration of geometric sensitivity and spatial gradient reconstruction for aeroelastic shape optimization, (10th AIAA Multidisciplinary Design Optimization Conference. No. AIAA 2014-0470 (January 2014), AIAA: AIAA National Harbor, MD)
[27] Kulkarni, M. D.; Cross, D. M.; Canfield, R. A., Discrete adjoint formulation for continuum sensitivity analysis, AIAA J., 54, 2, 758-766 (2016)
[28] Liu, S.; Canfield, R. A., Boundary velocity method for continuum shape sensitivity of nonlinear fluid structure interaction problems, J. Fluids Struct., 40, 284-301 (2013)
[29] Liu, S.; Canfield, R. A., Equivalence of continuum and discrete analytic sensitivity methods for nonlinear differential equations, Struct. Multidiscip. Optim., 48, 6, 1173-1188 (2013)
[30] Liu, S.; Canfield, R. A., Two forms of continuum shape sensitivity method for fluid structure interaction problems, J. Fluids Struct., 62, 46-64 (2016)
[31] Palacios, F.; Colonno, M. R.; Aranake, A. C.; Campos, A.; Copeland, S. R.; Economon, T. D.; Lonkar, A. K.; Lukaczyk, T. W.; Taylor, T. W.R.; Alonso, J. J., Stanford university unstructured (su2): an open-source integrated computational environment for multi-physics simulation and design, (51st AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition. No. AIAA 2013-0287, Grapevine, Texas (January, 7-10 2013))
[32] Pelletier, D.; Hay, A.; Etienne, S.; Borggaard, J., The sensitivity equation method in fluid mechanics, Eur. J. Comput. Mech./Rev. Eur. Méc. Numér., 17, 1-2, 31-61 (2008) · Zbl 1292.76037
[33] Thomas, J. L.; Diskin, B.; Rumsey, C. L., Towards verification of unstructured-grid solvers, AIAA J., 46, 12, 3070-3079 (2008)
[34] Vassberg, J. C.; Jameson, A., In pursuit of grid convergence for two-dimensional euler solutions, J. Aircr., 47, 4, 1152-1166 (August 2010)
[35] Villa, A., Convergence of weakly imposed boundary conditions: the one-dimensional hyperbolic case, SIAM J. Sci. Comput., 31, 4, 3116-3127 (2009) · Zbl 1242.65183
[36] Villa, A.; Barbieri, L.; Malgesini, R., Ghost cell boundary conditions for the euler equations and their relationships with feedback control, Commun. Appl. Ind. Math., 3, 1 (2012) · Zbl 1329.65199
[37] Zhang, Z.; Liuy, F.; Schusterz, D. M., An efficient euler method on non-moving cartesian grids with boundary-layer correction for wing flutter simulations, (Aerospace Sciences Meeting, Reno, Nevada (2006))
[38] Zienkiewicz, O. C.; Zhu, J. Z., The superconvergent patch recovery (spr) and adaptive finite element refinement, Comput. Methods Appl. Mech. Eng., 101, 1, 207-224 (1992) · Zbl 0779.73078
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.