An isogeometric numerical study of partially and fully implicit schemes for transient adjoint shape sensitivity analysis

Abstract

In structural design optimization involving transient responses, time integration scheme plays a crucial role in sensitivity analysis because it affects the accuracy and stability of transient analysis. In this work, the influence of time integration scheme is studied numerically for the adjoint shape sensitivity analysis of two benchmark transient heat conduction problems within the framework of isogeometric analysis. It is found that (i) the explicit approach ( \beta = 0) and semi-implicit approach with \beta<0.5 impose a strict stability condition of the transient analysis; (ii) the implicit approach ( \beta=1) and semi-implicit approach with \beta > 0.5 are generally preferred for their unconditional stability; and (iii) Crank–Nicolson type approach ( \beta=0.5) may induce a large error for large time-step sizes due to the oscillatory solutions. The numerical results also show that the time-step size does not have to be chosen to satisfy the critical conditions for all of the eigen-frequencies. It is recommended to use \beta \approx \mathrm{0.75} for unconditional stability, such that the oscillation condition is much less critical than the Crank–Nicolson scheme, and the accuracy is higher than a fully implicit approach.

Keywords

isogeometric shape optimization / design-dependent boundary condition / transient heat conduction / implicit time integration / adjoint method

References

[1] Li Q, Steven G P, Querin O M, Shape and topology design for heat conduction by evolutionary structural optimization. International Journal of Heat and Mass Transfer, 1999, 42(17): 3361–3371

[2] Xie G, Liu Y, Sunden B, Computational study and optimization of laminar heat transfer and pressure loss of double-layer microchannels for chip liquid cooling. Journal of Thermal Science and Engineering Applications, 2013, 5(1): 011004

[3] Sigmund O, Torquato S. Design of materials with extreme thermal expansion using a three-phase topology optimization method. Journal of the Mechanics and Physics of Solids, 1997, 45(6): 1037–1067

[4] Gao T, Zhang W, Zhu J, Topology optimization of heat conduction problem involving design-dependent heat load effect. Finite Elements in Analysis and Design, 2008, 44(14): 805–813

[5] Iga A, Nishiwaki S, Izui K, Topology optimization for thermal conductors considering design-dependent effects, including heat conduction and convection. International Journal of Heat and Mass Transfer, 2009, 52(11–12): 2721–2732

[6] Yaji K, Yamada T, Kubo S, A topology optimization method for a coupled thermal-fluid problem using level set boundary expressions. International Journal of Heat and Mass Transfer, 2015, 81: 878–888

[7] Xia Q, Xia L, Shi T. Topology optimization of thermal actuator and its support using the level set based multiple-type boundary method and sensitivity analysis based on constrained variational principle. Structural and Multidisciplinary Optimization, 2018, 57(3): 1317–1327

[8] Choi K K, Kim N H. Structural Sensitivity Analysis and Optimization 1: Linear Systems. New York: Springer, 2005

[9] Dems K, Rousselet B. Sensitivity analysis for transient heat conduction in a solid body-Part I: External boundary modification. Structural Optimization, 1999, 17(1): 36–45

[10] Dems K, Rousselet B. Sensitivity analysis for transient heat conduction in a solid body-Part II: Interface modification. Structural Optimization, 1999, 17(1): 46–54

[11] Haftka R T, Shore C P. Approximation Methods for Combined Thermal/Structural Design. NASA Technical Paper 1428. 1979

[12] Haftka R T. Techniques for thermal sensitivity analysis. International Journal for Numerical Methods in Engineering, 1981, 17(1): 71–80

[13] Greene W H, Haftka R T. Computational aspects of sensitivity calculations in transient structural analysis. Computers & Structures, 1989, 32(2): 433–443

[14] Gao Z Y, Grandhi R V. Sensitivity analysis and shape optimization for preform design in thermo-mechanical coupled analysis. International Journal for Numerical Methods in Engineering, 1999, 45(10): 1349–1373

[15] Haftka R T, Malkus D S. Calculation of sensitivity derivatives in thermal problems by finite differences. International Journal for Numerical Methods in Engineering, 1981, 17(12): 1811–1821

[16] Wang Z P, Turteltaub S, Abdalla M M. Shape optimization and optimal control for transient heat conduction problems using an isogeometric approach. Computers & Structures, 2017, 185: 59–74

[17] Michaleris P, Tortorelli D A, Vidal C A. Tangent operators and design sensitivity formulations for transient non-linear coupled problems with applications to elastoplasticity. International Journal for Numerical Methods in Engineering, 1994, 37(14): 2471–2499

[18] Tortorelli D A, Haber R B, Lu S C Y. Design sensitivity analysis for nonlinear thermal systems. Computer Methods in Applied Mechanics and Engineering, 1989, 77(1–2): 61–77

