• Journal of the Korean Society for Precision Engineering
• /
• v.20 no.8
• /
• pp.175-184
• /
• 2003

A method of the shape design sensitivity analysis based on the boundary integral equation formulation is presented for two-dimensional inhomogeneous thermal conducting solids with multiple domains. Shape variation of the external and interface boundary is considered. A sensitivity formula of a general performance functional is derived by taking the material derivative to the boundary integral identity and by introducing an adjoint system. In numerical analysis, state variables of the primal and adjoint systems are solved by the boundary element method using quadratic elements. Two numerical examples of a compound cylinder and a thermal diffuser are taken to show implementation of the shape design sensitivity analysis. Accuracy of the present method is verified by comparing analyzed sensitivities with those by the finite difference. As application to the shape optimization, an optimal shape of the thermal diffuser is found by incorporating the sensitivity analysis algorithm in an optimization program.

• Proceedings of the Korea Electromagnetic Engineering Society Conference
• /
• 2005.11a
• /
• pp.169-174
• /
• 2005

In this paper, we proposed a new algorithm of the inverse scattering for the reconstruction of unknown dielectric scatterers using the finite-difference time-domain method and the design sensitivity analysis. We introduced the design sensitivity analysis based on the gradient for the fast convergence of the reconstruction. By introducing the adjoint variable method for the efficient calculation, we derived the adjoint variable equation. As an optimal algorithm we used the steepest descent method and reconstructed the dielectric targets using the iterative estimation. To verify our algorithm we will show the numerical examples for the two-dimensional $TM^2$ cases.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1053-1057
• /
• 2004

This paper presents the adjoint variable design sensitivity analysis for thermal systems considering both conduction and convection heat transfer. Both nodal temperature and total heat flow are considered to be objective functions and design sensitivity formulas are derived for each case. For the case of convection heat transfer, the adjoint analysis is carefully proceeded to obtain a precise result. A topology optimization example is examined for a simple planar square plate in order to design a heat exchanger as verification.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.04a
• /
• pp.127-134
• /
• 2004

In this paper, using an adjoint variable method, we develop a design sensitivity analysis (DSA) method applicable to heat conduction problems in steady state. Also, a topology design optimization method is developed using the developed DSA method. Design sensitivity expressions with respect to the thermal conductivity are derived. Since the already factorized system matrix is utilized to obtain the adjoint solution, the cost for the sensitivity computation is trivial. For the topology design optimization, the design variables are parameterized into normalized bulk material densities. The objective function and constraint are the thermal compliance of structures and allowable material volume, respectively. Through several numerical examples, the developed DSA method is verified to yield very accurate sensitivity results compared with finite difference ones, requiring less than 0.3% of CPU time far the finite differencing. Also, the topology optimization yields physical meaningful results.

• Journal of computational fluids engineering
• /
• v.11 no.1 s.32
• /
• pp.57-66
• /
• 2006

An optimal shape design approach is presented for a subsonic S-shaped intake using aerodynamic sensitivity analysis. Two-equation turbulence model is employed to capture strong counter vortices in the S-shaped duct more precisely. Sensitivity analysis is performed for the three-dimensional Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method For code validation, the result of the flow solver is compared with experiment data and other computational results of bench marking test. To study the influence oj turbulence models and grid refinement on the duct flow analysis, the results from several turbulence models are compared with one another and the minimum number of grid points, which can yield an accurate solution is investigated The adjoint variable code is validated by comparing the complex step derivative results. To realize a sufficient and flexible design space, NURBS equations are introduced as a geometric representation and a new grid modification technique, Least Square NURBS Grid Approximation is applied With the verified flow solver, the sensitivity analysis code and the geometric modification technique, the optimization of S-shaped intake is carried out and the enhancement of overall intake performance is achieved The designed S-shaped duct is tested in several off-design conditions to confirm the robustness of the current design approach. As a result, the capability and the efficiency of the present design tools are successfully demonstrated in three-dimensional highly turbulent internal flow design and off-design conditions.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.35 no.4
• /
• pp.435-447
• /
• 1999

It is main objective of this approach to present a method to analyse stochastic design sensitivity for problems of structural dynamics with randomness in design parameters. A combination of the adjoint variable approach and the second order perturbation method is used in the finite element approach. An alternative form of the constant functional that holds for all times is introduced to consider the time response of dynamic sensitivity. The terminal problem of the adjoint system is solved using equivalent homogeneous equations excited by initial velocities. The numerical procedures are shown to be much more efficient when based on the fold superposition method: the generalized co-ordinates are normalized and the correlated random variables are transformed to uncorrelated variables, whereas the secularities are eliminated by the fast Fourier transform of complex valued sequences. Numerical algorithms have been worked out and proved to be accurate and efficient : they can be readily adapted to fit into the existing finite element codes whose element derivative matrices can be explicitly generated. The numerical results of two cases -2 dimensional portal frame for the comparison with reference and 3-dimensional frame structure - for the deterministic sensitivity analysis are presented.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.4
• /
• pp.297-303
• /
• 2014

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that $C^1$-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.3 s.174
• /
• pp.645-652
• /
• 2000

A method for shape design sensitivity analysis is proposed utilizing commercial software ANSYS for thermal conduction problems. While the sensitivity formula is derived analytically by introduing adjoint variable concept, sensitivity calculation in practice as well as the primal and adjoint solution of thermal conduction is performed using the ANSYS very easily. Since the formula always takes boundary integral form, sensitivity evaluation in ANSYS requires a little more addition of post-processing routine which involves evaluation of boundary variable from the obtained solution. Though the BEM has been used as a better tool for this purpose, the present study shows it can also be calculated using any kind of analysis code such as ANSYS since the formula is based on analytic nature. Therefore the present study provides a new and efficient way of optimization which was not possible before using commercial software. The usefulness of the method is illustrated via a weight minimization problem of thermal diffuser.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.33 no.1
• /
• pp.1-10
• /
• 2005

An aerodynamic shape optimization technique has been developed for helicopter rotor blades in hover based on a continuous adjoint method on unstructured meshes. The Euler flow solver and the continuous adjoint sensitivity analysis were formulated on the rotating frame of reference for hovering rotor blades. In order to handle the repeated evaluation of the design cycle efficiently, the flow and adjoint solvers were parallelized using a domain decomposition strategy. A solution-adaptive mesh refinement technique was adopted for the accurate capturing of the tip vortex. Applications were made for the aerodynamic shape optimization of Caradonna-Tung rotor blades and UH60 rotor blades in hover. The results showed that the present method is an effective tool to determine optimum aerodynamic shapes of rotor blades requiring less torque while maintaining the desired thrust level.