From CFD-Wiki
(Difference between revisions)
|
|
| Line 1: |
Line 1: |
| - | As its name means, gradient-based methods need the gradient of objective functions to design variables. The evaluation of gradient can be achieved by '''finite difference method''', '''linearized method''' or '''adjoint method'''. Both finite difference method and linearized method has a time-cost proportional to the number of design variables and not suitable for design optimization with a large number of design variables. Apart from that, finite difference method has a notorious disadvantage of subtraction cancellation and is not recommended for practical design application.
| |
| | | | |
| - | Suppose a cost function <math>J</math> is defined as follows,
| |
| - |
| |
| - | <math>J=J(U,\alpha)</math>
| |
| - |
| |
| - | where <math>U</math> and <math>\alpha</math> are the flow variable vector and the design variable vector respectively. <math>U</math> and <math>\alpha</math> are implicitly related through the flow equation, which is represented by a residual function driven to zero.
| |
| - |
| |
| - | <math>{R}(U(\alpha),\alpha)=0</math>
| |
| - |
| |
| - | The sensitivity of the cost function <math>J</math> with respect to the design variables <math>\alpha</math>, that is <math>\frac{D J }{D \alpha_{i}}</math>, is needed for design purpose. The following is three main methods to obtain this sensitivity.
| |
| - |
| |
| - | Finite difference method is the most straightforward approach, where the sensitivity is calculated through finite difference, using different cost function values corresponding to different design variable input
| |
| - |
| |
| - | <math>\frac{DJ}{D\alpha_{i}}=\frac{J(\alpha_{i}+\delta \alpha)-J(\alpha_{i})}{\delta \alpha}</math>
| |
| - |
| |
| - | The defect with this approach is that first of all, not efficient, because the computational cost is linearly proportional to the number of design variables which is practically too expensive. Second of all, the interval <math>\delta\alpha</math> is difficult to determine due to the concern of accuracy and machine error.
| |
| - |
| |
| - |
| |
| - | Linearized method:
| |
| - |
| |
| - | [[Adjoint method]]:
| |
| - |
| |
| - | ''Reference:Time-stepping for adjoint CFD codes from automatic differentiation, 2010, Faidon Christakopoulos, Dominic Jones and Jens D. Mueller''
| |
Latest revision as of 01:33, 6 January 2012