[19] Tortorelli D A, Haber R B. First-order design sensitivities for transient conduction problems by an adjoint method. International Journal for Numerical Methods in Engineering, 1989, 28(4): 733–752

[20] Wang Z P, Kumar D. On the numerical implementation of continuous adjoint sensitivity for transient heat conduction problems using an isogeometric approach. Structural and Multidisciplinary Optimization, 2017, 56(2): 487–500

[21] Kane J H, Kumar B L, Stabinsky M. Transient thermoelasticity and other body force effects in boundary element shape sensitivity analysis. International Journal for Numerical Methods in Engineering, 1991, 31(6): 1203–1230

[22] Jarny Y, Ozisik M N, Bardon J P. A general optimization method using adjoint equation for solving multidimensional inverse heat conduction. International Journal of Heat and Mass Transfer, 1991, 34(11): 2911–2919

[23] Kleiber M, Służalec A. Material derivative and control volume approaches to shape sensitivity analysis of nonlinear transient thermal problems. Structural Optimization, 1996, 11: 56–63

[24] Dorai G A, Tortorelli D A. Transient inverse heat conduction problem solutions via Newton’s method. International Journal of Heat and Mass Transfer, 1997, 40(17): 4115–4127

[25] Korycki R. Sensitivity analysis and shape optimization for transient heat conduction with radiation. International Journal of Heat and Mass Transfer, 2006, 49(13–14): 2033–2043

[26] Huang C H, Chaing M T. A transient three-dimensional inverse geometry problem in estimating the space and time-dependent irregular boundary shapes. International Journal of Heat and Mass Transfer, 2008, 51(21–22): 5238–5246

[27] Służalec A, Kleiber M. Shape optimization of thermo-diffusive systems. International Journal of Heat and Mass Transfer, 1992, 35(9): 2299–2304

[28] Gu Y X, Chen B S, Zhang H W, A sensitivity analysis method for linear and nonlinear transient heat conduction with precise time integration. Structural and Multidisciplinary Optimization, 2002, 24(1): 23–37

[29] Chen B, Tong L. Sensitivity analysis of heat conduction for functionally graded materials. Materials & Design, 2004, 25(8): 663–672

[30] Haftka R T, Grandhi R V. Structural shape optimization—A survey. Computer Methods in Applied Mechanics and Engineering, 1986, 57(1): 91–106

[31] van Keulen F, Haftka R T, Kim N H. Review of options for structural design sensitivity analysis, Part 1: Linear systems. Computer Methods in Applied Mechanics and Engineering, 2005, 194(30–33): 3213–3243

[32] Cho S, Ha S H. Isogeometric shape design optimization: Exact geometry and enhanced sensitivity. Structural and Multidisciplinary Optimization, 2009, 38(1): 53–70

[33] Qian X. Full analytical sensitivities in nurbs based isogeometric shape optimization. Computer Methods in Applied Mechanics and Engineering, 2010, 199(29–32): 2059–2071

[34] Nagy A P, Abdalla M M, Gürdal Z. Isogeometric sizing and shape optimisation of beam structures. Computer Methods in Applied Mechanics and Engineering, 2010, 199(17–20): 1216–1230

[35] Nagy A P, Abdalla M M, Gürdal Z. Isogeometric design of elastic arches for maximum fundamental frequency. Structural and Multidisciplinary Optimization, 2011, 43(1): 135–149

[36] Liu H, Yang D, Wang X, Smooth size design for the natural frequencies of curved Timoshenko beams using isogeometric analysis. Structural and Multidisciplinary Optimization, 2019, 59(4): 1143–1162

[37] Weeger O, Narayanan B, Dunn M L. Isogeometric shape optimization of nonlinear, curved 3D beams and beam structures. Computer Methods in Applied Mechanics and Engineering, 2019, 345: 26–51

[38] Nagy A P, IJsselmuiden S T, Abdalla M M. Isogeometric design of anisotropic shells: Optimal form and material distribution. Computer Methods in Applied Mechanics and Engineering, 2013, 264: 145–162

[39] Hirschler T, Bouclier R, Duval A, Isogeometric sizing and shape optimization of thin structures with a solid-shell approach. Structural and Multidisciplinary Optimization, 2019, 59(3): 767–785

[40] Lian H, Kerfriden P, Bordas S. Implementation of regularized isogeometric boundary element methods for gradient-based shape optimization in two-dimensional linear elasticity. International Journal for Numerical Methods in Engineering, 2016, 106(12): 972–1017

[41] Lian H, Kerfriden P, Bordas S. Shape optimization directly from CAD: An isogeometric boundary element approach using T-splines. Computer Methods in Applied Mechanics and Engineering, 2017, 317: 1–41

[42] Wang C, Xia S, Wang X, Isogeometric shape optimization on triangulations. Computer Methods in Applied Mechanics and Engineering, 2018, 331: 585–622

