1 / 28

Solving Inverse Design Problems Involving Radiant Enclosures through Gradient-Based Optimization

Solving Inverse Design Problems Involving Radiant Enclosures through Gradient-Based Optimization. K. J. Daun, J. R. Howell, and D. P. Morton Department of Mechanical Engineering The University of Texas at Austin. Radiant Enclosure Design Problems.

ami
Télécharger la présentation

Solving Inverse Design Problems Involving Radiant Enclosures through Gradient-Based Optimization

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Solving Inverse Design Problems Involving Radiant Enclosures through Gradient-Based Optimization K. J. Daun, J. R. Howell, and D. P. Morton Department of Mechanical Engineering The University of Texas at Austin

  2. Radiant Enclosure Design Problems • Radiant enclosure design problems are encountered in diverse industrial settings • Examples: • Paint drying in the automotive industry • Baking ovens in the food industry • RTP of semiconductor wafers • Optimal geometry of solar collector-concentrators

  3. General Form of the Design Problem Heater Surface qs = ? qs = ? • Enclosure contains a heater surface and a design surface • Goal is to find the enclosure geometry and heater settings that produces the desired qs and T over the design surface • In order for the problem to be well-posed, only one BC can be specified on each surface! • Can be solved using forward, inverse, or optimization design methodologies qs = ? qs = ? qs = ? qs = ? Design Surface qsDS, TDS are known

  4. Forward Design Methodology • Either TDS or qsDS is specified over each surface • Assume TDSis specified over the design surface • The designer • Guesses heater settings and enclosure geometry • Calculates qsover the design surface • Adjusts design and repeats analysis if necessary • Requires many iterations, and the final solution quality is limited

  5. Inverse Design Methodology • Both TDSand qsDS are specified over the design surface, while the heater settings are unspecified • The inverse design problem is in its explicit form, which is ill-posed • The resulting set of ill-conditioned equations, Ax = b, must be solved using regularization methods • Finds a solution in few iterations, but it is often non-physical • This method is only applicable to certain types of enclosure design problems

  6. Optimization Design Methodology • TDS is specified over the design surface • Define an objective function, F(F), and design parameters, F, that control enclosure geometry and heater settings • Optimal design is found by minimizing F(F), usually through gradient-based minimization • Much more efficient than forward “trial-and-error” method • Easy to implement design constraints by restricting the domain of F

  7. Gradient-Based Minimization • Objective function local minimum is found iteratively • At the kth step: • Check if Fk = F*, usually by checking if |F(Fk)| < ecrit • Choose search direction, pk • Choose step size, ak,usually by minimizing f(ak) = F(Fk+akpk) • Take a step, Fk+1 = Fk + akpk

  8. Gradient-Based Minimization • Search direction is chosen based on 1st and 2nd-order curvature at F(Fk): • Search direction choices:

  9. Classes of Enclosure Design Problems Heater Surface q1= ? q2= ? q3= ? q4= ? q5= ? q1 = ? q2 = ? q3 = ? q4 = ? Design Surface T = Tspec, q = qspec Design Surface T = Tspec, qs = qspec Heater Surface q1= ? q4(t)= ? q1(t)= ? q2(t)= ? q3(t)= ? q5(t)= ? q2= ? Heater Surface q3= ? q4= ? q5= ? V Design Surface T = Tspec, q = qspec Design Surface T(t) = Tspec, q(t) = qspec

  10. Inverse Boundary Condition Design Heater Surface qs= ? qs= ? qs= ? qs= ? qs= ? Design Surface T = Tspec, qs = qspec • Goal is to find heat flux distribution over the heater surface • Enclosure geometry is fixed • Solved using both inverse and optimization design methodologies • Inverse design methodologies by Oguma and Howell (1995), Harutunian et al. (1995) • Optimization design methodology by Daun et al. (2003)

  11. Governing Radiosity Equations = + Radiosity Emitted Radiation Reflected Incident Radiation = + Incoming Energy Outgoing Energy

  12. Solution of Governing Equations • Calculating F(F) and 2F(F) requires 1st- and 2nd-order heat flux sensitivity, e.g. • Radiosity sensitivities are solved by post-processing the radiosity solution • Solve Ax = b, xi = qoi(F) • Solve Ax = b, xi = qoi(F)/Fp • Solve Ax = b, xi = 2qoi(F)/FpFq • Heat flux sensitivities are found from radiosity sensitivities

  13. Example Problem: Inverse BC Design F9 F10 F11 F12 F12 F11 F10 F9 F8 F8 F7 F7 Heater Surface e = 0.9 F6 F6 F5 F5 F4 F4 Design Surface qstarget = 2 W/m2 Eb = 1 W/m2 e = 0.5 F3 F3 F2 F2 F1 F1 qs = 0W/m2 qs = 0W/m2

  14. Example Problem: Inverse BC Design 12 sT4(u) [W/m2] 10 8 F9 F9 6 F1 F1 4 F4 F4 F3 F3 2 0 0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 u

  15. Example Problem: Inverse BC Design 2 sT4(u) [W/m2] (specified) 1.5 1 0.5 0 -0.5 -1 qs(u, F*) [W/m2] -1.5 -2 -2.5 -3 0.95 0.65 0.7 0.75 0.8 0.85 0.9 0.95 u

  16. Optimization of Transient Problems Heater Surface q5(t)= ? q4(t)= ? q1(t)= ? q2(t)= ? q3(t)= ? V Design Surface T(t) = Tspec(t), q(t) = qspec(t) • Solved using: • Gradient-based optimization (Fedorov et al., 1998) • Inverse design methodology (Ertürk et al., 2001) • Nonlinear controls with regularization (Gwak and Masada, 2002)

  17. Governing Equations d(u) + = - Stored Energy Incoming Energy Outgoing Energy • Enclosure surface discretized into N elements, and time domain discretized into Nt time steps • Results in a non-linear system of equations that must be solved at each time step

  18. Transient Optimization • Objective function is defined as • Heater settings are controlled by cubic splines, where t = t/tmax. • Objective function sensitivities found analytically, through direct differentiation

  19. Example Problem: Transient Optimization 1 2 3 4 5 6 7 8 Heater surface T1 = 300 K e1= 1 T2 = 1000 K e2 = 1 Design surface Refractory surface V t =0 t =1

  20. Example Problem: Transient Optimization 1 3 4 5 6 7 8 2 1 0.9 0.8 0.7 0.6 qs(F)/sT24 0.5 0.4 0.3 0.2 0.1 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 t

  21. Example Problem: Transient Optimization 1.2 x = 0 x = 1 1 Design Surface x = 1 0.8 Target T(F, t)/Ts 0.6 x = 0 0.4 0.2 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 t

  22. Geometric Optimization F2 F1 qs(F3) qs(F4) Heater Surface qs(F5) qs(F6) qs(F7) Design Surface T = Tspec, q = qspec • Both heater settings and geometry can be designed • Only optimization methods have been used • Deterministic optimization for diffuse-gray enclosures (infinitesimal-area method) • Stochastic optimization for enclosures with specular surfaces (Monte Carlo method)

  23. Enclosures Containing Specular Surfaces ` • The exchange factor, Fji, is estimated using Monte Carlo analysis: • The uncertainty in Fji induces a random error in qs(F) and F(F) Energy in from other sources Energy leaving by emission Energy in from surrounding surfaces – = = -

  24. Stochastic Optimization • Uses Kiefer-Wolfowitz scheme based on steepest-descent, used when an unbiased estimate of F(F) is unavailable • At the kth iteration, • Check if Fk = F* • Set • Set ak = a0/ka, 0 a  1 • Fk+1 = Fk + akpk • The gradient is estimated by central finite difference, with

  25. Example Problem: Geometric Optimization F2 Specular Surface qs3 = 0 W/m2 r3 = 1 F1 Specular Surface qs2 = 0 W/m2 r2 = 1 Heater Surface qs1 = 1 W/m2 e1 = 1 u Design Surface Eb4 = 0 W/m2, qstarget = -1 W/m2, e4 = 1

  26. Example Problem: Geometric Optimization F2 F1 F2 F1 Initial Enclosure Geometry Minimization Path Optimal Enclosure Geometry

  27. Example Problem: Geometric Optimization qs(u, F0) (initial) qs(u, F*) (optimal)

  28. Conclusions • Optimization is used to design many types of radiant enclosures • Solves the inverse design problem implicitly, through iteration • More efficient than the forward design method, and usually produces better solutions • More straightforward than inverse design (regularization) method, and easier to implement design constraints

More Related