×

The effect of shocks on second order sensitivities for the quasi-one-dimensional Euler equations. (English) Zbl 1367.76032

Summary: The effect of discontinuity in the state variables on optimization problems is investigated on the quasi-one-dimensional Euler equations in the discrete level. A pressure minimization problem and a pressure matching problem are considered. We find that the objective functional can be smooth in the continuous level and yet be non-smooth in the discrete level as a result of the shock crossing grid points. Higher resolution can exacerbate that effect making grid refinement counter productive for the purpose of computing the discrete sensitivities. First- and second-order sensitivities, as well as the adjoint solution, are computed exactly at the shock and in its vicinity, and are compared to the continuous solution. It is shown that in the discrete level the first-order sensitivities contain a spike at the shock location that converges to a delta function with grid refinement, consistent with the continuous analysis. The numerical Hessian is computed, and its consistency with the analytical Hessian is discussed for different flow conditions. It is demonstrated that consistency is not guaranteed for shocked flows. We also study the different terms composing the Hessian and propose some stable approximation to the continuous Hessian.

MSC:

76L05 Shock waves and blast waves in fluid mechanics
76M12 Finite volume methods applied to problems in fluid mechanics
76N15 Gas dynamics (general theory)

Software:

HE-E1GODF; TAPENADE
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] J. Betts, W. Huffman, D. Young, An investigation of algorithm performance for aerodynamic design optimization, Boeing Technical Report BCSTECH-94-061, Boeing Computer Services, 1994.; J. Betts, W. Huffman, D. Young, An investigation of algorithm performance for aerodynamic design optimization, Boeing Technical Report BCSTECH-94-061, Boeing Computer Services, 1994.
[2] Sherman, L.; Taylor, A.; Green, L.; Newman, P., First- and second-order aerodynamic sensitivity derivatives via automatic differentiation with incremental iterative methods, J. Comput. Phys., 129, 307-331 (1996) · Zbl 0933.76070
[3] D. Ghate, M. Giles, Efficient Hessian calculation using automatic differentiation, in: AIAA Paper 2007-4059.; D. Ghate, M. Giles, Efficient Hessian calculation using automatic differentiation, in: AIAA Paper 2007-4059.
[4] Arian, E.; Ta’asan, S., Analysis of the Hessian for aerodynamic shape optimization: inviscid flow, Comput. Fluids, 28, 853-877 (1999) · Zbl 0969.76015
[5] Arian, E.; Vatsa, V., A preconditioning method for shape optimization governed by the Euler equations, Int. J. Comput. Fluid D, 12, 17-27 (1999) · Zbl 0962.76607
[6] Coleman, T.; Garbow, B.; More, J., Software for estimating sparse Hessian matrices, ACM Trans. Math. Software, 11, 363-377 (1985) · Zbl 0584.65010
[7] W. Jou, W. Huffman, D. Young, R. Melvin, M. Bieterman, C. Hilmes, F. Johnson, Practical considerations in aerodynamic design optimization, in: AIAA Paper 1995-1730.; W. Jou, W. Huffman, D. Young, R. Melvin, M. Bieterman, C. Hilmes, F. Johnson, Practical considerations in aerodynamic design optimization, in: AIAA Paper 1995-1730.
[8] Papadimitriou, D.; Giannakoglou, K., The continuous direct-adjoint approach for second order sensitivities in viscous aerodynamic inverse design problems, Comput. Fluids, 38, 1539-1548 (2009) · Zbl 1242.76300
[9] Arian, E.; Iollo, A., Challenges in aerodynamic optimization, (Hafez, M.; Oshima, K.; Kwak, D., Computational Fluid Dynamics Review 2010 (2010), World Scientific), 447-467
[10] Zervogiannis, T.; Papadimitrou, D.; Giannakoglou, K., Total pressure losses minimization in turbomachinery cascades using the exact Hessian, Comput. Methods Appl. Mech. Eng., 199, 2697-2708 (2010) · Zbl 1231.76247
[11] T. Matsuzawa, M. Hafez, Treatment of shock waves in design optimization via adjoint equation approach, in: AIAA Paper 1998-2537.; T. Matsuzawa, M. Hafez, Treatment of shock waves in design optimization via adjoint equation approach, in: AIAA Paper 1998-2537.
[12] Bardos, C.; Pironneau, O., A formalism for the differentiation of conservation laws, C.R. Math., 335, 839-845 (2002) · Zbl 1020.35052
[13] Bardos, C.; Pironneau, O., Data assimilation for conservation laws, Methods Appl. Anal., 12, 103-134 (2005) · Zbl 1114.35122
[14] Arian, E.; Iollo, A., Analytic Hessian derivation for the quasi-one-dimensional Euler equations, J. Comput. Phys., 228, 476-490 (2009) · Zbl 1162.76049
[15] Giles, M.; Pierce, N., Analytic adjoint solutions for the quasi-one-dimensional Euler equations, J. Fluid Mech., 426 (2001) · Zbl 0967.76079
[16] Toro, E., Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction (2009), Springer · Zbl 1227.76006
[17] Saad, Y.; Schultz, H., Gmres: A general minimal residual algorithm for solving nonsymmetric linear systems, SIAM J. Sci. Stat. Comput., 7, 856-869 (1986) · Zbl 0599.65018
[18] L. Hascoaet, V. Pascual, Tapenade 2.1 user’s guide, INRIA, 2004.; L. Hascoaet, V. Pascual, Tapenade 2.1 user’s guide, INRIA, 2004.
[19] A. Iollo, M. Salas, S. Ta’asan, Shape optimization governed by the Euler equations using an adjoint method, Technical Report 93-78, ICASE, 1993.; A. Iollo, M. Salas, S. Ta’asan, Shape optimization governed by the Euler equations using an adjoint method, Technical Report 93-78, ICASE, 1993.
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.