[43] Wang Z P, Poh L H, Dirrenberger J, Isogeometric shape optimization of smoothed petal auxetic structures via computational periodic homogenization. Computer Methods in Applied Mechanics and Engineering, 2017, 323: 250–271

[44] Wang Z P, Poh L H. Optimal form and size characterization of planar isotropic petal-shaped auxetics with tunable effective properties using IGA. Composite Structures, 2018, 201: 486–502

[45] Kumar D, Wang Z P, Poh L H, Isogeometric shape optimization of smoothed petal auxetics with prescribed nonlinear deformation. Computer Methods in Applied Mechanics and Engineering, 2019, 356: 16–43

[46] Wang Y, Benson D J. Geometrically constrained isogeometric parameterized level-set based topology optimization via trimmed elements. Frontiers of Mechanical Engineering, 2016, 11(4): 328–343

[47] Wang Y, Xu H, Pasini D. Multiscale isogeometric topology optimization for lattice materials. Computer Methods in Applied Mechanics and Engineering, 2017, 316: 568–585

[48] Xie X, Wang S, Xu M, A new isogeometric topology optimization using moving morphable components based on R-functions and collocation schemes. Computer Methods in Applied Mechanics and Engineering, 2018, 339: 61–90

[49] Lieu Q X, Lee J. Multiresolution topology optimization using isogeometric analysis. International Journal for Numerical Methods in Engineering, 2017, 112(13): 2025–2047

[50] Hou W, Gai Y, Zhu X, Explicit isogeometric topology optimization using moving morphable components. Computer Methods in Applied Mechanics and Engineering, 2017, 326: 694–712

[51] Liu H, Yang D, Hao P, Isogeometric analysis based topology optimization design with global stress constraint. Computer Methods in Applied Mechanics and Engineering, 2018, 342: 625–652

[52] Hao P, Yuan X, Liu C, An integrated framework of exact modeling, isogeometric analysis and optimization for variable-stiffness composite panels. Computer Methods in Applied Mechanics and Engineering, 2018, 339: 205–238

[53] Guo Y, Ruess M. Nitsche’s method for a coupling of isogeometric thin shells and blended shell structures. Computer Methods in Applied Mechanics and Engineering, 2015, 284: 881–905

[54] Cai S Y, Zhang W H, Zhu J, Stress constrained shape and topology optimization with fixed mesh: A B-spline finite cell method combined with level set function. Computer Methods in Applied Mechanics and Engineering, 2014, 278: 361–387

[55] Zhang W, Zhao L, Gao T, Topology optimization with closed B-splines and Boolean operations. Computer Methods in Applied Mechanics and Engineering, 2017, 315: 652–670

[56] Wang Y, Wang Z P, Xia Z, Structural design optimization using isogeometric analysis: A comprehensive review. Computer Modeling in Engineering & Sciences, 2018, 117(3): 455–507

[57] Xia L, Xia Q, Huang X, Bi-directional evolutionary structural optimization on advanced structures and materials: A comprehensive review. Archives of Computational Methods in Engineering, 2018, 25(2): 437–478

[58] Meng L, Zhang W, Quan D, From topology optimization design to additive manufacturing: Today’s success and tomorrow’s roadmap. Archives of Computational Methods in Engineering, 2019 (in press)

[59] Kaminski W. Hyperbolic heat conduction equation for materials with a nonhomogeneous inner structure. Journal of Heat Transfer, 1990, 112(3): 555–560

[60] Xia L, Breitkopf P. Recent advances on topology optimization of multiscale nonlinear structures. Archives of Computational Methods in Engineering, 2017, 24(2): 227–249

[61] Wang Z P, Turteltaub S. Isogeometric shape optimization for quasi-static processes. International Journal for Numerical Methods in Engineering, 2015, 104(5): 347–371

[62] Xia Q, Shi T, Liu S, A level set solution to the stress-based structural shape and topology optimization. Computers & Structures, 2012, 90: 55–64

[63] Xia Q, Shi T, Xia L. Stable hole nucleation in level set based topology optimization by using the material removal scheme of BESO. Computer Methods in Applied Mechanics and Engineering, 2019, 343: 438–452

[64] Reddy J N, Gartling D K. The Finite Element Method in Heat Transfer and Fluid Dynamics. Boca Raton: CRC Press, 2001

[65] Bergheau J M, Fortunier R. Finite Element Simulation of Heat Transfer. Hoboken: John Wiley & Sons, 2013

[66] Carter W C. Lecture Notes on Mathematics for Materials Science and Engineers. MIT 3.016, 2012

[67] Ho Lee D, Man Kwak B. Shape sensitivity and optimization for transient heat diffusion problems using the BEM. International Journal of Numerical Methods for Heat & Fluid Flow, 1995, 5(4): 313–326